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ABSTRACT 

This  thesis  examines  the  pool  boiling  heat  transfer 
characteristics  of  HCFC-124  and  HCFC-124/oil  mixtures  with  up 
to  10%  (by  weight)  miscible  alkylbenzene  oil.  One  smooth  and 
4  enhanced  tubes  were  tested:  a  19  and  26  low  integral -fin 
tube  (GEWA-K) ;  a  modified  finned  tube  (TURBO-B) ;  and  a  porous 
coated  tube  (HIGH  FLUX) .  The  tests  were  carried  out  using  the 
procedure  used  for  CFC-114  at  the  same  saturation  temperature 
of  2.2  °C.  This  allowed  for  direct  comparison  of  the  pool 
boiling  heat  transfer  characteristics  between  the  two 
refrigerants . 

The  smooth  and  GEWA-K  19  fin  per  inch  tube  performance  in 
pure  HCFC-124  and  HCFC-124/oil  mixtures  ranged  between  10  to 
50%  better  than  in  pure  CFC-114  and  CFC-114/oil  mixtures  for 
all  heat  fluxes.  The  HIGH  FLUX  and  TURBO-B  tubes  were  similar 
in  performance.  With  pure  HCFC-124,  the  finned  tubes  typically 
provided  enhancements  in  the  heat  transfer  coefficient  between 
2  and  3  times  that  of  a  smooth  tube.  The  HIGH  FLUX  and  TURBO-B 
surfaces  typically  provided  additional  enhancements  2  times 
that  of  the  finned  tubes.  With  the  addition  of  oil,  the  heat 
transfer  increased  from  the  smooth  and  finned  tubes  but 
decreased  from  the  HIGH  FLUX  and  TURBO-B  tubes.  The  HIGH  FLUX 
and  TURBO-B  tubes  therefore  exhibited  enhancements  less  than 
the  finned  tubes  at  high  oil  concentrations  and  high  heat 
fluxes . 
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NOMENCLATURE 

A  area 

Ab  tube  outside  surface  area  of  active  boiling  section 

Ac  cross- sectional  area  of  tube 

cp  specific  heat 

D  diameter 

Di  tube  inside  diameter 

DQ  outside  diameter  of  unenhanced  ends 

D1  diameter  at  the  position  of  the  thermocouple 

D2  diameter  of  test  section  to  the  base  of  fins 

g  gravitational  acceleration 

h  heat  -  transfer  coefficient 

hf  latent  heat  of  vaporization 

I  current 

Is  current  reading  by  AC  Current  Sensor 

k  thermal  conductivity  of  liquid  refrigerant 

kc  thermal  conductivity  of  copper 

L  active  boiling  tube  length 

Lu  non-boiling  length  of  the  test  tube 

Nu  Nusselt  number 

p  tube  outside  wall  perimeter 

Pr  Prandtl  Number 

Q  heat  -  transfer  rate  from  boiling  surface 

Qf  heat  -  transfer  rate  through  one  non-boiling  end 

Qh  heat  -  transfer  rate  from  cartridge  heater 

q  heat  flux 

Ra  Rayleigh  number 

T  temperature 

Ta  average  wall  temperature  at  the  thermocouple  location 

T  critical  temperature 


Tf  film  temperature 

Tn  temperature  of  the  respective  (n)  thermocouple 

location 

Tsat  saturation  temperature 

Two  outer  wall  temperature  of  the  boiling  test  tube 

V  voltage  across  the  cartridge  heater 

Vg  voltage  reading  by  AC -DC  true  RMS  converter 

a  thermal  diffusivity 

(3  volumetric  thermal  expansion  coefficient 

6  uncertainty  in  measurement  and  calibration 

8  superheat  (Two -Tsat) 

\i  dynamic  viscosity 

v  kinematic  viscosity 

p  density 

cr  surface  tension 
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I .  INTRODUCTION 

In  an  ongoing  effort  to  eliminate  the  use  of 
chlorof luorocarbons  (CFCs),  which  have  been  linked  to  the 
depletion  of  the  earth's  protective  ozone  layer,  the  United 
States  Navy  is  looking  into  using  recently  developed 
alternative  refrigerants.  There  are  2  principal  alternative 
refrigerants:  hydro f luorocarbons  (HFCs) ,  and 
hydrochlorof luorocarbons  (HCFCs) .  HFCs  contain  no  chlorine  and 
have  zero  ozone  depletion  potential  (ODP) .  HCFCs  do  contain 
chlorine,  but  with  the  attached  hydrogen  atom  to  the  CFC 
structure,  the  molecule  is  less  stable  and  the  ozone  damaging 
chlorine  is  released  in  the  lower  atmosphere  prior  to  reaching 
the  ozone  layer.  The  ODP  of  HCFC-124  is  about  2%, 
significantly  less  than  the  70%  ODP  for  CFC- 114.  Figure  1.1, 
provided  by  Baehr  [Ref.  1],  displays  the  advantages  of  HCFC- 
124  over  CFC- 114  with  respect  to  both  ozone  depletion 
potential  (ODP)  and  global  warming  potential  (GWP) . 

The  US  Navy  is  currently  using  CFC- 114  (R-114)  in  all  of 
its  larger  capacity  centrifugal  air  conditioning  units  (>100 
tons).  The  search  for  a  'drop-in'  replacement  has  identified 
HCFC-124  (R-124)  as  a  possible  alternative  to  R- 114 .  'Drop-in' 
here  means  operating  at  the  same  evaporating  temperature  as 
existing  systems.  Therefore  the  operating  parameters  of  R-124 
under  similar   conditions  as   present   R-114  chiller  systems 
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must  be  taken  into  account.  For  example,  the  shell  must  be 
capable  of  withstanding  higher  pressures  (11  psig  for  R-124  at 
2.2  °C  compared  to  0  psig  for  R-114  at  2.2  °C) .  The  shutdown 
system  pressure  is  also  a  significant  factor  when  the  system 
is  subjected  to  higher  temperatures  from  a  surrounding  engine 
room  or  perhaps  a  casualty  condition.  A  temperature  of  100  °F 
would  result  in  a  system  vapor  pressure  of  31  psig  for  R-114 
compared  to  68  psig  for  R-124.  If  the  heat  exchanger  shell 
can  handle  the  increased  system  pressure,  hardware 
modifications  should  be  limited  to  alterations  such  as 
compressor  impeller  size  and  changes  to  sealing  materials  to 
maintain  leak-tight  integrity. 

If  hardware  modifications  can  be  made,  then  the  other 
major  consideration  is  the  heat  transfer  characteristics  of  R- 
124  compared  to  those  of  R-114.  This  is  not  such  a  simple 
comparison  as  the  two  fluids  may  exhibit  varying  boiling 
characteristics  for  each  type  of  enhanced  surface.  In 
addition,  the  oil  used  for  each  fluid  is  different  since  each 
must  be  miscible  in  its  particular  refrigerant  and  this  may 
significantly  affect  heat  transfer  performance. 

To  evaluate  and  compare  the  heat  transfer  performance  of 

R-124  with  that  of  R-114,  the  following  thesis  objectives  were 

established: 

1.  Install  a  new  fiberglass -strengthened  evaporator  and 
condenser  section  into  existing  apparatus  to  withstand 
the  higher  saturation  pressure  associated  with  R-124. 


2.  Collect  pool  boiling  heat  transfer  data  of  pure  R-114  and 
R-114/oil  mixtures  over  a  range  of  increasing  and 
decreasing  heat  fluxes  for  the  newly  manufactured  tubes 
and  check  their  repeatability  with  existing  R-114  data 
using  the  newly  modified  apparatus. 

3.  Modify  existing  data  reduction  program  to  account  for  the 

thermophysical  properties  of  R-124. 

4.  Collect  pool  boiling  heat  transfer  data  of  pure  R-124  and 
R-124/oil  mixtures  over  a  range  of  increasing  and 
decreasing  heat  fluxes  using  a  smooth  tube  and  four 
enhanced  surface  tubes . 

5.  Study  effect  of  oil  on  boiling  hysteresis  for  a  smooth 
and  porous  coated  tube. 

6.  Compare  single  tube  heat  transfer  characteristics  of  R- 
124  with  previously  obtained  R-114  data  taken  by  Sugiyama 
[Ref .  2] . 


II.   POOL  BOILING  OF  REFRIGERANTS 

A.   POOL  BOILING  CURVE 

The  pool  boiling  curve,  as  provided  by  Bar- Cohen  [Ref .3] , 
is  shown  in  Figure  2.1.  The  applicable  regions  within  the 
curve  where  single  tube  pool  boiling  heat  transfer  tests  of  R- 
114  and  R-124  refrigerants  were  conducted  are  reviewed. 
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Figure  2.1   Pool  Boiling  Curve 


1.   Natural  Convection 

This  region  (a-b)  occurs  during  the  increasing  heat 
flux  phase  where  there  are  no  vapor  bubbles  being  produced  at 
the  tube  surface.  In  this  region,  heat  transfer  is  due  solely 
to  single-phase  free  convection  effects.  Several  correlations 
are  available  to  predict  heat  transfer  performance  in  this 
region,  including  the  Churchill  and  Chu  [Ref.  4]  correlation 
shown  below  in  equation  (1) . 

h=k/D0[.6+.3  87(g-/3-e-D03/u-a)  1/6/  [1+ ( 0  .  559/Pr)  9/16]  8/27]  2    (1) 


2 .   Incipient  Boiling 

The  transition  from  the  natural  convection  region  to 
the  nucleate  boiling  region  is  the  subject  of  numerous 
studies.  The  superheat  (Twall-Tsat)  at  which  bubble  inception 
(point  b)  occurs  is  referred  to  as  the  onset  of  nucleate 
boiling  (ONB)  or  incipient  boiling.  Considerable  reduction  in 
tube  wall  temperature  occurs  at  this  point,  resulting  in  a 
considerable  reduction  in  superheat  noted  as  ATex  in  Figure 
2.1.  This  'jump'  in  performance  only  occurs  during 
increasing  heat  flux  and  is  referred  to  as  the  thermal 
overshoot.  A  hysteresis  effect  occurs  when  decreasing  heat 
flux  data  (from  point  c  to  point  a)  are  compared  to  the 
initial  ascending  path  (a-b-c).  Several  factors  believed  to 
influence  the  thermal  overshoot  at  the  ONB  for  highly-wetting 


liquids  were  studied  by  Bar- Cohen  [Ref .  3] .  Probable 
influential  factors  included:  vibration,  surface  tension, 
enhanced  commercial  surfaces,  pressure,  time  between  test 
runs,  oil  mixture,  temperature  attained  at  end  of  decreasing 
heat  flux  phase  of  the  cycle,  and  variation  with  cool -down 
procedure.  Because  of  these  numerous  influences  on  the  boiling 
incipience  superheat,  past  studies,  including  Bar- Cohen  [Ref. 
3] ,  reveal  highly  variable  and  widely  distributed  superheats 
at  ONB .  Furthermore,  the  process  of  first  nucleation  is 
inherently  stochastic,  leading  to  variations  in  behavior  from 
one  day  to  the  next . 

3.   Nucleate  Boiling  Region 

Once  nucleate  boiling  commences,  during  increasing 
heat  flux  (point  c) ,  sustained  improvement  in  heat  transfer  is 
realized  throughout  the  entire  fully-developed  nucleate 
boiling  region  (c-d) .  As  vapor  bubbles  form  at  nucleation 
sites  and  separate  from  the  surface,  dramatic  increases  in 
heat  transfer  are  attained.  There  are  several  factors 
affecting  where  this  region  is  located  on  the  pool  boiling 
curve,  the  most  important  of  which  is  the  number  of  active 
sites  on  the  surface,  as  may  occur  on  different  boiling 
surfaces.  Thome  [Ref.  5]  provides  detailed  insight  into 
enhanced  boiling  heat  transfer  from  various  commercially 
available  surfaces.  For  plain  tubes,  several  correlations  have 
been  proposed  to  predict  boiling  performance.   Two  such 


correlations  were  provided  by  Stephan  and  Abdelsalam  [Ref .  6] 
and  Chongrungreong  and  Sauer  [Ref.  7]  . 

B.   SINGLE  TUBE  BOILING  EXPERIMENTS  AT 
NAVAL  POSTGRADUATE  SCHOOL 

Heat  transfer  measurements  were  collected  by  McManus 
[Ref .8]  for  single  tube  pool  boiling  of  R-114  at  a  saturation 
temperature  of  13.8°C  with  up  to  6%  oil  concentration.  The 
boiling  tubes  in  these  tests  were  heated  internally  by  warm 
water.  HIGH  FLUX,  TURBO-B  and  smooth  tubes  were  tested  in  this 
manner.  McManus  [Ref.  8]  reported  enhancements  (compared  to 
the  smooth  tube)  in  pure  R-114  as  high  as  14.6  and  6.4  for  the 
HIGH  FLUX  and  TURBO-B  tubes  respectively,  while  performance 
decreased  to  enhancements  of  7.0  and  4.9  respectively  with  a 
6%  oil  concentration. 

Pool  boiling  heat  transfer  coefficients  in  pure  R-114  and 
R-114/oil  mixtures  up  to  10%  (by  weight)  were  measured  by 
Wanniarachchi  et  al  [Ref.  9]  for  electrically-heated  smooth 
and  enhanced  tubes.  These  tests  were  run  at  a  saturation 
temperature  of  2.2°C.  At  30  kW/m2  Wanniarachchi  et  al  [Ref.  9] 
reported  enhancement  results  for  the  HIGH  FLUX,  Thermoexcel -E, 
Thermoexcel-HE,  and  GEWA-T  tubes  of  9.1,  8.2,  6.8,  and  4.4 
respectively.  Also  reported  were  the  enhancements  for  these 
tubes  with  3%,  and  10%  oil  concentration. 

Sugiyama  [Ref.  2]  conducted  extensive  pool  boiling  tests 
with  R-114  and  R-114/oil  mixtures  using  a  large  variety  of 
enhanced  surface  tubes.  For  the  re-entrant  cavity  tubes  with 
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pure  R-114,  Sugiyama  [Ref .  2]  reported  enhancements  of  6.4  and 
6.2  for  the  HIGH  FLUX  and  TURBO-B  tubes  respectively  at  35 
kW/m;.  With  10%  oil  these  enhancements  were  reduced  to  3 . 8  and 
5.6  respectively.  This  degradation  in  performance  was 
partially  attributed  to  the  clogging  of  surface  pores  with 
oil.  For  GEWA-K  26  and  40  fpi  tubes,  enhancements  of  2.4  and 
3.1  were  reported  with  pure  R-114  at  35  kW/m:.  With  10%  oil 
these  same  tubes  were  surprisingly  reported  to  have 
enhancements  of  3.2  and  4.6  respectively.  This  indicated  an 
improved  heat  transfer  performance  with  addition  of  oil  in  the 
case  of  finned  tubes,  believed  to  be  caused  by  foaming  action 
created  within  the  pool  that  was  most  prevalent  at  higher  oil 
concentrations  (3-10%)  and  heat  fluxes  (>30  kW/m2)  . 

The  effect  of  oil  on  the  incipience  of  nucleate  boiling 
for  pure  R-114  and  R-114/oil  mixtures  was  studied  by  Memory 
and  Marto  [Ref.  10]  with  smooth,  finned  and  re-entrant  cavity 
tubes.  They  observed  a  significantly  lower  heat  flux  at 
incipience  for  the  HIGH  FLUX  and  TURBO-B  tubes  compared  with 
the  smooth  and  finned  tubes  in  pure  R-114.  However,  consistent 
delay  in  the  boiling  incipience  of  the  HIGH  FLUX  and  TURBO-B 
tubes  was  exhibited  in  the  presence  of  3  and  10%  oil 
concentrations.  At  these  oil  concentrations,  boiling 
incipience  for  the  smooth  and  finned  tubes  showed  no 
dependence  with  oil.  Also  emphasized  was  the  inherent 
randomness  and  variation  connected  with  temperature  overshoot 
associated  with  boiling  incipience. 


C.   COMPARISON  OF  R-114  TO  R-124 

1.   Property  Comparison  of  R-114/R-124 

Table  1  is  provided  to  show  direct  comparison  of  the 
properties  of  R-114  and  R-124  at  Tsat=2.2°C 

Table  1.  R-114 /R-124  PROPERTY  COMPARISON 


Property 

R-114 

R-124 

Psat/  (kPa) 

96 

177 

p(liq) / (kg/m2) 

1526 

1429 

p (vapor) / (kg/m3) 

7.5 

11.3 

cp(liq)/ (J/kg-K) 

932 

1060 

hfg/(kJ/kg) 

134.7 

158.9 

k(liq) / (W/m-K) 

0.0701 

0.0744 

H (liq) / (g/m- s) 

0.4449 

0.3375 

a/ (N/m) 

0.0136 

0.0129 

2.   Need  for  Oil  Mixtures 

The  oil  mixtures  used  are  to  simulate  the  actual 
environment  in  refrigeration  systems  which  can  collect  lube 
oil  within  a  unit's  flooded  evaporator  to  concentrations  as 
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high  as  10%.  The  oils  used  for  R-124  and  R-114  are  both 
miscible  within  their  respective  refrigerants.  However,  with 
R-114,  a  York-C  mineral  oil  was  used,  while  with  R-124,  an 
alkylbenzene  oil  was  used. 
3.   Need  for  R-124  Data 

There  is  presently  no  data  on  single  tube  pool  boiling 
heat  transfer  with  R-124  or  R-124/oil  mixtures.  Tests  must 
therefore  be  conducted  with  R-124  to  parallel  R-114  results 
obtained  to  date. 


11 


III.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

A.   SYSTEM  DESCRIPTION 

Details  of  the  relocation  and  early  modifications  to  the 

original  apparatus  are  given  by  Sugiyama  [Ref.  2]  and  Lake 

[Ref .  11] .   To  accommodate  the  increased  operating  pressures 

associated  with  R-124,  stronger  evaporator  and  condenser  units 

were  installed.  The  data  reduction  program  DRP71  was  modified 

(DRPGB)  to  include  the  thermophysical  properties  of  R-124. 

A  schematic  diagram  representing  the  experimental 
apparatus  is  shown  in  Figure  3.1  and  comprises  the  following 
components:  1)  an  evaporator  tee  used  to  boil  the  refrigerant; 
(2)  a  condenser  tee  used  to  condense  the  refrigerant;  (3)  a 
reservoir  for  liquid  refrigerant;  (4)  a  refrigerant/oil 
subsystem; (5 )  a  refrigeration  cooling  support  system;  (6)  a 
vacuum  pump;  (7)  a  data  acquisition  and  instrumentation 
system;  (8)  an  aluminum/plexiglas  framework  housing  items  (1) 
to  (4)  . 

The  general  operation  of  the  system  was  closed  loop  with 
vapor  generated  in  the  evaporator  by  the  heated  tube  located 
within  the  refrigerant  pool  of  the  evaporator.  The  vapor 
passed  through  an  aluminum  'L-shaped'  tube  to  the  condenser 
where  the  vapor  was   condensed  by  a  chilled  water/ethylene 
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Figure  3.1  Schematic  of  Apparatus 
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glycol  mixture  circulated  through  the  condenser  coil  by  an  8 
gpm  turbine  type  positive  displacement  pump.  The  R-114  or  R- 
124  condensate  was  then  gravity  fed  back  to  the  evaporator 
through  a  return  line.  R-502  and  R-12  refrigeration  systems 
maintained  the  water/ethylene  glycol  mixture  in  the  sump 
between  -12  and  -18  degrees  Celsius.  Sump  cooling  was 
accomplished  using  a  counter  flow  heat  exchanger  by 
recirculating  the  sump  coolant  (using  an  identical  positive 
displacement  pump)  over  the  R-502  side  of  the  heat  exchanger. 
To  further  maintain  this  coolant  temperature,  the  R-12  system 
was  located  directly  within  the  water/ethylene  glycol  sump. 

The  addition  of  oil  to  the  refrigerant  within  the 
evaporator  was  accomplished  by  using  a  graduated  cylinder  to 
measure  the  precise  quantity.  This  simply  drained  into  the 
evaporator  by  gravity. 

Accurate  horizontal  positioning  of  the  boiling  tube  was 
critical  and  accomplished  by  supporting  each  end  of  the  tube 
by  Teflon  inserts.  The  tube  was  sealed  using  0- rings.  For 
safety,  a  relief  valve  on  the  aluminum  pipe  between  the 
evaporator  and  condenser  was  set  at  a  gage  pressure  of  49  psi. 

B.   BOILING  TEST  SECTION 
1.   Evaporator 

The  evaporator  was  a  Pyrex  glass  'T'  section  coated 
with  a  continuous  winding  of  fiberglass  filaments  impregnated 
with  a  modified  epoxy  resin.  This  upgraded  evaporator  shown  in 
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Figure  3.2  was  designed  for  a  working  pressure  of  50  psi  gage 
and  initially  hydrostatically  safety  tested  to  75  psi  gage. 
The  original  glass  was  designed  for  a  working  pressure  of  30 
psi  gage.  The  inside  surface  was  the  same  as  before  and 
maintained  the  advantages  of  minimizing  nucleate  boiling  at 
the  surface  as  well  as  being  corrosion  resistant.  The  two  ends 
of  the  evaporator  were  fitted  with  a  cast-iron  flange  and 
gasket  assembly  as  shown  in  Figure  3.2.  They  mated  with 
aluminum  flanges  at  each  end  of  the  evaporator,  which  in  turn 
contained  the  thermocouple  housings,  oil  entry  connection  and 
Teflon  inserts  containing  the  evaporator  tube. 
2 .   Evaporator  tubes 

Figure  3.3  shows  a  schematic  of  a  boiling  tube.  These  were 
accurately  held  in  place  using  Teflon  inserts  and  sealed  by  0- 
rings .  A  4 -stud  strong-back  secured  the  Teflon  inserts  inside 
the  aluminum  flanges. 

Five  different  tubes  were  tested  within  the  apparatus. 
Dimensions  of  each  tube  are  provided  in  Table  5  of  Appendix  C 
while  an  expanded  view  of  the  enhanced  surfaces  are  shown  in 
Figure  3.4.  The  tubes  used  were: 

1 .  Smooth  tube 

2.  GEWA-K  (19  fpi)  finned  tube 

3.  GEWA-K  (26  fpi)  finned  tube 

4.  HIGH  FLUX  (porous  coated)  enhanced  tube 

5.  TURBO-B  (deformed  surface)  enhanced  tube 
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a)  Corning  Pyrex  Glass  Evaporator  (D  x  Q.  =  402x51  mm, 
L  =  178  mm,  Li=  127  mm) 
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b)  Cast  Iron  Flange  and  Gasket  (d}  =  190mm,  d0  =  210mm, 
L  =  14mm,  A  =  21°  ) 


Figure  3.2  Sketch  of  Pyrex  Glass  Vessels  and  Flanges 


16 


^§^ 


Figure  3.3   Schematic  of  boiling  tube 
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Pigure   3.4   Surface   Profiles   of   Enhanced  Tubes 
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The  heater  used  within  each  boiling  tube  was  nominally  a 
1  kW,  240-Volt  stainless-steel  cartridge  heater,  measuring 
6.35  mm  in  outside  diameter  and  203.2  mm  in  length.  The  actual 
heated  length  of  the  heater  was  190  mm.  The  HIGH  FLUX  and 
GEWA-K  26  fpi  tubes  were  previously  manufactured  while  the 
TURBO- B,  GEWA-K  19  fpi  and  smooth  tubes  required  assembly 
prior  to  testing  with  R-124  or  R-114.  The  manufacturing 
process  of  these  tubes  is  described  below. 

The  heater  was  placed  inside  a  copper  sleeve  (Figure  3.5 
provides  dimensions  of  smooth  tube  sleeve)  with  a  tight 
mechanical  fit  of  .001  inch  clearance.  The  copper  sleeve  was 
machined  with  6  grooves  (1.27mm  x  1.27mm)  on  its  outer 
diameter  of  various  lengths.  Thermocouples  were  laid  in  each 
of  these  grooves  and  the  edges  of  the  grooves  were  peened  over 
with  a  punch.  Any  resulting  imperfections  were  smoothed  oui 
with  250  grit  emery  paper.  The  thermocouple  wires  and  heater 
leads  all  exited  the  copper  sleeve  at  the  same  end.  The  copper 
sleeve  containing  the  heater  and  thermocouples  was  then  placed 
on  2  aluminum  V-blocks  and  the  outer  surface  of  the  sleeve  was 
covered  with  a  flux  and  solder  paste  (50/50  tin/lead);  great 
care  was  taken  to  ensure  complete  coverage  of  the  thermocouple 
wires  lying  in  the  grooves.  The  thermocouple  and  heater  leads 
were  then  passed  through  the  enhanced  boiling  tube  (the  outer 
diameter  of  the  copper  sleeve  was  machined  to  fit  tightly  into 
each  respective  tube) .  To  prevent  shorting  between  the  heater 
and  outer  tube  shell,  a  high  temperature  RTV  silicone  rubber 
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was  applied  to  the  heater  and  thermocouple  leads  at  the  point 
where  they  leave  the  sleeve.  The  copper  sleeve  was  then 
tapped  gently  into  the  tube  and  positioned  accurately  so  as  to 
line  up  precisely  with  the  8  inch  active  length  of  the 
enhanced  surface.  This  enhanced  surface  length  was  different 
from  the  actual  heated  length  which  was  only  7.5  inches. 

The  entire  assembly  was  then  placed  horizontally  in  an 
oven  and  heated  to  between  400°F  and  500°F.  This  was  the 
temperature  'window'  at  which  the  flux  melted,  but  the 
thermocouple  leads  did  not.  The  temperature  of  the  tube  was 
monitored  by  attaching  the  thermocouples  to  a  simple 
temperature  indicator.  After  stabilizing  the  oven  temperature 
between  these  temperatures  and  ensuring  that  the  solder  had 
melted,  the  tube  was  taken  out  and  quenched  with  a  damp  cloth 
at  one  end  (as  close  to  the  end  as  possible)  where  the 
thermocouple  and  heater  leads  protrude.  The  oven  was  then 
placed  vertically  and  the  tube  reinserted  from  below  keeping 
the  lower  end  quenched  while  holding  both  tube  and  cloth  with 
a  pair  of  channel  locks.  When  the  temperature  increased  to 
>400°F,  the  solder  again  began  to  melt.  Additional  40/60  tin- 
lead  solder  was  then  applied  by  hand  from  the  end  which 
protrudes  from  the  top  of  the  oven.  Remelting  and  addition  of 
extra  solder  in  this  way  ensured  a  uniform  circumferential 
layer  of  solder  and  helped  to  prevent  formation  of  air 
pockets.  The  tube  was  then  removed  from  the  oven  and  quenched 
over  its  whole  length  with  a  damp  cloth  while  in  a  vertical 
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position.  It  was  then  left  to  cool  in  air  for  24  hours.  This 
procedure  differed  from  earlier  manufactured  tubes  in  that  the 
previous  copper  sleeves  had  4  thermocouples  leaving  each  end 
of  the  sleeve.  Trying  to  repeat  this  previous  procedure 
accurately  led  to  large  discrepancies  in  the  thermocouple 
readings  (>20°C)  as  heater  power  was  increased.  It  was  felt 
this  was  due  to  air  pockets  created  in  the  fabrication 
procedure  (the  tube  could  never  be  placed  vertically  in  the 
oven  with  leads  protruding  from  both  ends) .  The  above 
procedure  was  much  simpler  and  led  to  excellent  uniformity  in 
the  thermocouple  values  at  all  heater  powers. 

C.   CONDENSER  SECTION 

The  condenser  was  also  a  fiberglass  coated  Corning  Pyrex 
glass  Tee  identical  in  size  and  shape  to  the  evaporator  shown 
in  Figure  3.2.  The  refrigerant  was  condensed  on  a  3/8  inch 
outside  diameter  copper  tube  that  was  coiled  within  the 
condenser.  The  condenser  was  mounted  vertically.  An  O-ring 
sealed  aluminum  flange  mounted  at  the  top  of  the  condenser 
contained  not  only  the  water  ethylene/glycol  coolant 
inlet/outlet,  but  also  a  3/8  inch  outside  diameter  copper  tube 
vent,  which  was  connected  either  to  a  vacuum  pump  via  valve 
Vll  or  to  a  charging  cylinder  via  valve  V10.  The  bottom  of  the 
condenser  was  also  capped  with  an  O-ring  sealed  aluminum 
flange  and  contained  the  gravity  feed  return  line  to  the 
evaporator  via  valve  V5 .  The  refrigerant  vapor  from  the 
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evaporator  entered  the  condenser  through  an  aluminum  L- shaped 
tube,  connected  to  the  condenser  by  O-ring  and  gasketed 
flange.  A  pressure  gage  ranging  from  3  0  inches  mercury  vacuum 
to  150  psig  was  mounted  on  the  L- shaped  tube;  this  was  used 
only  as  a  rough  visual  check  of  pressure  within  the  facility. 

D.  OIL  ADDITION  SECTION 

To  provide  oil  addition  into  the  refrigerant,  a 
cylindrical  aluminum  reservoir,  6  inches  in  diameter  and  6 
inches  high,  and  a  graduated  cylinder  355  mm  in  length  and 
25.4  mm  in  diameter  is  used  as  shown  in  Figure  3.1.  The  oil 
used  with  R-124  was  a  miscible  alkylbenzene  oil  while  the  oil 
used  with  R-114  was  a  York-C  mineral  oil,  also  miscible.  The 
graduated  cylinder  was  filled  from  the  reservoir  via  valve  V2  , 
while  the  reservoir  was  replenished  through  valve  V3  at  the 
top.  Oil  entered  the  evaporator  from  the  graduated  cylinder  by 
gravity  feed  via  valve  VI.  Oil  concentrations  (by  weight)  of 
0,1,2,3,6,  and  10%  were  used  for  each  of  the  tubes. 

E.  COOLING  SYSTEM 

1.   Water/Ethylene- Glycol  Coolant  Tank 

At  the  base  of  the  apparatus,  a  5.35  cubic  foot  tank 
stored  30  gallons  of  ethylene-glycol/water  mixture. 
Recirculation  and  discharge  ports  were  connected  to  the  tank 
together  with  a  thermocouple  to  monitor  sump  temperature. 
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2.  Refrigeration  Plants 

To  provide  the  cooling  of  the  water/ethylene  glycol 
mixture,  a  1/2  ton  R-502  and  1/4  ton  R-12  refrigeration  system 
were  installed  near  the  sump  as  seen  in  Figure  3.1.  Both 
plants  contained  a  hermetically  sealed  compressor  assembly, 
an-air  cooled  condenser,  receiver,  filter-dryer,  pressure 
regulator,  temperature  control  switch,  and  thermostatic 
expansion  valve.  The  R-502  plant  (Figure  3.6)  used  a  counter 
current  heat  exchanger  while  the  R-12  plant  used  a  coiled 
evaporator  placed  within  the  sump.  The  sump  temperature  was 
maintained  at  approximately  -17  °C  when  both  these  systems 
were  used.  Each  plant  was  controlled  by  a  temperature  control 
switch  and  thermostatic  expansion  valve. 

3 .  Pump  and  Control  Valve 

An  8  gpm,  positive  displacement  pump  was  installed 
with  a  1  inch  diameter  PVC  pipe  connected  to  the 
water/ethylene  glycol  sump.  The  discharge  of  the  pump  was 
piped  to  the  condenser  via  valve  VC  as  shown  in  Figure  3.1.  A 
by- pass  valve  V9  was  placed  upstream  of  the  control  valve  VC 
to  avoid  overloading  the  pump  when  control  valve  VC  was  closed 
and  provides  for  mixing  of  the  return  coolant  from  the 
condenser.  At  high  heat  fluxes,  the  by-pass  valve  needed  to  be 
closed  down  to  increase  flow  of  coolant  through  the  condenser. 
An  identical  pump  discharged  the  ethylene  glycol/water  coolant 
through  the  R-502  counter- current  heat  exchanger. 
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Figure  3.6   R-502  Refrigerant  plant 
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F.  REFRIGERANT  RESERVOIR 

The  reservoir  used  to  store  the  R-114  or  R-124  was  made 
from  an  aluminum  cylinder,  9  inches  in  diameter  and  10  inches 
in  height.  A  sight  glass  attached  on  the  reservoir  allowed  for 
the  liquid  level  of  refrigerant  to  be  monitored.  The  reservoir 
was  connected  to  the  ' L'  shaped  tube  through  valve  V7  and  to 
the  liquid  line  into  the  evaporator  through  valve  V6  as  shown 
in  Figure  3.1. 

G.  FRAME 

To  house  the  entire  assembly,  an  aluminum  frame  was 
constructed  measuring  1070.1  mm  x  509.8  mm  x  610.1  mm.  With  a 
view  to  keeping  the  apparatus  in  an  isothermal  atmosphere,  as 
well  as  for  safety  considerations,  the  sides  were  covered  with 
12 . 7  mm  thick  plexiglas  with  hinges  mounted  on  the  2  sides  to 
provide  easy  access  to  the  assembly.  The  top  of  the  frame  was 
covered  by  plywood  while  the  bottom  was  covered  by  an  aluminum 
plate.  The  entire  frame  rested  on  the  coolant  sump  tank.  The 
front  side  of  plexiglas  had  small  holes  to  accommodate  the 
valve  stems  of  valves  VI  through  V8 ,  so  they  can  be  operated 
while  maintaining  the  integrity  of  the  enclosed  frame. 

H.   INSTRUMENTATION 

1.   Power  Measurement 

The  boiling  tube  heater  was  powered  by  a  240  volt  AC 
source,  adjusted  to  0-220  volts  and  0-5  amps  by  a  variac 
control.   Power  input  to  the  heater  was  measured  with  an  AC 
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Figure  3.7   Data  acquisition/control  unit 
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current  inductive  sensor  and  a  voltage  sensor.  The  output  was 
a  proportional  signal  in  volts  having  been  converted  through 
an  AC -DC  true  R.M.S.  converter.  Both  the  AC  current  sensor  and 
the  AC-DC  true  R.M.S.  converter  were  connected  to  the  data 
acquisition/control  unit  as  shown  in  Figure  3.7. 
2 .   Temperature  Measurement 

The  primary  means  of  obtaining  temperatures  was  the 
use  of  30  gage  copper- constantan  thermocouples.  Specific 
temperatures  recorded  were  from  6  thermocouples  in  the  boiling 
tube  (8  in  the  earlier  manufactured  tubes)  ,  3  thermocouples  in 
the  evaporator  pool  (housed  within  separate  sealed 
thermocouple  wells  within  the  evaporator  pool),  and  one 
thermocouple  to  monitor  the  water/ethylene  glycol  sump  tank. 
A  Hewlett-Packard  3497A  data  acquisition  system  controlled  by 
a  Hewlett-Packard  9826  computer  read  these  thermocouple 
measurements . 

I.   DATA  ACQUISITION  AND  REDUCTION  PROGRAM 

Data  collection  and  reduction  program  DRP71,  used  by 
Sugiyama  [Ref .2]  for  use  with  R-114  was  modified  for  use  with 
R-124.  The  thermophysical  properties  of  saturated  R-124  were 
obtained  from  Ref  Prop  software  [Ref.  12]  and  were  used  to 
generate  equations  for  use  in  the  calculations  of  the 
reduction  program.  Appendix  A  contains  plots  of  each  property 
versus  temperature  and  the  respective  equation.  The  modified 
program  was  renamed  DRPGB. 
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The  Hewlett-Packard  9826  computer  reduced  all  the 
thermocouple,  current,  and  R.M.S.  voltage  readings  using  the 
DRPGB  program.  The  program  was  run  through  the  keyboard 
allowing  regulation  of  the  data  acquisition  unit.  A  printout 
of  the  reduced  data  was  provided  by  a  Hewlett-Packard  Inkjet 
printer.  Table  2  lists  the  channel  assignments  for  the  sensor 
inputs  to  the  data  acquisition  system. 

J.   PROCEDURE  FROM  KEYBOARD 

The  following  procedure  outlines  the  procedure  used  to 
collect  and  process  data  from  the  Hewlett  Packard  data 
acquisition  system: 


1.  Selected  the  'taking  data'  option  after  initial 
loading  of  program  DRPGB. 

2.  Selected  the  refrigerant  being  tested,  tube  heating 
mode,  (ie  electrically  or  water  heated)  and 
thermocouple   type  from  DRPGB. 

3.  Assigned  name  to  file  of  data  to  be  processed  and 
stored. 

4.  Identified  and  entered  any  defective  thermocouples  by 
channel  #  1-8,  noting  that  the  new  tubes  have  only  6 
thermocouples  and  channels  7  &  8  were  always 
' defective' . 

5.  Selected  test  tube  type. 

6.  Set  desired  saturation  temperature  (2.22°C)  for 
present  tests. 

7.  Set  desired  heat  flux  desired  using  variac  adjustment. 

8.  Attained  desired  saturation  temperature  by  adjusting 
flow  of  coolant  through  condenser  coils  with  control 
valve  VC. 
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9.  Once  saturation  temperature  was  achieved,  waited  5 
minutes  for  steady  state  conditions  prior  to  taking 
data . 

10.  Prompted  data  acquisition  unit  to  scan  all  the 
channels  listed  in  Table  2.  All  channel  readings  were 
made  in  volts  and  stored  in  user  specified  fields. 

11.  Calculated  temperature  and  power  from  these  voltages. 

12.  Calculated  heat  transfer  rate  from  cartridge  heater. 

13.  Calculated  the  average  wall  temperature  of  the  test 
tube  and  calculated  wall  superheat  (Twall-  Tsat) . 

14.  Calculated  physical  properties  of  R-124  using  property 
correlations  f rom  RefProp  [Ref.12]  at  film  temperature 
(Twall+Tsat) /2. 

15.  Calculated  the  natural  convection  heat- transfer 
coefficient  of  R-124  from  the  unheated  ends  of  the 
test  tube. 

16.  Computed  heat  losses  from  unheated  ends  of  the  test 
tube . 

17.  Calculated  heat  flux  from  heated  length  of  the  test 
tube . 

18.  Calculated  the  heat- transfer  coefficient  from  the 
heated  length  of  the  test  tube. 

19 .  Stored  heat  flux  and  wall  superheat  for  each  data  set 
in  user  specified  fields. 

20.  Plotted  data  on  available  software. 


Sample  calculations  of  the  above  procedure  are  given  in 
Appendix  C. 
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TABLE  2.   HP  3497A  CHANNEL  ASSIGNMENTS 


Channel 

Channel  Assignments 

00-07 

tube  wall  temperatures 

08-10 

pool  liquid  temperature 

11 

sump  temperature 

20 

RMS  voltage 

21 

current  sensor 
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IV.  EXPERIMENTAL  PROCEDURE 

A.   ASSEMBLY  PREPARATION 

1.  Vacuum  test  of  the  apparatus 

After  placing  the  desired  boiler  cube  in  position 
within  the  evaporator,  the  system  was  evacuated  down  to  a 
pressure  of  29  inch  Hg .  vacuum  using  a  portable  vacuum  pump 
via  valves  Vll  and  V8  in  Figure  3.1.  The  system  was  then 
isolated  for  approximately  30  minutes  to  check  for  leak- tight 
integrity  (indicated  by  any  significant  vacuum  loss)  .  If  a 
leak  was  detected,  a  positive  air  pressure  (20  psig)  was 
placed  on  the  system  via  valves  V10  and  V8 .  A  soap/water 
solution  was  applied  on  the  sealing  areas  to  identify  any 
leaks  (if  there  is  a  leak,  bubbles  are  seen) .  If  identified, 
the  leak  source  was  corrected  and  the  system  integrity 
retested. 

2.  Charging  apparatus  with  refrigerant 

After  system  integrity  was  ensured,  the  evaporator  was 
filled  with  refrigerant  to  a  level  of  20  mm  above  the  test 
tube.  The  refrigerant  was  drained  by  gravity  from  the 
reservoir  into  the  evaporator  through  valves  V6  and  V5  .  Valve 
V7  was  also  open  to  equalize  pressure  in  the  system.  The 
reservoir  was  then  isolated  from  the  system  with  the  closure 
of  valves  V5  through  V7 . 


32 


3.   Degassing  and  Data  Acquisition  Channel  Check 

Prior  to  connecting  all  thermocouple  leads,  the 
resistance  between  the  heater  leads  and  outer  tube  wall  was 
checked  with  a  simple  multimeter.  Readings  were  usually 
greater  than  1  Mt2  .  Each  channel's  output  was  checked  using 
program  SETUP71.  Any  problem  thermocouples  were  either  fixed 
or  excluded  from  the  calculations  in  program  DRPGB .  Before 
energizing  the  heater  and  starting  the  tests,  the  evaporator 
pool  was  slowly  cooled  to  a  saturation  pressure  of  2.22  °C  by 
allowing  a  low  coolant  flow  through  the  condenser.  Once 
completed,  the  heating  surface  flux  was  set  to  approximately 
90  kW/m2  to  degas  the  refrigerant;  any  non- condensible  gases 
were  collected  in  the  condenser  and  vented  off  using  the 
vacuum  pump  for  approximately  one  minute.  The  apparatus  was 
then  secured  for  a  period  of  approximately  8  hours  prior  to 
operation  to  allow  the  surface  to  cool  and  become  fully  wetted 
in  the  refrigerant. 

B.   OPERATIONAL  PROCEDURE 

The  following  procedure  was  used  to  operate  the  facility 
(this  follows  closely  the  procedure  used  by  Sugiyama  [Ref .2] ) : 

1.  The  R-502  and  R-12  refrigeration  units  were  utilized 
until  the  temperature  of  the  water- ethylene  glycol  sump 
reached  a  minimum  of  -  8  °C  before  starting  tests  (with 
continued  cooling  to  -17  °C) . 

2.  The  data  acquisition  unit,  computer  and  variac  control 
panel  were  switched  on. 
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3.  The  computer  program  SETUP  was  loaded  and  run.  All  data 
acquisition  channels  were  rechecked.  The  refrigerant 
pool  temperature  was  then  cooled  down  slowly  to  2.2  °C 
by  circulating  coolant  through  the  condenser. 

4.  The  data  acquisition/reduction  program  DRPGB  was  loaded 
and  run. 

5.  The  desired  heat  flux  setting  for  the  test  tube  was 
input  to  the  program  DRPGB  and  the  variac  was  adjusted 
to  attain  that  heat  flux.  Starting  heat  flux  was 
between  500  and  600  W/m2 , 

6.  Control  valve  VC  was  adjusted  to  control  coolant  flow 
through  the  condenser  to  maintain  a  constant  saturation 
temperature  of  2.22  °C  +  0.1  °C . 

7.  At  each  heat  flux  setting,  the  thermocouple  values  were 
allowed  to  stabilize  for  about  5  minutes  prior  to 
taking  data. 

8 .  The  following  data  were  then  taken  and  used  in  the 
acquisition/reduction  program  DRPGB:  pool  temperature; 
test  tube  wall  temperatures;  sump  temperature;  current 
and  voltage  sensor  readings. 

9.  Data  were  taken  at  each  desired  heat  flux  setting  and 
always  at  the  same  saturation  temperature  (the  coolant 
flow  in  the  condenser  was  adjusted  to  ensure  this) .  At 
each  heat  flux  setting,  program  DRPGB  generated  a 
permanent  printout  listing:  wall  temperatures  of  the 
test  tubes;  refrigerant  pool  temperatures;  sump 
temperature;  wall  superheat;  heat  -  transfer  coefficient 
and  the  corrected  heat  flux  (ie.  to  account  for  heat- 
transfer  lost  from  each  end) . 

10.  For  each  heat  flux  setting,  steps  6  through  9  were 
repeated.  The  superheat  (Twall-Tsat)  and  the  respective 
heat  flux  were  plotted  on  a  log-log  scale. 

11.  Table  3  provides  a  summarized  listing  of  all  data  runs. 
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TABLE  3 .  LISTING  OF  DATA  RUNS 


Data  File 

Refrigerant 

Tube  Type 

Purpose 

D0821SM0 

R-114 

Smooth 

repeatability 

D0824SM0 

R-114 

Smooth 

repeatability 

D1015SM0 

R-124 

Smooth 

Data 

D1017SM0 

R-124 

Smooth 

Data 

D1223SM0 

R-124 

Smooth 

Data 

D1222SM0 

R-124 

Smooth 

Data 

D1221SM0 

R-124 

Smooth 

Data/samp  calc 

D1220SM0 

R-124 

Smooth 

uncert  analys 

D1219SM0 

R-124 

Smooth 

Data 

D1217SM0 

R-124 

Smooth 

Data 

D1016SM4 .44* 

R-124 

Smooth 

Data  Tsat4.4°C 

D1114SM0 

R-124 

Smooth 

Data 

D1114SM1 

R-124 

Smooth 

Data 

D1114SM2 

R-124 

Smooth 

Data 

D0108SM3 

R-124 

Smooth 

Data 

D0107SM3 

R-124 

Smooth 

Data 

D0106ASM3 

R-124 

Smooth 

Data 

D0106SM3 

R-124 

Smooth 

Data 

D0105SM3 

R-124 

Smooth 

Data 

D0104SM3 

R-124 

Smooth 

Data 

D1227SM3 

R-124 

Smooth 

Data 

D1116SM3 

R-124 

Smooth 

Data 

D1116SM6 

R-124 

Smooth 

Data 

D0112ASM10 

R-124 

Smooth 

Data 

D0112SM10 

R-124 

Smooth 

Data 

D0111ASM10 

R-124 

Smooth 

Data 

D0111SM10 

R-124 

Smooth 

Data 
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Continuation  of 

Table  3. 

Data  File 

Refrigerant 

tube  type 

Purpose 

D0110SM10 

R-124 

Smooth 

Data 

D0109ASM10 

R-124 

Smooth 

Data 

D0109SM10 

R-124 

Smooth 

Data 

D1117SM10 

R-124 

Smooth 

Data 

D1011190 

R-114 

GK-19 

Data 

D1010190 

R-114 

GK-19 

Data 

D0901190 

R-114 

GK-19 

Data 

D0907193 

R-114 

GK-19 

Data 

D09091910 

R-114 

GK-19 

Data 

D1105190 

R-124 

GK-19 

Data 

D1118190 

R-124 

GK-19 

Data 

D1118191 

R-124 

GK-19 

Data 

D1118192 

R-124 

GK-19 

Data 

D1119193 

R-124 

GK-19 

Data 

D1120196 

R-124 

GK-19 

Data 

D11211910 

R-124 

GK-19 

Data 

D0803260 

R-114 

GK-2  6 

Repeatability 

D1021260 

R-124 

GK-26 

Data 

D1122260 

R-124 

GK-26 

Data 

D1122261 

R-124 

GK-26 

Data 

D1122262 

R-124 

GK-2  6 

Data 

D1123263 

R-124 

GK-26 

Data 

D1123266 

R-124 

GK-26 

Data 

D11252610 

R-124 

GK-26 

Data 

D11242610 

R-124 

GK-26 

Data 

D1004TB0 

R-114 

TURBO  - B 

Data 

D1006TB0 

R-114 

TURBO  - B 

Repeatability 
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Continuation  of 

Table  3. 

Data  File 

Refrigerant 

Tube  Type 

Purpose 

D1002GYX0 

R-114 

GEWA-YX 

Data 

D1110TB0 

R-124 

TURBO -B 

Data 

D1110TB0 

R-124 

TURBO -B 

Data 

D1111TB1 

R-124 

TURBO -B 

Data 

D1111TB2 

R-124 

TURBO  - B 

Data 

D1112TB3 

R-124 

TURBO  - B 

Data 

D1112TB6 

R-124 

TURBO  - B 

Data 

D1113TB10 

R-124 

TURBO  - B 

Data 

D0119AHF0 

R-124 

HIGH  FLUX 

Data 

D0118BHF0 

R-124 

HIGH  FLUX 

Uncert  Analys 

D0118HF0 

R-124 

HIGH  FLUX 

Data 

D0114HF0 

R-124 

HIGH  FLUX 

Data 

D0113AHF0 

R-124 

HIGH  FLUX 

Data 

D0113HF0 

R-124 

HIGH  FLUX 

Data 

D1108HF0 

R-124 

HIGH  FLUX 

Data 

D1126HF0 

R-124 

HIGH  FLUX 

Data 

D1126HF1 

R-124 

HIGH  FLUX 

Data 

D1126HF2 

R-124 

HIGH  FLUX 

Data 

D0125HF3 

R-124 

HIGH  FLUX 

Data 

D0122HF3 

R-124 

HIGH  FLUX 

Data 

D0121CHF3 

R-124 

HIGH  FLUX 

Data 

D0121HF3 

R-124 

HIGH  FLUX 

Data 

D0120AHF3 

R-124 

HIGH  FLUX 

Data 

D0120HF3 

R-124 

HIGH  FLUX 

Data 

D0119AHF3 

R-124 

HIGH  FLUX 

Data 

D1128HF3 

R-124 

HIGH  FLUX 

Data 

D1128HF6 

R-124 

HIGH  FLUX 

Data 
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Continuation  of 

Table  3. 

Data  File 

Refrigerant 

Tube 

Type 

Purpose 

D1130HF10 

R-124 

HIGH 

FLUX 

Data 

D0128HF1C 

R-124 

HIGH 

FLUX 

Data 

D0127CHF10 

R-124 

HIGH 

FLUX 

Data 

D0127BHF10 

R-124 

HIGH 

FLUX 

Data 

D0127HF10 

R-124 

HIGH 

FLUX 

Data 

D0126BHF10 

R-124 

HIGH 

FLUX 

Data 

D0126AHF10 

R-124 

HIGH 

FLUX 

Data 

D0126HF10 

R-124 

HIGH 

FLUX 

Data 

The  tabulated  data  was  filed  using  the  following  file  name 
system: 

Example  (D0126AHF10) : 

1.  First  letter  D  simply  refers  to  data 

2.  The  following  4  characters  are  always  the  date;  in  this 
case  Jan  26. 

3 .  The  next  letter  represents  order  of  the  day  test  was 
completed  on  the  respective  date; 

no  letter-  1st  run 
A-  2nd  run 
B-  3rd  run 
etc  . 

4.  The  following  2  letters  represented  the  tube  used; 

SM=  smooth 

19=  19  fpi  GEWA-K 

26=  2  6  fpi  GEWA-K 

HF=  HIGH  FLUX 

TB=  TURBO -B 

5 .  The  last  1  or  2  numbers  represent  the  percentage  of 
oil;  in  this  case  10%. 

*   this  one  run  was  pure  R-124  tested  at  4.4  °C  vice  2.2  °C 
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V.   RESULTS  AND  DISCUSSION 

A.      REPRODUCTION  OF  EXPERIMENTAL  DATA 

The  first  tests  conducted  were  to  verify  the 
reproducibility  of  the  single  tube  data  taken  by  Sugiyama 
[Ref.  2]  with  R-114.  A  GEWA-K  26  fin  per  inch  (fpi)  tube 
comparison  with  R-114  is  given  in  Figure  5.1  and  shows  very 
similar  results  with  previous  data.  Uncertainty  bands  for  both 
wall  superheat  and  heat  flux  are  given  on  the  figure. 
Consequently,  at  low  heat  flux  (<  1000  W/m2)  ,  there  is  more 
scatter  due  to  the  lower  wall  superheats  and  hence  higher 
uncertainty.  Overall,  the  agreement  between  the  two  data  sets 
is  satisfactory. 

In  Figure  5.2,  an  increasing  heat  flux  data  run  using 
a  smooth  tube  is  compared  with  Sugiyama' s  [Ref.  2]  data.  It 
can  be  seen  that  in  the  low  heat  flux,  natural  convection 
region  and  the  high  heat  flux,  nucleate  boiling  region  the 
data  agreement  is  good.  In  the  mid- range  of  heat  flux,  (3 
kW/m2<  q  <  15  kW/m2 )  ,  the  data  varies  somewhat,  depending  upon 
how  the  tube  nucleates.  Curve  #1  nucleates  at  a  relatively  low 
wall  superheat  (12.5  °C)  while  curve  #2  (both  taken  by 
Sugiyama  [Ref.l])  nucleates  at  a  high  wall  superheat  of  about 
35  °C  (point  C) .   The  present  data  (curve  #3)  exhibits  partial 
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26  fpi  tube  (decreasing  flux) 
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nucleation  (see  point  A)  where  3  of  the  six  thermocouples  had 
a  low  wall  superheat  (<  14  °C)  and  the  remaining  3  had  a  wall 
superheat  (>  18  °C) .  Observing  the  tube  surface,  it  could  be 
clearly  seen  that  half  the  tube  was  boiling  while  the  other 
half  was  in  the  natural  convection  mode.  At  point  B,  all  six 
thermocouples  were  within  0.9  °C  of  each  other,  indicating 
complete  nucleation  over  the  entire  surface.  At  this  point 
(and  above),  agreement  with  the  previous  data  is  excellent. 
The  discrepancies  in  this  mid  heat  flux  region  demonstrate  the 
random  behavior  of  nucleate  boiling  incipience;  this  is 
further  investigated  later  in  this  chapter.  After  all  sites 
nucleated  (Point  B)  ,  the  standard  deviation  was  0  .  9  °C  with  an 
average  thermocouple  temperature  of  12.9°C.  The  point  of 
incipience  and  its  random  occurrence  at  different  levels  of 
superheat  will  be  examined  later  in  this  chapter. 

Sugiyama  [Ref  .2]  did  not  conduct  tests  on  a  19  fpi  GEWA- 
K  tube.  This  is  one  of  the  most  common  tubes  used  in 
commercial  evaporators.  Therefore,  base-line  data  for  the  19 
fpi  GEWA-K  tube  with  R-114  and  R-114/oil  mixtures  were  taken 
during  this  thesis  before  evacuating  and  cleaning  the 
apparatus,  and  filling  with  R-124.  The  R-114  data  obtained  was 
incorporated  in  the  R-124/R-114  comparison  of  the  19  fpi  tube 
(Figure  5  .  12)  . 

Figure  5.3  shows  the  repeatability  of  decreasing  heat 
flux  data  for  the  newly  constructed  19  fpi  GEWA-K  tube  using 
R-124.  The  good  agreement  is  an  indication  of  the  consistency 
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Figure  5.3  Repeatability  Comparison  of  R-124  Data  with  0%  Oil 
for  19  fpi  GEWA-K  Tube  (decreasing  flux) 
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of  the  data  reduction  program  DRPGB,  the  apparatus  and  the 
operator.  Decreasing  heat  flux  runs  provide  better  comparison 
as  temperature  overshoot  and  hysteresis  effects  are  avoided. 

B.   POOL  BOILING  HEAT-TRANSFER  COEFFICIENT  CORRELATIONS 

Among  the  correlations  available  to  predict  pool  boiling 
heat  transfer  coefficients  for  a  smooth  tube,  most  can  be 
expressed  as  heat  flux  as  a  function  of  wall  superheat.  The 
pure  R-124  smooth  tube  data  is  fitted  very  closely  by  the 
equation: 

q=17.0  aT3-0  (1) 

A  correlation  in  terms  of  fluid  physical  properties  and 
tube  geometry  is  given  by  Chongrungreong  and  Sauer  [Ref .  7] . 
For  R-124  at  177  kPa,  this  reduces  to: 

q=45.8  aT2-32  (2) 

Another  simpler  correlation  based  only  on  system  pressure 
is  also  given  by  Chongrungreong  and  Sauer  [Ref.  7] : 

q=76.9  aT2-22  (3) 

Stephan  and  Abdelsalem  [Ref.  6]  present  a  statistical 
correlation  for  four  different  groups  of  fluids.  For 
refrigerants,  their  correlation  is  presented  in  the  form: 

q-  (CIaT)3-92  (4) 

where  CI  is  given  as  a  function  of  refrigerant  type  and  system 
pressure.  Unfortunately,  as  R-124  is  a  relatively  new 
refrigerant,  their  graphical  representation  of  CI  does  not 
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include  R-124.  However,  the  present  data  could  be  fitted  to  an 
equation  of  the  form  given  in  equation  (4) : 

q=(1.34  AT) 3-92  (5) 

Hence,  at  a  pressure  of  1  atm,  Cl  =  1.34. 

Figure  5.4  shows  the  present  R-124  smooth  tube  data  for 
decreasing  heat  flux  along  with  the  above  correlations.  At 
lower  heat  fluxes  the  prediction  of  the  Chongrungreong  and 
Sauer  [Ref.  7]  correlation  (equation  (2)  above)  exhibits 
considerably  closer  prediction  than  the  other  correlations.  In 
the  more  complex  nucleate  boiling  region,  the  Stephan  and 
Abdelsalem  [Ref.  6]  correlation  and  the  Chongrungreong  and 
Sauer  [Ref.  7]  correlation  (equation  (3))  display  a  fair 
prediction  of  the  slope. 

C.   SMOOTH  TUBE  DATA  IN  REFRIGERANT/OIL  MIXTURES 

Although  tests  were  conducted  with  oil  concentrations  of 
0,1,2,3,6,  &  10%,  only  data  with  0,3  and  10%  have  been  plotted 
in  most  of  the  figures  for  clarity.  Figures  5.5  and  5.6  show 
the  heat  transfer  performance  of  pure  R-124,  R-124/3%  oil,  and 
R- 124/10%  oil  for  a  smooth  tube  during  increasing  and 
decreasing  heat  flux  respectively.  It  can  be  seen  that  at  the 
highest  heat  fluxes,  as  oil  concentration  is  increased  to  3%, 
the  heat  transfer  is  enhanced  by  a  maximum  of  17%  over  pure  R- 
124.  For  further  increases  in  oil  concentrations  to  10%, 
performance  drops  off.  This  drop-off  in  thermal  performance 
has  also  been   observed   by  Jensen  and  Jackman  [Ref.  13]  who 
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Figure  5.5  Performance  Comparison  for  Boiling  R-124/  0%,  3%, 
&  10%   Mixtures  for  Smooth  Tube  (increasing  flux) 
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Figure  5.6  Performance  Comparison  for  Boiling  R-124/  0%,  3%, 
&  10%  Mixtures  for  Smooth  Tube  (decreasing  flux) 
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attributed   the   decline   to   a   steeper   oil   concentration 
gradient,  inhibiting  bubble  growth. 

When  oil  is  mixed  with  R-124,  two  competing  mechanisms 
occur:  (1)  Oil  accumulates  in  the  vicinity  of  the  vapor  bubble 
surface  as  the  more  volatile  refrigerant  is  evaporated; 
(2)  diffusion  of  this  accumulated  oil  back  into  the  bulk 
mixture  establishes  an  equilibrium  concentration  gradient 
near  the  departing  bubble.  Mechanisms  (1)  and  (2)  should  act 
to  inhibit  bubble  growth  and  lead  to  an  overall  decrease 
in  the  heat  transfer  with  oil  addition  compared  with  pure 
refrigerant.  It  therefore  appears  that  there  is  an  additional 
mechanism  that  has  a  significant  impact  on  the  thermal 
performance  of  pool  boiling  of  refrigerants  with  small 
quantities  of  oil.  For  R-124/oil  mixtures  of  3%  and  greater, 
significant  foaming  in  the  evaporator  was  observed  at  mid  to 
high  heat  fluxes  (>  10000  kW/m2) .  It  is  possible  that  a  third 
mechanism  which  could  be  affecting  heat  transfer  is  this 
foaming  that  occurred  within  the  bulk  R-124/oil  mixture. 
Udomboresuwan  and  Mesler  [Ref .14]  reported  enhancement  in  pool 
boiling  in  the  presence  of  foam.  They  attributed  this 
enhancement  to  the  foam  causing  the  liquid/vapor  interface  to 
be  closer  to  the  tube  wall.  This  would  create  a  thin  liquid 
film  between  a  vapor  bubble  and  the  heated  tube  surface, 
thereby  providing  large  heat- transfer  coefficients.  Higher 
heat  fluxes  and  greater  oil  concentrations  tend  to  increase 
the  foaming,  which  could  explain  why  this  enhancement  is  seen 
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more  predominately  at  high  heat  fluxes.  However,  as  oil 
concentration  increases,  eventually  the  oil  concentration 
gradient  becomes  so  steep  that  mechanisms  (1)  and  (2)  tend  to 
dominate  and  heat  transfer  performance  drops  off  (ie.  the 
foaming  enhancement  effect  is  offset) . 

At  low  heat  fluxes,  the  scatter  in  the  data  is  attributed 
to  uncertainty.  For  the  increasing  heat  flux  data  (Fig  5.5), 
it  appears  that  the  addition  of  oil  tends  to  delay  the  point 
at  which  boiling  commences.  Due  to  the  statistical  nature  of 
incipience,  however,  one  should  be  wary  of  drawing  much 
conclusion  from  this  figure.  The  effect  of  oil  on  incipience 
is  discussed  in  greater  detail  later  in  this  chapter. 

Figure  5.7  compares  smooth  tube  decreasing  heat  flux  data 
from  this  thesis  for  R-124  with  Sugiyama  [Ref.  2]  R-114  data 
at  oil  concentrations  of  0%,  3%,  and  10%.  At  all  levels  of 
heat  flux,  the  R-124  data  consistently  demonstrate  improvement 
in  heat  transfer  over  R-114.  This  ranges  from  50%  at  low  heat 
flux  to  15%  at  high  heat  flux.  The  heat  flux  from  a  boiling 
surface  is  proportional  to  bubble  diameter  (db)  ,  bubble 
departure  frequency  (fb),  number  of  active  nucleation  sites 
per  unit  area,  vapor  density  (p) ,  and  latent  heat  of 
vaporization  (hf  ) .  Barthau  [Ref.  15]  concluded,  after 
determination  of  db,  fb,  and  active  nucleation  site  density 
for  a  smooth  copper  tube  in  R-114  using  optical  techniques, 
that  as  saturation  pressure  was  increased,  db  decreased,  while 
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fb,  nucleation  site  density  and  heat  transfer  increased.  While 
conducting  the  tests  at  2.2  °C  the  corresponding  saturation 
pressure  of  11  psig  for  R-124  is  greater  than  the  pressure  of 
0  psig  for  R-114.  This  suggests  that  the  nucleation  site 
density  was  greater  for  R-124  than  for  R-114  at  the  same 
saturation  temperature.  Added  to  this  is  the  higher  latent 
heat  of  vaporization,  vapor  density,  and  thermal  conductivity 
for  R-124  at  this  saturation  temperature.  The  R-114/R-124 
property  comparison  at  2.2  °C  saturation  temperature  is 
summarized  in  Table  1.  It  is  therefore  reasonable  to  expect  a 
higher  heat  transfer  coefficient  for  R-124  at  the  tested 
saturation  temperature  of  2.2  °C . 

D.   BOILING  PERFORMANCE  OF  19  AND  2  6  FPI  GEWA-K  TUBES 

Figures  5.8  and  5.9  show  the  heat  transfer  performance  of 
a  19  fpi  GEWA-K  tube  in  pure  R-124,  R-124/3%  oil,  and  R- 
124/10%  oil  during  increasing  and  decreasing  heat  flux 
respectively.  There  is  a  consistent  increase  in  thermal 
performance  for  oil  concentrations  of  3%  and  10%  over  pure  R- 
124  at  high  heat  fluxes.  This  increase  is  substantially 
greater  than  that  obtained  with  the  GEWA-K  26  fpi  tube  for  the 
same  oil  mixtures,  as  shown  in  Figures  5.10  and  5.11  for 
increasing  and  decreasing  flux  respectively.  Surprisingly  for 
both  tubes,  there  seems  to  be  a  significant  influence  of  oil 
in  the  natural   convection   region  (Figs  5.8  &  5.10) .   There 
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Figure  5.8  Performance  Comparison  for  Boiling  R-124/  0%,3%,& 
10%  Mixtures  for  19  fpi  GEWA-K  Tube  (increasing  flux) 
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Figure  5.9  Performance  Comparison  for  Boiling  R-124/0%, 3%, 10% 
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Figure  5.10     Performance  Comparison  for  Boiling  R-124/   0%,3%, 
&    10%   Oil   Mixtures    for   26    fpi    GEWA-K   Tube    (increasing   flux) 
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tends  to  be  a  decrease  in  performance  for  increasing  oil 
concentration.  This  was  not  so  obvious  with  Sugiyama  [Ref.  2] 
R-114  data.  Once  again,  the  point  of  incipience  seems  to  be 
slightly  delayed  with  increasing  oil  percentage.  The  thermal 
performance  of  the  finned  tube  is  better  than  the  smooth  tube, 
evidenced  by  thermal  enhancements  of  2.1  to  2.9  with  the  19 
fpi  tube  and  1.3  to  2.1  with  the  26  fpi  tube.  However,  it 
should  be  noted  that  the  heat  flux  calculation  does  not 
account  for  the  increased  surface  area  of  the  finned  tube,  but 
rather  uses  a  root  diameter.  Some  of  the  enhancement  is 
therefore  due  to  the  increased  surface  area  of  the  fins  (Figs. 
5.24-5.26  compare  heat  transfer  coefficients  with  smooth 
tube)  ,  with  each  fin  behaving  essentially  as  the  copper 
surface  of  a  smooth  tube  in  terms  of  thermal  performance. 

Figures  5.12  and  5.13  compare  R-124  and  R-114  decreasing 
heat  flux  data  for  the  19  and  26  fpi  GEWA-K  tubes  respectively 
for  0%,  3%,  and  10%  oil  concentrations.  The  19  fpi  tube  shows 
considerable  improvement  over  the  entire  heat  flux  range  while 
the  2  6  fpi  tube  shows  improvement  only  in  the  low  heat  flux 
region.  The  improvement  of  the  lower  fin  density  tube  can  be 
speculated  to  be  due  in  large  part  to  better  circulation  and 
mixing  between  the  fins.  For  both  tubes  at  low  oil 
concentrations  and  high  heat  flux,  the  foaming  effect  provides 
improvement  in  heat  transfer  similar  to  that  seen  with  the 
smooth  tube. 
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E.   BOILING  PERFORMANCE  OF  HIGH  FLUX  AND  TURBO -B  TUBES 

Figures  5.14  and  5.15  show  the  heat  transfer  performance 
of  a  HIGH  FLUX  tube  in  pure  R-124,  R-124/3%  oil,  and  R-124/10% 
oil  during  increasing  and  decreasing  heat  flux  respectively. 
Uncertainty  analysis  (Appendix  D)  was  performed  on  the  HIGH 
FLUX  tube  and  typical  uncertainty  bands  are  also  given.  The 
uncertainty  is  largest  at  low  wall  superheat  and  low  heat  flux 
since  the  uncertainty  is  of  the  order  of  the  measured  values. 
It  can  be  seen  that  the  heat  transfer  performance  decreases 
for  all  heat  fluxes  as  oil  concentration  increases.  This 
decrease  in  performance  is  due  most  likely  to  the  oil  becoming 
trapped  within  the  porous  coating.  As  the  refrigerant  is 
evaporated,  a  very  steep  concentration  gradient  is  developed 
within  the  pores;  oil  diffusion  is  also  retarded  because  of 
the  presence  of  the  interconnected  pores;  this  is  particularly 
severe  at  high  heat  fluxes  with  10%  oil  mixture,  where  a  very 
large  drop-off  in  performance  is  seen.  In  the  natural 
convection  region,  the  effect  of  oil  is  not  as  pronounced  as 
at  high  flux,  exhibiting  only  minor  reduction  in  heat  transfer 
compared  to  the  pure  refrigerant.  Any  advantageous  effect  of 
foaming  is  negated  due  to  the  pores  remaining  impenetrable  to 
the  foam.  Figure  5.16  compares  R-124  and  R-114  decreasing  heat 
flux  data  for  the  HIGH  FLUX  tube  for  0%,3%,  and  10%  oil 
concentrations.  The  trend  of  the  data  for  both  fluids  is  very 
similar.  With  such  a  high  density  of  nucleation  sites  already 
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Figure  5.14   Performance  Comparison  for  Boiling  R-124/0%,3% 
&  10%  Oil  Mixtures  for  HIGH  FLUX  Tube  (increasing  flux) 
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Figure  5.15   Performance  Comparison  for  Boiling  R- 124/0%, 3% 
&  10%  Oil  Mixtures  for  HIGH  FLUX  Tube  (decreasing  flux) 
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present  on  the  surface,  there  does  not  appear  to  be  a 
significant  thermal  enhancement  at  higher  operating  pressures 
as  seen  with  the  finned  and  smooth  tubes. 

Figures  5.17  and  5.18  show  the  heat  transfer  performance 
of  a  TURBO-B  tube  in  pure  R-124,  R-124/3%  oil,  and  R-124/10% 
oil  mixtures  during  increasing  and  decreasing  heat  flux 
respectively.  The  results  of  oil  addition  are  similar  to  the 
HIGH  FLUX  tube.  The  notable  exception  is  that  at  high  oil 
concentrations  and  high  heat  fluxes,  the  heat  transfer 
performance  of  the  TURBO-B  tube  does  not  drop  off  as  rapidly 
as  the  HIGH  FLUX  tube.  Again,  only  a  small  effect  of  oil  in 
the  natural  convection  region  was  observed. 

Figure  5.19  compares  R-124  and  R-114  decreasing  heat  flux 
data   for   the   TURBO-B   tube   for   0%,   3%,   and   10%   oil 
concentrations.   Again,   this  comparison  is  similar  to  the 
results  already  discussed  with  the  HIGH  FLUX  tube. 

F.   EFFECT  OF  R-124 /OIL  MIXTURES  ON  INCIPIENCE  OF  NUCLEATE 
BOILING 

Separate  tests  were  conducted  to  analyze  the  incipience  of 

nucleate  boiling  or  the  onset  of  nucleate  boiling  (ONB)  in 

pure  R-124  and  how  this  was  affected  by  the  presence  of  oil. 

These   tests  were   taken  in  accordance  with  the  general 

procedures  outlined  in  the  preceeding  chapter.   However, 

smaller  heat  flux  increments  were  applied  to  the  test  tube 

during  the  natural  convection  region  in  order  to  obtain  the 

most  accurate  wall  superheat  at  the  ONB.   Also,  during  these 
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Figure  5.17   Performance  Comparison  for  Boiling  R-124/0%,3%, 
&  10%  Oil  Mixtures  for  TURBO-B  Tube  (increasing  flux) 
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Figure    5.18      Performance   Comparison   for   Boiling   R-124/0%,3%, 
&  10%  Oil  Mixtures  for  TURBO-B  Tube    (decreasing  flux) 


66 


(M      M      (M 

iii 

C     C     I 


1  5..2T 

H  ST    3 

0>  3    CO 

3  « 
CO 


# 


u  - 

a.    U 


■a 
o 


§8 


£   5 


<">•   -   — 


y  y 

X    a. 

U    CJ 


I     ♦    o    o 


I  I  1   I    I     I      I         I  I  I  I  I   I    I     I      I         I  1 1  I   I   I    I     I      I         I 


s 

© 
© 
o 


o 
o 
© 
o 


< 

E 


■a 

a 
1 


m 
v 

X 


© 
© 
© 


11 

u 

3 

■j 

X 

3 

5 

H 

ca 

• 

oc 

3> 

CO 

cN 

C£ 

i 

3 

U  b 

u* 

i« 

U 

s2 

— • 

Cfl 

13 

cy 

— 

s 

3 

* 

Tf 

.2 

cN 

2 

u 

S3 

— 

O 

u 

^ 

^m 

i^ 

■-m 

fH 

o 

u 

9* 

— 

o 

u 

t/j 

ha 

T3 

#" 

<3 

a 

« 

s 

*9< 

0 

•M 

U 

i 

u 

H 

U 

• 

o 

m 

** 

© 
© 


(M-rvW/M)/1u»P!JJ»°D    J9jsu«Ji-iB»H 


67 


incipience  tests  the  refrigerant  pool  was  cooled  down  at 
approximately  half  the  rate  of  the  prior  increasing  and 
decreasing  heat  flux  runs.  Each  incipience  test  was  conducted 
with  a  minimum  of  8  hours  between  between  re- testing  to  allow 
adequate  time  for  the  R-124  pool  to  settle  and  for  the 
nucleation  sites  to  de-activate.  The  maximum  time  allowed  for 
the  aging  of  the  tube  surface  was  48  hours.  Similar  tests  were 
repeated  seven  times  on  the  smooth  and  HIGH  FLUX  tubes.  Figure 
5.20  summarizes  the  data  on  a  probability  chart.  Due  to  the 
statistical  nature  of  ONB,  You  et  al  [Ref.  16]  proposed  this 
chart  as  the  best  means  of  presenting  incipient  data.  The 
seven  incipient  wall  superheats  for  each  condition  are  given 
for  the  smooth  and  HIGH  FLUX  tube  at  0%,  3%,  and  10%  oil 
concentrations  as  the  percent  probability  of  nucleation 
occurring.  For  example,  the  smooth  tube  in  pure  refrigerant 
never  nucleated  at  a  wall  superheat  less  than  12°C.  Also  the 
data  show  that  for  a  smooth  tube  in  pure  R-124,  nucleation 
occurred  four  times  at  a  wall  superheat  greater  than  20  °C, 
and  three  times  at  a  lower  wall  superheat;  this  corresponds  to 
a  probability  of  nucleation  occurring  at  a  superheat  of  20  °C 
between  3/7  (43)%  and  4/7  (57)%.  As  the  superheat  further 
increases,  the  probability  increases  also  until  at  a  superheat 
of  28°C  there  is  a  100%  probability  of  nucleation.  Bar-Cohen 
[Ref.  3]  also  reported  widely  varying  R-113  superheats  when 
'departure   from   convection'  or   nucleation   occurred.   The 
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Figure  5.20   Probability  of  Nucleation  Occurring  with  Smooth 
and  HIGH  FLUX  Tube  for  0%,3%,  &  10%  Oil  Mixture 
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probability  chart  of  Figure  5.20  can  be  made  more  accurate  by 
plotting  a  greater  number  of  incipient  superheats  for  each 
tube . 

For  the  smooth  tube,  there  is  a  distinct  shift  to  higher 
superheats  with  increases  in  oil  concentration.  This 
represents  a  delay  in  incipience  caused  by  the  addition  of 
oil.  Memory  and  Marto  [Ref.  10]  conducted  similar  tests  with 
R-114  and  R-114/oil  mixtures  and  reported  the  oil 
concentration  to  have  no  systematic  effect  on  the  point  of 
incipient  boiling  for  the  smooth  tube.  Wanniarachchi  et  al 
[Ref.  9]  reported  an  increase  in  the  wall  superheat  needed  to 
initiate  boiling  for  the  smooth  tube  in  R-114  with  10%  oil, 
noting  that  3%  oil  caused  no  increase.  Memory  and  Marto  [Ref. 
10]  precicted  the  effect  of  oil  on  superheat  needed  to 
nucleate  from  a  vapor- filled  cavity  by  incorporating  the 
physical  properties  of  the  refrigerant/oil  mixture  into  a 
nucleation  parameter  (N)  given  by  Marto  and  Lepere  [Ref.  17] : 

N=a,Tsat/^vapor'hfg 

If  the  parameter  were  to  increase,  a  delay  in  the  incipient 
point  with  increasing  oil  concentrations  would  appear  likely. 
However,  it  is  not  clear  how  the  oil  effects  this  parameter; 
with  the  smooth  tube  in  R-124/oil  mixtures  it  resulted  in  a 
delay  in  incipience,  while  the  HIGH  FLUX  tube  exhibited  no 
distinct  delay  in  incipience.  Also  the  HIGH  FLUX  tube  has  a 
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whole  range  of  sites  that  can  be  easily  activated  while  the 
smooth  tube  only  tends  to  have  very  small  sites  which  may  move 
outside  the  range  of  sites  that  can  be  activated  when  a 
refrigerant/oil  mixture  is  used. 

In  the  case  of  pure  refrigerant,  it  does  appear  that  with 
the  relatively  few  incipient  tests  conducted,  the  wall 
superheats  at  incipience  for  the  smooth  and  HIGH  FLUX  tubes 
are  larger  on  average  for  R-114  than  for  R-124.  This  initial 
result  would  support  the  approximately  50%  greater  nucleation 
parameter  of  R-114  over  R-124.  This  subject  requires  further 
study. 

Figure  5.21  shows  an  expanded  view  of  the  smooth  tube's 
heater  sleeve  rolled  out  to  show  the  exact  locations  of  each 
thermocouple  (longitudinal  and  circumferential)  imbedded  on 
the  surface  of  the  sleeve.  Several  data  sets  were  taken  for 
the  smooth  tube  and  then  repeated  after  reversing  the  position 
of  the  tube  within  the  pool.  This  procedure  was  conducted  to 
see  if  there  were  any  noticeable  trends  in  the  thermocouple 
readings  which  may  be  an  indication  of  local  hot  or  cold 
regions  within  the  pool.  However,  this  reversal  of  the  tube 
had  no  consistent  effect  on  the  thermocouple  readings  or  the 
point  at  which  nucleation  occurred.  This  appears  to  underline 
the  uniformity  of  thermocouple  output  within  the  pool.  A  more 
detailed  survey  of  this  data  needs  to  be  done. 
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Figure  5.21   Expanded  View  of  Thermocouple  Locations  within 
the  Heater  Sleeve  of  Smooth  Tube  (#1  at  12  O'clock) 
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G.   SUMMARY  OP  THE  R-124  DATA 

Figures  5.22  and  5.23  summarize  the  effect  of  all  oil 
concentrations  ( 0 , 1 , 2 , 3 , 6 , &10%)  in  R-124  for  the  five  tubes  at 
heat  fluxes  of  25  and  90  KW/m2  respectively;  25  kW/m2 
represents  a  typical  operating  heat  flux.  For  each  tube,  the 
ratio  of  the  heat  -  transfer  coefficient  at  a  given  oil 
concentration  divided  by  that  for  the  same  tube  in  pure  R-124 
is  given.  At  low  heat  flux,  the  smooth  tube  and  the  19  fpi 
tube  give  an  increase  in  performance  as  oil  is  added,  up  to 
concentrations  of  6%.  The  smooth  tube  does  the  best  in  this 
respect,  giving  an  increase  in  heat  transfer  of  about  23% 
(compared  to  a  smooth  tube  in  pure  R-124)  at  6%  oil 
concentration.  The  HIGH  FLUX  and  TURBO- B  tubes  show  a  decrease 
in  performance  as  oil  is  added,  while  the  26  fpi  tube  shows  no 
change.  At  the  highest  oil  concentration  (10%),  performance 
for  all  tubes  is  worse  than  that  at  6%,  although  the  smooth 
tube  still  indicates  better  heat  transfer  than  when  no  oil  was 
added,  which  is  rather  surprising.  Similar  trends  were 
obtained  for  R-114,  although  at  10%  oil,  all  tubes  had  poorer 
heat  transfer  than  with  no  oil.  At  higher  heat  flux,  general 
effects  are  similar,  although  now  the  19  fpi  tube  shows  the 
best  performance  of  about  35%  at  6%  oil  concentration. 
Interestingly,  between  3  and  6%  oil  concentrations,  all  tubes 
except  the  HIGH  FLUX  tube  show  improved  heat  transfer  over  the 
same  tube  in  pure  refrigerant.  Although  below  3%  oil  there  is 
an  unexplained  'dip'  in  the  heat  transfer  coefficient  ratio 
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with  the  TURBO -B  tube.  The  HIGH  FLUX  surface  has  very  poor 
performance  at  10%  oil  concentration,  being  only  35%  of  that 
in  the  pure  fluid. 

Figures  5.24,  5.25  and  5.26  compare  the  five  tubes 
directly  for  oil  concentrations  of  0,  3,  and  10%  respectively. 
From  these  figures,  enhancements  (ie.  improvement  in  heat 
transfer  compared  to  the  smooth  tube)  can  be  calculated  at  any 
desired  heat  flux.  These  are  provided  in  Table  4  at  heat 
fluxes  2  5  and  9  0  kW/m2 .  The  advantages  of  enhanced  surfaces 
are  clearly  apparent.  The  smooth  tube  always  has  the  poorest 
performance  (as  expected)  and  is  used  as  a  base-line  with 
which  the  other  tubes  are  compared.  Between  the  finned  tubes, 
the  19  fpi  tube  consistently  outperforms  the  26  fpi  tube  at 
all  heat  fluxes  and  oil  concentrations:  enhancements  range 
between  2.1-2.9  for  the  19  fpi  tube  and  1.3-2.1  for  the  26  fpi 
tube.  Indeed,  at  the  highest  heat  fluxes  at  10%  oil 
concentration,  the  19  fpi  tube  is  the  best  tube  tested 
(enhancement  of  2.6).  The  GEWA-K  19  fpi  tube  to  smooth  tube 
area  ratio  is  2.4,  while  the  GEWA-K  26  fpi  tube  to  smooth  tube 
area  ratio  is  3.2.  The  HIGH  FLUX  and  TURBO -B  tubes  are  the 
best  tubes  over  most  conditions  of  heat  flux  and  oil 
concentration,  with  both  having  fairly  equal  performance.  The 
exception  is  the  HIGH  FLUX  tube  at  high  heat  flux  and  oil 
concentration,  where  performance  drops  off  rapidly,  as  already 
mentioned.  Over  practical  heat  flux  ranges  (15-30  kW/m2) 
therefore,  the  porous  coated  and  modified  fin  tubes  are  the 
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best  for  all  oil  concentrations.  Although  the  TURBO-B  and  HIGH 
FLUX  tube  performance  is  degraded  at  10%  oil  concentration, 
in-service  oil  separators  should  maintain  a  3%  or  less 
concentration  in  an  actual  evaporator  system.  At  this  level, 
the  TURBO-B  and  HIGH  FLUX  tube  are  superior  at  normal  chiller 
heat  flux  between  10  and  30  KW/m2 . 

A  point  of  interest  worth  noting  was  that  the  Viton  and 
neoprene  0- rings  used  experienced  considerable  swelling  of 
approximately  30%  by  volume  when  subjected  to  R-124.  This  was 
enough  to  prevent  their  re-use  when  the  Teflon  inserts  were 
removed  to  change  the  tube.  Also  noted  was  that  the  0- rings 
only  swelled  when  taken  out  of  the  system.  When  undisturbed, 
they  maintained  a  good  seal  (this  might  not  be  true  if  the 
tube  were  sliding  past  the  0-ring) .  Also  significant,  after 
removed  from  use  within  the  apparatus  for  60  days,  both  types 
of  0- rings  did  shrink  back  to  their  original  size. 

Table  4.   ENHANCEMENT  SUMMARY  FOR  R-124 


Heat   Flux 

2  5  kW/m2 

9  0  kW/m2 

Oil  Concentration 

0% 

3% 

10% 

0% 

3% 

10% 

Smooth 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

19  fpi 

2.9 

2.7 

2.8 

2.1 

2.4 

2  .6 

26  fpi 

2.1 

1.7 

1.8 

1.3 

1.3 

1.5 

HIGH  FLUX 

6.2 

4.9 

4.8 

2.9 

2  .3 

1  .  1 

TURBO  -  B 

6.0 

4.3 

3.7 

2.4 

2.2 

1.9 
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VI.   CONCLUSIONS 


1.  For  the  smooth  and  finned  tubes,  significant  improvements 
in  heat  transfer  (>50%)  were  obtained  in  pure  R-124 
compared  with  pure  R-114  at  the  same  saturation 
temperature.  This  was  attributed  to  the  higher  saturation 
pressure  of  R-124  activating  more  nucleation  sites  and  the 
higher  latent  heat  of  vaporization,  vapor  density,  and 
thermal  conductivity  of  R-124  over  R-114. 

2.  For  the  HIGH  FLUX  and  TURBO-B  surfaces,  performance  in  the 
two  refrigerants  was  similar,  since  the  number  of  active 
nucleation  sites  is  already  high. 

3.  For  pure  R-124,  the  HIGH  FLUX  tube  gave  the  largest 
enhancement  of  any  tube  over  the  whole  range  of  heat  flux. 

4.  For  the  smooth  and  finned  tubes,  significant  improvements 
in  the  heat  transfer  (>50%)  were  also  obtained  in  R- 
124/oil  mixtures  when  compared  with  R-114/oil  mixtures  at 
the  same  saturation  temperature.  This  is  also  assumed  to 
be  due  to  the  higher  saturation  pressures  associated  with 
R-124. 

5.  Increasing  the  oil  concentration  in  R-124  by  up  to  6%  for 
the  smooth  and  finned  tubes  improved  heat  transfer 
performance  compared  to  the  pure  refrigerant.  For  example 
at  a  heat  flux  of  at  90  kW/irr ,  performance  for  the  smooth 
and  finned  tubes  increased  by  18%  and  35%  respectively.  At 
higher  oil  concentrations,  performance  deteriorated. 

6.  Any  increase  in  oil  concentration  in  R-124  for  the  HIGH 
FLUX  or  TURBO-B  tubes  led  to  a  decrease  in  performance 
such  that  at  high  heat  fluxes  and  high  oil  concentrations, 
the  HIGH  FLUX  tube  gave  only  35%  of  its  performance  in 
pure  R-124  and  the  19  fpi  tube  gave  the  best  overall 
performance . 

7.  At  chiller  design  heat  fluxes  (10-30  kW/m2)  and  oil 
concentrations  of  3%  maximum,  the  HIGH  FLUX  and  TURBO-B 
surfaces  still  provide  the  best  heat  transfer  which  is 
about  the  same  as  with  R-114.  Therefore,  similar  or 
slightly  better  performance  is  expected  with  R-124  when 
used  as  a  'drop- in'  replacement  for  R-114  in  a  refrigerant 
evaporator  operated  at  the  same  saturation  temperature. 

8.  Addition  of  oil  delays  the  point  of  incipience  (occurs  at 
a  higher  wall  superheat)  for  the  smooth  tube.  For  the  HIGH 
FLUX  tube,  oil  addition  has  no  measurable  effect  on  the 
wall  superheat  at  incipience. 
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VII.   RECOMMENDATIONS 

1.  The  remaining  enhanced  tubes  should  be  tested  in  R-124  at 
0%,3%,  and  10%  oil  concentrations  and  compared  to  the  five 
tubes  above  and  their  prior  performance  in  R-114  and  R- 
114/oil  mixtures. 

2.  Incipience  tests  performed  with  finned  tubes  and  TURBO-B 
tubes  should  be  conducted  to  compare  the  influence  of  oil 
on  wall  superheat. 

3  .  The  influence  of  the  variation  in  time  between  test  runs 
on  boiling  nucleation  (ie.  onset  of  nucleate  boiling)  must 
be  investigated  more  thoroughly. 

4.  The  single  tube  apparatus  should  be  modified  to  test  an 
array  of  two  horizontal  boiling  tubes  and  compare  these 
results  with  the  R-114  tube  array  data  previously 
recorded. 

5.  A  prototype  chiller  presently  using  R-114  should  be  tested 
with  R-124  to  monitor  the  performance. 
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APPENDIX  A.   THERMOPHYSICAL  PROPERTIES  OF  R-124 

The  following  thermophysical  properties  of  saturated  R-124 
are  plotted  versus  temperature  in  degrees  Celsius.  They  were 
generated  from  REFPROP  [Ref .  12]  . 


1. 

Saturation  Pressure   / (kPa) 

2  . 

Liquid  Density   /(kg/m3) 

3  . 

Latent  Heat  of  Vaporization   / (kJ/kg) 

4. 

Liquid  Specific  Heat   / (kJ/kg  K) 

5. 

Liquid  Viscosity   / (mP)    note:  1  mP=  lxlO"4 

kg/m  s 

6. 

Liquid  Thermal  Conductivity   / (W/m  K) 
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APPENDIX  B.   REPRESENTATIVE  DATA  SET 


Oate  :     19  Nov  1992 

NOTE:  Program  name  :  ORPSB 
Ois*  number   *  01 

Oate  :     19  Nov  1992 

NOTE:  Program  name  :  0RPG8 
Disk  number   -  03 
New  file  name:  Oil  191 93 
TC  are  defective  at  locations   7      8 
Tube  Number:    13 

Oata  Set  Number  -    I   Bulk  Oil  Z  «   3.0 
TIME:      15:40:01 

TC  No:     1       2       3       4       5      6       7      8 
Temp  :    3.23    3.38    3.31    3.20    3 . 05    3.32  -99.99  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 
3.25    2.14    2.14    2.18   10.84   -13.3 
Thetab   Htube      Qdo 
1.063   6.475E*02   6.879E*02 

Oata  Set  Number  -    2   Bulk  Oil  X  »   3.0 
TIME:      15:42:22 

TC  No:     I       2       3       4       5      S       7       8 

Temp  :  3.58    3.80    3.72    3. SI    3.45   3.54  -99.99  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

3.53  2.22    2.24    2.25   10.91   -13.2 

Thetab  Htube      Qdp 

1.377  6.422E<-02   8.34SE+02 

Oata  Set  Number  •    3   Bulk  Oil  Z  •   3.0 
TIME:       15:45:20 

TC  No:     I       2       3       4       5      5       7      8 

Temp  :  3.8S    4.35    4. 35    4.19    3.37   3.95  -99.99  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

4.09  2.09    2.21     2.17   10.83   -13.1 

Thetab  Htube      Qdp 

1.921  6.S23E»02   I.2S3E+03 

Oata  Set  Number  «    4   8ulk  Oil  Z  -   3.0 

TIME:  15:47:12 

TC  No:  1       2      3      4      5      6      7      8 

Temp  :  4.10   5.03    5.02    4.87    4.45    4.23  -99.99  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

4.52  2.15    2.27    2.27   10.92   -13.0 

Thetab  Htube      Qdp 

2.350  5.986E+02   I.S42E+03 
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Data  Set  Number  -    5   8ulk  Oil  X  • 

TIME:  15:49:09 

TC  No:  t       2       3       4 

Temp  :  4.30    5.75   5 . 7S    5.44 
Twa     Tliqd    T 1  i  q  d2   Tvapr   Psat 

5.07  :.0S    2.20    2.23   10. 3$ 

Thetab  Htube      Qdp 

2.345  7.332E+02   2.22SE+03 


3.0 

5      6      7 
4.58    4.5b  -99.99 
Tsump 


3 
99.99 


3ulk  Oil  x  =   3.0 


Cata  Set  Number  ■    5 
TIME:      15:51 : 30 
TC  No:     1       2       3       4 

Temp  :  4.40    6.62    6.60    6.21 
Twa     Tliqd    Tliqd2   Tvapr   Psat 

5.61  2.12     2.16    2.28   10.93   -12. 

Thetab  Htube      Qdp 

3.331  9.I54E+02   3.049E*03 


5      6      7      8 
5.12    4.73  -99.99  -99.99 
Tsump 


Data  Set  Number  -    7   Bulk  Oil  X  =■   3.0 
TIME:      15:54:01 

TC  No:     1       2       3       4       5       6 

Temp  :  4.75    7.39    7.27   6.75   5.39    B.2J 
Twa    Tliqd    Tliqd2   Tvapr   Psat    Tsunp 

6.11  2.15    2.21     2.23   10.93   -12.7 

Thetab  Htube      Qdp 

3.831  1.031E+03   3.951E+03 


7 
■99.99 


■99.99 


Oata  Set  Number  -    8   Bulk  Oil  X  =   3.0 

TIME:  15:56:43 

TC  No:  1       2       3       4       5       6       7 

Temp  :  5.14    8.28   8.15   7.38   5.55    5.59  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

5.56  2.24    2.27    2.35   11.00   -12.6 

Thetab  Htube      Qdp 

4.302  I.2S9E+03   S.413E+03 


•99.99 


Oata  Set  Number  -    9   Bulk  Oil  X  -   3.0 
TIME:      15:59:50 

TC  No:     I       2       3       4       5       6       7 

Temp  :  4.78    6.39    8.17    6.85    4.88    4.45  -99.99 
Twa    Tliqd    Tliqd2   Tvapr   Psat    Tsunp 

5.97  2.29    2.24    2.32   10.97   -12.5 

Thetab  Htube      Qdp 

3.650  Z.1SIE+03   7.849E*03 


■99.99 


Oata  Set  Number  -   10   Bulk  Oil  X  »   3.0 
TIME:      16:02:33 

TC  No:     I       2       3       4       5       6       7       3 

Temp  :  4.83   5.18    5.13    4.84    4.91    4.53  -99.99  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

4.86  2.30    2.25    2.23   10.93   -12.4 

Thetab  Htube      Qdp 

2.533  4.150E+03   1.072E+04 

Oata  Set  Number  -   II   Bulk  Oil  X  -   3.0 
TIME:      16:17:38 

TC  No:     I       2       3       4       5       6       7       8 

Temp  :  5.14    5.55    5.59   5.22    5.25    4.77  -39.99  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsunp 

5.20  2.20    2.12    2.15   10.82   -11.9 

Thetab  Htube      Qdp 

3.044  5.093E+03   I.552E+04 
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Oat*  Set  Nu«oer  »   12   Bulk  Oil  X  -   3.3 
TIME:       16:23:31 

rc  No:     1      2      3      4      5      5      7      3 

Temp  :  5. "9   5.34   6.23    5.90    6.31    5.24  -39.99  -99.99 

Tua  Tliqa    Tliqa2   Tvapr   Psat    Tsump 

5.83  2.26    2.19    2.23   10.96   -M.9 

Thetab  Htube      Qdp 

3.633  6.3296*33   2.2976*34 

Oata  Set  Nunoer  •   13   Bulk  Oil  X  -   3.0 
TIME:      16:23:15 

TC  No:     I       2       3       4       5       6      7      g 

Temp  :  6.51    7.18    7.13    6.69    6.85    5.99  -99.99  -99.99 

Tua  Tliqd    Tliqd2   Tvapr   Psat    Tsump 

6.60  2.31     2.26    2.23   10.93   -11.8 

Thetab  Htube      Qdp 

4.316  7.794E+33   3.364E*04 

Oata  Set  Number  -   14   Bulk  Oil  X  -   3.0 
TIME:       16:24:49 

TC  No:     I       2       3       4       5       6       7       9 

Temp  :  7.40    7.84    7.85    7.36    7.58    6.49  -99.99  -99.99 
Tua     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

7.22  2.35    2.27    2.32   13.97   -11.3 

Thetab  Htube      Qdp 

4.895  9.726E*03   4.76lE*04 

Oata  Set  Number  -   15   Bulk  Oil  X  »   3.0 
TIME:      16:29:15 

TC  No:     I       2       3       4       5       6       7      8 

Temp  :  8.17    8.32    8.52    8.M    8.32    7.09  -99.99  -99.99 
Tua     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

7.78  2.29    2.23    2.28   10.93   -11.7 

Thetab  Htube      Qdp 

5.501  i.300E*04   7.ISIE*04 

Oata  Set  Number  -   16   Bulk  Oil  X  -   3.0 
TIME:       16:33:55 

TC  No:  1       2       3       4       5       6       7      8 

Temp  :  8.79   8.96    9.33    8.97    9.12    7.59  -99.99  -99.99 
Tua     Tliqd    TUqd2   Tvapr   Psat    Tsump 

8.37  2.20    2.14    2.20   10.86   -11.4 

Thetab  Htube      Qdp 

6.169  1.5S3E*04   9.5776*04 

Oata  Set  Number  -   17   Bulk  Oil  X  •   3.0 
TIME:      16:36:37 

TC  No:  I       2       3       4       5       6       7      8 

Temp  :  7.44    7.74    7.92    7.59    7.77    6.56  -99.99  -99.99 
Tua     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

7.24  2.18    2.M     2.17   10.83   -1  1  .4 

Thetab  Htube      Qdp 

5.070  l.230E*04   6.235E+04 

Oata  Set  Number  -   18   Bulk  Oil  X  -   3.0 
TIME:       16:40:03 

TC  No:     I       2       3       4       5       6       7      8 

Temp  :  6.63   7.33    7.12    6.77    6.99    5.96  -99.99  -99.99 
Tua     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

6.57  2.24    2.18    2.24   10.90   -1  1  .4 

Thetab  Htube      Qdp 

4.329  9.65~E*03   4.180E*34 
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Data  Set  Number  ■   '9   Bulk  Oil  X    •   2.8 
TI"E:       16:43:34 

TC  Sc:     :      Z      3      *      5      5      7 

Tei"P  :  6.09   6.65    5.56    B.27    6.41    5.54  -99.99 
Twa     Tliqd    TiiqdZ   T.ap"-   Put    Tjuwp 

6.14  2.27     2.20     2.24   10.90   -II. S 

Thetab  Htube       Qdp 

3.699  7.8596*03   3.06~E*04 


a 

•99.99 


Cata  Set  Number  =   20   Bulk  Oil  X    =   3.0 

TIME:  16:45:59 

TC  No:  1       2       3       4       5       6       7       8 

Tenp  :  5.62    6.13   6.13    5.32    5.92    5.19  -99.99  -99.99 
Tua     Tliqd    Tliqd2   Tvapr   Psat    T5unp 

5.70  2.23    2.14     2.19   10.85   -11.5 

Thetab  Htube       Qdp 

3.506  6.754E+03   2.368E+04 


Data  Set  Number  -   21   Bulk  Oil  X    =   3.0 

TIME:  16:48:34 

TC  No:  1       2       3       4       5       6       7 

Temp  :  5.17    5.73    5.71    5.36    5.40    4.85  -99.99 
Twa     Tliqd    TiiqdZ   Tvapr   Psat    Tsump 

5.30  2.31     2.22     2.25   10.91   -11.7 

Thetab  Htube      Qdp 

3.045  5.S5SE+03   1.692E+04 


8 
■99.99 


Data  Set  Number  •   22   Bulk  Oil  X  •   3.0 

TIME:  16:51:59 

TC  No:  1       2       3       4       5       6       7 

Temp  :  4.69   5.20   5.15    4 . 3S    4.92    4.44  -99.99 
Twa     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

4.82  2.27    2.17    2.21   10.87   -11.8 

Thetab  Htube       Odp 

2.613  4.525E+03   1.I82E*04 


8 

■99.99 


Oata  Set  Number  ■   23   Bulk  Oil  X   •      3.0 
TIME:      16:54:33 

TC  No:     1       2       3       4       5       6       7 

Temp  :  4.32    4.77    4.63    4.45    4.55    4.21  -99.99 
T*a     Tliqd    TiiqdZ   Tvapr   Psat    Tsump 

4.45  2.30    2.21     2.27   10.92   -11.9 

Thetab  Htube      Qdp 

2.186  3. 6536*03   7.985E*03 


•99.99 


Oata  Set  Number  -   24   Sulk  Oil  t  -   3.3 
TIME:      16:S7:5S 

TC  No:     I       2       3       4       5       6       7 

Temp  :  3.94    4.31    4.23    4.10    4.11    3.91  -99.99 
Tuia     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

4.08  2.24    2.15    2.21   10.87   -12.0 

Thetab  Htube      Qdp 

1.863  2.839E+03   5.289E+03 


■99.99 


Oata  Set  Number  -   25   Bulk  Oil  X  -   3.0 
TIME:      17:01:00 

TC  No:     I       2       3       4       5       6       7 

Temp  :  3.56    3.91    3.94    3.75    3.74    3.60  -99.99 
Tua     Tliqd    Tliqd2   Tvapr   Psat    Tsump 

3.72  2.22     2.08     2.14   10.80   -12.1 

Thetab  Htube       Qdp 

1.577  2.16SE+03   3.4I9E+03 


■99.99 
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0«t«   Set    Nu*t>*r     -      25      Su .  ■    0.1   Z   ■■      3-9 
TIflC:  17:93:15 

TC    Mo:  I  2  3  4  S  S  7  8 

Te-o    :  3. 56        3.83        3.77       3.74        3.69        3. 61    -99.99    -99.99 
T««           i:  .3d         TUgoZ      '-«ur      DJat  'i-"0 

3.69  2.35  2.21  2.27      13.92      -  i  2 . 2 

Thetab  HtuOe  Qap 

I .429  1 .899E*83      2.S"8E'93 

Cat*    Set    Number    -      27      Bulk    Oil    X    •      3.9 
TI«£;  I7:96:9S 

TC   No;  I  2  3  4  5  6  7  3 

Te«p    :  3.41         3.63        3.56        3.52        3.48        3.4S    -99.99    -99.99 
Tw«          Tliqd         Tliqd2      '-tor      Psat         Tju«p 

3.59  2.35  2.29  2.24      19.89      -12.3 

Thetab  Htube  Odp 

1.269  l.S7lE*93      1. 9816*93 

Oete  Set   Number   •     28     Sulk   Oil   t  •     3.9 
TIHE:  17:98:94 

TCNo:l  2  3  4  5  6  7  8 

Tanp    :        3.24        3.42        3.34        3.34        3.39        3.32    -99.99    -99.99 
Twa  Tliqd        Tl»qd2      Tvapr     P»at        Tju«o 

3.32       2.26  2.15         2.14      19.89      -12.4 

Thetab      Htube  Qdo 

1. 185      i.iS9£*93      l.373E*93 

Oat*   Sat    Number  •     29  Bulk   Oil    X   •  3.9 
TIHE:              I7:I9:SS 

TC   No:           1               2  3  4  5              6               7               8 

Te«p    :  3. IS  3.23  3.18       3. 16  3.12        3.29   -99.99   -99.99 

*»«          Tliqd  Tliqd2  Tvapr     P»at  Tsu«p 

3.17  2.28          2.2S  2.23      19.89  -I2.S 
Theteb  Htube               Odp 

.937      t.9S6E*93      9.893E*92 

Oat*   Set    Number    -      39     Bulk   Oil   t   •      3.9 
TIHE:  17:13:34 

TC    Mo'-  1  2  3  4  5  6  7  9 

Te»o    :  3.29        3.49        3.36        3.33        3.29        3.28    -99.99    -99.99 
T«a          Tliqd        TUqd2      Tvaor     Ptat        Tiunp 

3.39  2.21  2.98         Z.H     19.77      -12.6 

Thetab  Htube  Qdo 

1.295  l.246E»93      i.S8i£*93 

Oat*   Set   Nu«ber  •      31  Bulk   Oil    X   •  3.9 
TIHE:               17:16=99 

TC   No:           I  2  3  4              5              6               7               8 

Te«o    =  3. IS  3.2S  3.29       3.16  3.13        3.23    -99.99   -99.99 

T»e         Tliqd  TUqd2  T»epr     P»at        T»u«p 

3.18  2.39  2.27  2.28      18.93      -i2.7 
Thetab  Htube               Odo 

.997     9.9226*92      9.99i£*92 

Oata  Set  Number  •   32  Bulk  Oil  X  •   3.9 

TIME:       17:17:39 

TC  No.     I      2      3      4      5      6      7      8 

Te«p  :  2.92    3.99   2.92   2.89   2.91    3.93  -99.99  -99.99 
Twa    Tiiqq    Tiiqd2   Tvapr  P»et    Tsunp 

2.94  2. IS    2.17    2.15   19.81   -12.7 

Thetab  Htube       Odp 

.796  8. 8886*92   7.9746*92 

Oata  Set  Nu«ber  •   33   Bulk  Oil  X  •   3.8 

TIH£:  17:  19:  16 

TC  No:  I       2      3      4      5      6       7       8 

Te»p  :  2.92    2 . 9S    2.88   2.86    2.87    2.99  -93  59  -99.99 
T«a    TUqd    Tiiqd2   Tvapr  Psat    Tju«o 

2.91  2.17    2.19    2.16   19.82   -12.8 

Thetab  Htube      Odp 

.7S3  7.7*se.92   S.829E*82 
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APPENDIX  C.  SAMPLE  CALCULATIONS 

Data  file  number  D1221SM0  was  used  to  conduct  the  sample 
calculations.  The  saturation  temperature  was  2.28  °C  with  a 
single  smooth  tube  heat  flux  of  31290  W/m2 . 

A.  TEST-TUBE  DIMENSIONS 
DQ=. 015875  m 
Di=.0127  m 

D1=. 01245  m 
L=.2032  m 
Lu=.0762  m 

B.  MEASURED  PARAMETERS 

V=141.95  volts 
1=2.34  amps 
Tl=16.84  °C 
T2-17.15  °C 
T3=16.42  °C 
T4=16.22  °C 
T5=15.66  °C 
T6=15.48  °C 
Tsat=2.28  °C 
kc=344  W/m-K 
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C.  OUTER  WALL  TEMPERATURE  OP  THE  BOILING  TUBE 


p=7r-D0  =  7r-  .  015875  m  =.04987  m 


Ac  =  tt(D02-D12)  /4  =  7r(  (  .015875)  2-  (  .0127)  2)  /4 
Ac=7.13xl0~5    m2 


Qh=VI=141.95-2.34=    332.16    W 


T,T    =ETn/n=E  Tn/6=16.3    °C 

a  vy  ^ 


T„o=Tavg-i  Qh-  (ln(D0/D1)/(27r-L-kc)) 

Two=16.3-I  332.16-  ( In  (.  015875/  .  01245  )/(  2tt  •  L  •  kc) 


Two=16-12    °C 


^=Two-Tsat  =  16.12-2..28  =  13.84    °C 


D.   PROPERTIES  OP  R-124  AT  FILM  TEMPERATURE 

The  thermophysical  properties  of  R-124  were  obtained  from 
REFPROP   [Ref .  12]  and  are  shown  in  Appendix  A. 

Tf=(Two+Tsat)/2=(16.12+2.28)/2=  9.2  °C 

/x  =  3.4593-  (.0426)Tf+  (  3  .  9485x10  "4  )  Tf2  -  (4  .  19  3x10  '6  )  Tf3  + 
(2.0709xlO"8)Tf4 
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/z  =  3  .459  3-  (  .042  6)  -9.2+  (  3  .  9485x10  "4  )  ■  (9.2)  2-  (4. 19  3x10  "6)  ■  (9.2) 3+ 

(2 . 0709X10"8)  ■  (9  .2) 4 
fi  =  3  .  098xl0"4    N-s/m2 


p=143  4.8-  (2.8619)Tf  -  ( 6 . 7267x10 " 3 ) Tf2  -  ( 7 . 2852x10 "5  )  Tf3 
p=1434.8-  (2.8619)  -9.2-  ( 6 . 726  7x10 "3  )  •  (9.2) 2-  ( 7 . 2852x10 "5  )  ■  (9.2) 
p=1407.84    kg/m3 


u=/i/p  =  3.098xl0"4/l407.84  =    2.201x10"7        m2/s 


k=7. 5 191x10 "2-  (3.543  6xl0"4)Tf-  ( 1 . 9  545x10 ~7) Tf2  + (3 . 183  5x10 "9  )  Tf3 


k=7.5191xl0"2-  (3.5436xl0~4)  ■  (9.2)  -  (1.9545x10 


(9.2)"+ 


(3.1835xl0~9)  ■  (9.2) 3 
k=7.192xl0"2    W/m-K 


c    =1.0  542+ (2.14  05xlO"3)Tf+ ( 1 . 070  9x10 ~5 ) Tf2  -  ( 6 . 4721x10 ~8 ) Tf3 
( 1. 4324x10 "9)Tf4+ (4.136xlO"11)Tf5 


c   =1.0  542+ (2. 14 05x10 "3)  ■  (9.2)  +  ( 1 . 0709x10 "5)  •  (9.2)2-  (6.4  721 


xlO"8)  •  (9.2) 3-  (1. 4324x10 ~9)  ■  (9.2)4+ (4.136x10 


-11 


9.2) 


c   =1074.74    J/kg-K 


a=k/ (p-c    ) =    7. 192x10 "2/ (14  07.84-1074.74 
a=4.753x!0~8   m2/s 
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0«- (Ap/AT) /p=2.094xl0~3     (1/K) 
Pr=    u/a=4.63 

E.   HEAT -FLUX  CALCULATION 

The  average  natural  convection  heat- transfer  coefficient 
at  the  non-boiling  ends  of  the  test  tube  is  calculated  using 
the  Churchill  and  Chu  [Ref.  4]  correlation: 

For  simplification  let  1  represent  the  following  in  the 
correlation  : 

R=  [  (g-0-Do3- 6 ■ tanh(m-Lu)  ) / ( v • a • Lu ■ m) ] 1/6 
then 

h=k/D0[.6+.3  87[»/ ( 1+ ( 0 . 559/Pr) 9/16)  8/27]  ]2 

where 

m=[  (h-p)/(kc-Ac)  ]1/2 

and  an  iterative  process  is  used  in  program  DRPGB  to  compute 

the  value  of  h  by  starting  with  the  value  for  h  of  190 

(W/m2,K) .  The  resulting  values  for  h  and  m  were: 

h=265.43  (W/m2-K) 

m=23.23  (1/m) 

therefore : 

Qf= lh-p-kc-Ac) 1/2-  6-  tanh(m  -Lu)=7.44  W 
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F.   HEAT  FLUX  THROUGH  ACTIVE  BOILING  LENGTH 


Q=Qh-2 • Qf= 3 32. 16- (2 ■ 7.44) =317.2  8  W 


Ab=7r-D0'  L=7f0. 015875-  0 . 2032=1 . 0134x10 "2  m2 


q=Q/Ab=317.2  8/1.0134x!0"2=31, 308  W/m' 


h=q/0  =  31, 308/13 .84  =  2262.1  W/m2  •  K 


The   following   results   were   produced   by 
acquisition  and  reduction  program  DRPGB: 
q=31290  W/m2 
0=13.84  °C 
h=2262  W/m2-°K 


the   data 


TABLE  5.   DIMENSIONS  OF  BOILING  TUBES  TESTED 


TUBE 

Dl  mm 

D2  mm 

Di  mm 

Do  mm 

L   mm 

Lu  mm 

W/moK 

smooth 

12.4 

15.9 

12  .7 

15.9 

203.2 

76.2 

344 

GEWA-K 
19  fpi 

10.1 

12.7 

11.1 

12.7 

203.2 

76.2 

344 

GEWA-K 
26  fpi 

10.1 

12  .7 

11.1 

12.7 

203.2 

76.2 

344 

TURBO -B 

11.6 

13.8 

11.8 

15.8 

203  .2 

76.2 

398 

HIGH 
FLUX 

12.9 

15.8 

13  .2 

15.8 

203  .2 

76.2 

45  * 

*   HIGH  FLUX  tube  porous  coating  is  on  a  copper/nickel  base 
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APPENDIX  D.   UNCERTAINTY  ANALYSIS 

Using  methods  of  uncertainty  analysis  outlined  by  Kline 
and  McClintlock  [Ref .  18]  and  performed  by  Sugiyama  [Ref .  2] , 
four  data  points  are  analyzed.  The  data  points  chosen  were  for 
a  Smooth  tube  and  HIGH  FLUX  tube.  For  each  tube,  there  was  a 
data  point  at  a  high  heat  flux  and  one  at  a  low  heat  flux.  The 
following  is  a  sample  uncertainty  analysis  calculation  from 
data  file  number  D1220SM0  with  a  Smooth  tube,  0%  oil,  at  a 
heat  flux  of  86,000  W/m2 .  The  results  of  the  remaining  sample 
points  are  listed  in  Table  6.  Uncertainties  are  given  as  a 
percentage  of  the  calculated  parameter,  i.e.  a  relative 
uncertainty. 

A.   UNCERTAINTY  IN  HEAT - TRANSFER  RATE 

Is=1.99  volts      Vg=9.325  volts 

Estimated  uncertainties  in  the  measured  quantities  I  and  V  : 
51=+  0.025  amps 
<5V=+  0.05  volts 


1=1.9182 • Is=3 .82  amps 

V=25 -Vs=233 . 13  volts        Qh=VI 

where:   6=  uncertainty  in  measurement 
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6Qh/Qh=  ((6V/VS)2+(6I/IS)2)1/2 

5Qh/Qh=  ((.05/9.325)2+(.02  5/1.99)2)1/2 

5Qh/Qh=  1.3  7  % 

B.  UNCERTAINTY  IN  SURFACE  AREA 

Ab=    7T-DO-L 

DQ=    15.88     (mm)  <5DQ=    .1     (mm) 

L=    203.2     (mm)  6L=       .1     (mm) 

6Ab/Ab=     ((6D0/D0)2    +  (6L/L)2)1/2 

<5Ab/Ab=     ((.1/15.88)2  +     (  .  1/203  .  2  )  2  )  1/2 
<5Ab/Ab=    .63    % 

C.  UNCERTAINTY    IN   WALL    SUPERHEAT 


AT=    Two-Tsat  Tsat  =  2.23    °C  6Tsat=     .01    °C 

Two  =    Tavg-Qh  tm  (Do/Di)  /  (2  •  7T  •  L  ■  kc)  ] 
Tn=    thermocouple    readings 

Tavg=  ETn/n 
Tl=21.63    °C 

T2=21.41    °C 

T3=20.63  °C 
T4=20.94  °C 
T5=20.06  °C 
T6=20.40  °C 
Tavg=2  0.8  5    °C 
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standard  deviation=  (  (E  (Tn-Tavg)  2)  /n)  1/2 
standard  deviation=  0.55  °C 


The  logarithmic  term  in  the  equation  for  Two  can  be 
neglected  for  the  uncertainty  analysis  as  it  is  very  small  in 
comparison  to  the  standard  deviation. 
Two=Tavg=20.85  °C 
6Two=standard  deviation= . 55  °C 
AT=18.62  °C 

6AT/AT=(  (5Two/AT)2+(6Tsat/AT)2)  1/2 
5AT/AT= (  (.55/18. 62) 2+(. 01/18.  62 )2)1/2 
6AT/AT=2.9  5  % 

D.   UNCERTAINTY  IN  HEAT  FLUX 

q=  (Qh-2-Qf)/Ab 

Qh=VI     Qh=890.56  W         6Qh=12.2  W 

assuming  the  same  proportion  in  the  uncertainty  for  Qf : 

q=86000  W/m2      Ab=  0.010137  m2 

Qh-2 -Qf=  871.81  W 

Qf=  9 .38  W         6Qf=0.128  W 

6q/q= [ (6Qh/ (Qh-2Qf ) ) 2+ (26Qf/ (Qh-2Qf ) ) 2+ (6Ab/Ab) 2] 1/2 

6q/q= [ (12. 2/  (890.  56-2 -9.  38)  )2+  (2-  (.12 8)/  (890.  56-2-9. 38)  )2  + 

(6.39xl0"5/.010137)2]  1/2 
6q/q=  1.54  % 
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E.   UNCERTAINTY  IN  BOILING  HEAT -TRANSFER  COEFFICIENT 

h=q/AT 

6h/h  =  [ (6q/q)2+(6AT/AT)2]  1/2 
<5h/h=  [  (  .0154)2+( .0295)2] 1/2 
<5h/h=  3.3  % 

TABLE  6.   UNCERTAINTY  ANALYSIS  OF  FOUR  DATA  POINTS 


Parameters 

Smooth  tube 
4019  W/m2 

Smooth  tube 
86000  W/m2 

High  Flux 
2244  W/m2 

High  Flux 
88860  W/m2 

6Qh/Qh  (%) 

3.52 

1.37 

8.11 

1.36 

aAb/Ab  (%) 

.63 

.63 

.634 

.634 

6AT/AT  (%) 

5.20 

2.95 

17.4 

12.5 

6q/q   (%) 

1.29 

1.54 

8  .41 

1.5 

6h/h   (%) 

5.4 

3  .3 

19.3 

12  .6 

All  uncertainty  analysis  calculations  were  conducted 
during  decreasing  heat  flux  runs.  The  wall  superheat 
uncertainty  (6AT/AT)  for  the  HIGH  FLUX  tube  was  significantly 
greater  than  that  of  the  smooth  tube  due  to  a  greater 
deviation  in  the  thermocouple  temperature  readings. 
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APPENDIX  E.  SETUP  PROGRAM 

The  following  setup  program  was  used  in  the  preparation  of 
the  apparatus  prior  to  commencing  testing.  The  program  was 
written  in  Hewlett-Packard  Basic  5.0  for  both  the  Hewlett- 
Packard  9300  and  9852A  series  data  acquisition/control  unit. 
The  setup  program  performs  the  following: 

1.  Monitor  the  Sump  temperature 

2.  Monitor  the  evaporator  liquid  temperature 

3.  Measurement  and  readout  of  all  thermocouple  channels. 

4.  Measurement  and  readout  of  the  power  supplied  to  the 
tube  cartridge  heater. 

5.  Measurement  and  readout  of  the  power  supplied  to  the 
auxiliary  heaters  if  used. 
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1  i  PROGRAM:  SETUP 

2  i  DATE:   AU6UST  3.1991 

4  i  PROGRAMMER:  l_T  DEAN  SUGIYAMA 

5  i  MOOIFIED  Br  LANNIE  LAKE  JAN  22,  1992 

10  COM  /Cc/  C( 7 ) 

20  DATA  0.I00S509I  .2572?  .  94369  , -757345 . 9295 ,78025535 . 8  I  ,-9247486589 ,B . 97SSSE* 

1 1  .-2.5S192E*13 

21  OATA    3.94078EH4 

22  REAO   C(  ♦  ) 

23  ON  KEY  1  ,15  60T0  27 
25  PRINTER  IS  1 

27  PRINT 

28  PRINT 

30  PRINT  USING  "4X,  '"SELECT  OPTION"' 

31  PRINT  USIN6  '6X  ,"0=MONITOR  SUMP"" 

32  PRINT  USIN6  " 6X , ' " l =MONl TOR  LIQUID'" 

33  PRINT  USING  '6X,"2=CHECK  THERMOCOUPLES"" 

34  PRINT  USIN6  '6X."3=CHECK  MAIN  HEATER"" 

35  PRINT  USIN6  '6X,"4-CHECK  AUX  HEATERS'" 

36  PRINT  USING  '6X,"S  =  EXIT  PROGRAM'" 

37  PRINT  USIN6  '4X,"N0TE:  KEY  I  =  ESCAPE"" 

38  BEEP 

40  INPUT  Ido 

41  IF  Ido;5  THEN  Ido-5 

42  IF  ldo  =  0  THEN  50 

43  IF  Ido-I  THEN  155 

44  IF  Ido-2  THEN  173 

45  IF  Ido-3  THEN  195 

46  IF  Ido-4  THEN  195 

47  IF  Ido-5  THEN  23! 

48  PRINT 
49' 

50  PRINT 

51  PRINT  "SUMP  TEMPERATURE  DES  C  " 

53  PRINT 

54  OUTPUT    709i"AR    AF 1  1    ALII    L>R5" 
60  OUTPUT    709, 'AS    SA" 

70  Sum»0 

80  FOR    J-l    TO   5 

90  ENTER  709iE 

100  Sum=Sun+E 

!  10  NEXT  J 

120  Eave»Sum/5 

1  30  Ter»p»FNTvsv(Eave  ) 

140  PRINT    USING    "4X .MDD.DD"  iTenp 

141  BEEP 

142  PRINT 
150  UAIT  5 
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!5I 

GOTO  50 

152' 

!55 

PRINT 

156 

PRINT  "LIQUID  TEMFERATURE  OES  C 

158 

PRINT 

159 

OUTPUT  709; 'AR  AF08  AL09  ^RS" 

160 

Sun  =  0 

151 

FOR  1-1  TO  2 

is: 

OUTPUT  709  i "AS  SA" 

153 

ENTER  709;E 

164 

5un»Sun*E 

155 

NEXT  I 

166 

Eave»Suw/2 

167 

Temp-FNT*sv<  Eave  i 

168 

PRINT  USING  -4X  .MDO.DO'iTemp 

169 

BEEP 

170 

WAIT  5 

171 

60T0  155 

172' 

173 

PRINT 

176 

PRINT  "CHANNEL      TEMPERATURE 

177 

OUTPUT  709i'AR  AF00  ALII  VRS" 

178 

FOR  I-l  TO  12 

179 

OUTPUT  709;  "AS  SA" 

180 

Suf»  =  0 

181 

FOR  J-1  TO  5 

182 

ENTER  709iE 

183 

Su«"=Sun*E 

184 

NEXT  J 

!85 

Ea*e-Su«.'  5 

186 

Te*c"FNTvsv(Eave ) 

187 

PRINT  TABi 3 >;  I  :TAB< IS >;Te«p 

188 

NEXT  I 

189 

BEEP 

190 

WAIT  5 

191 

GOTO  173 

I94i 

195 

PRINT 

196 

OUTPUT  709; "AR  AF20  AL22  vRS- 

197 

FOR  1-1  TO  3 

198 

Output  709; "AS  sa- 

199 

5un  =  0 

200 

FOR  J»1  TO  5 

201 

ENTER  709;E 

202 

Su«»Sumf E 
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233  NEXT  J 

204  IP  1  =  1  THEN  VoIfSum/S 

205  IF  1-2  AND  Iao-3  THEN 

209  PRINT  -MAKE  SURE  VOLTAGE  BOX  IS  SET  TO  MAIN  HEATERS" 

210  Amp«Sum/5 
21  I  ENO  IF 

212  IF  1-3  AND  Ido-4  THEN 

216  PRINT  "MAKE  SURE  VOLTAGE  BOX  IS  SET  TO  Aux  HEATERS' 

2 1 7  Awp*Sufi/5 
2!9  END  IF 

219  NE<T  I 

220  Anp»A8S( Anp» I  . 91 82  ) 

221  Vol  t=ABS( Vol t «2S  ) 

222  Power-Vol t »Anp 

223  Resistence-Vol t /Amp 

224  PRINT 

225  BEEP 

226  PRINT  -VOLTAGE(V)  CURRENT(A)  RES ISTENCE< oh«s I  POWER(W)" 

227  PRINT 

228  PRINT  USING  - 1 X  ,5( MDOOO .DO ,4X )"  i  Vol  t  ,A«p  ,  Res  i  stence  .Power 

229  WAIT  5 

230  6OT0  195 

231  BEEP 

232  PRINT 

233  PRINT  'THAT'S  ALL  FOLKS'  ' 

234  ENO 

235  DEF    FNTvsvl V  ) 

236  COn    /Cc/    C( 7  ) 

237  T«C<0> 

238  FOR    1-1    TO    7 

239  T-T+C  I  >»V*I 

240  NEXT    I 

24!  T  =  T*6.S26397E-2»T.'  3.761  I  99E-3-T»S  .  0689Z59E-5  > 

242  RETURN  T 

250  FNEND 
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APPENDIX  F.   PROGRAM  DRPGB 

The  data  acquisition  and  reduction  program  DRPGB  is 
written  in  Hewlett-Packard  Basic  5.0  for  the  HP  9300  series 
computer  and  listed  on  the  following  pages. 
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10  I 
20  I 
30   ' 

41    ! 

50 

60 

70 

80 

90 

100 

I  10 

120 

130 

140 

150 

151 

160 

170 

180 

190 

Et(  19 

200 

210 

220 

230  > 

240 

250 

250 

270 

280 

290 

300  i 

310  i 

320  i 

330  I 

340 

350  i 

351 

360 

370 

380 

390 

400 

41  0 

420 

430 

431 

432 

433 

440 
450 
460 
470 
430 
490 
500 
510 

520 
530 
540  i 

550  ! 


FILE  NAME:  DRPG8 
DATE:       October  19,  1984 
REVISED:    July  11,  1991  BY  Doan  Sugiyema 

REVISED:    Aug   '8.  1992  BY  George  Bertsch  (including  R-124  properties) 
COM  /I  dp/  I  dp 
PRINTER  IS  I 
CALL  Select 

INPUT  "WANT  TO  SELECT  ANOTHER  OPTION  <l-Y,0-N>?"  "sel 
IF  Isel-I  THEN  GOTO  70 
BEEP 
BEEP 

PRINTER  IS  1 

PRINT  "DATA  COLLECTION/REPROCESSING  COMPLETED" 
END 

SUB  Main 
COM  /Iprop/  Ift 
COM  /lap/    Idp 
COM  /Cc/  C(7)  ,Ical 
COM  /Ulil/  D2  ,Di  ,Do  ,L  ,Lu  ,Kcu 

DIMEnf(12).T(12),D1a(13),D2a(l3),Dia(l3)  ,Doa( 13  >  ,La( 13  )  ,Lu*( 13  )  ,Kcua(  13) 
)  ,Tn$(4  U 151 

DATA  0.10086091  .25727.94369,-767345.8295,78025595.81 
DATA  -9247486589. 6. 97688E+1 1  ,-2 .66 1 92E+ 1 3  ,3 . 94078E+ 1 4 
READ  C( • > 

DATA  "Snoot h"  ,  "High  Flux "  , " Thermoexe 1-E  "  , " Thermoexel-HE" 
DATA  Smooth, High  F lux  ,Turbo-B  .High  Flux  Mod,TurDo-B  Mod 
READ  Tn$( •  ) 
PRINTER  IS  701 
BEEP 

IF  Idp-4  THEN  "PRINTER  IS  1 
IF  Idp-4  THEN  GOTO  2660 

INPUT  "ENTER  MONTH,  DATE  AND  TIME  ( MM : DD : HH : MM  :  SS  ) "  ,Dat eS 
OUTPUT  709i"TD"sDateS 
OUTPUT  709i"TD" 
ENTER  709iDate* 
PRINT 

PRINT  "  Month,  date  and  tine  :";Date$ 

PRINT  "  Date  : "  ,DATES( TIMEDATE  ) 

PRINT 

PRINT  USIN6  '10X.--NOTE:  Program  name  :  ORPGB""" 
BEEP 

INPUT  "ENTER  DISK  NUMBER" ,Dn 
PRINT  USIN6  *1GX  ."Disk  number   -  ■•,ZZ"tDn 
BEEP 

INPUT  "ENTER  INPUT  MODE  ( 0-3054A  ,  1 -FILE  ) "  ,  Im 
BEEP 
Ift-0 

INPUT  "SELECT  FLUID  < 0-1  1 4 , 1  -  1 24  ) "  ,  I f t 
BEEP 

INPUT  "SELECT  HEATIN6  MODE  (0-ELEC;  1 -WATER >".  Ihm 
BEEP 

INPUT  "ENTER  THERMOCOUPLE  TYPE  ( 0-NEW . 1 -OLD >".  leal 
IF  Im-0  THEN 
BEEP 

INPUT  "GIVE  A  NAME  FOR  THE  RAW  OATA  FILE" ,D2_f i leS 
PRINT  USIN6  "IBX, *"Naw  file  name:  " " , I 4A" ; D2_f l let 
Sizel-20 

CREATE  BDAT  D2_f l le*  ,Si ze I 
ASSI6N  »File2  TO  D2_fileS 

OUMMY  FILE  UNTIL  Nrun  KNOWN 
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560  OI_f  ileS-'DUMMV 

578  CREATE   BOAT   D1 J 1  let  .Sizel 

580  ASSIGN   "Filel    TO   Dl.fileS 

590  OUTPUT   CFilel  iDatet 

see  if  ihn-e  then 

sie  beep 

62e  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  (0-OEFAULT,  4  MAX.  >-. ldtc 

S3e  IF  Idtc-0  THEN 

6*e  Ldtci-e 

550  Ldtc2-e 

651  Ldtc3-0 

S52  Ldtc4-0 

66e  PRINT  USING  "ISX/'No  defective  TCs  exist*  "' 

S7B  END  IF 

S80  IF  Idte-I  THEN 

690  BEEP 

700  INPUT  'ENTER  DEFECTIVE  TC  LOCATION  ( ) -8  ) "  .Ldtc I 

710  PRINT  USIN6  *1BX,""TC  is  defective  at  location  •-.DD'iLdtcl 

720  Ldtc2-0 

721  Ldtc3-0 

722  Ldtc4-0 
730  END  IF 

740  IF  Idtc-2  THEN 

750  BEEP 

760  INPUT  "ENTER  DEFECTIVE  TC  LOCATIONS  < I -8  ) "  .Ldtc I  ,Ldtc2 

770  PRINT  USIN6  "16X,"TC  are  defective  at  locations  ' * ,DD ,4X  ,DO"  ;Ldtc  I  ,Ldtc2 

77)  Ldtc3-0 

772  Ldtc4-0 

780  ENO  IF 

781  IF  Idtc-3  THEN 

782  BEEP 

783  INPUT  'ENTER  DEFECTIVE  TC  LOCATION  <  1 -8 >' .Ldtc 1  ,Ldtc2  ,Ldtc3 

784  PRINT  USIN6  "IBX/'TC  is  defective  at  location  " "  ,3(  4X  ,D0  ) "  1  Ldtc  1  ,Ldtc2  ,Ld 
tc3 

785  Ldtc4-0 

786  END  IF 

787  IF  ldtc-4  THEN 

788  BEEP 

789  INPUT  'ENTER  OEFECTIVE  TC  LOCATION  (  1  -8 ) ■ .Ldtc 1  ,Ldtc2  ,Ldtc3  ,Ldtc4 

790  PRINT  USING  -16X.--TC  is  defective  at  location  ■ "  ,4(  4X ,OD )" iLdtcl  ,Ldtc2 .Ld 
tc3  ,Ldtc4 

792  END  IF 

793  IF  Idtc>4  THEN 
800  BEEP 

810  PRINTER  IS  I 

820  BEEP 

930  PRINT  'INVALID  ENTRY* 

840  PRINTER  IS  701 

850  SOTO  610 

860  END  IF 

870  END  IF 

880  OUTPUT  BFilel iLdtcl  ,Ldtc2  ,Ldtc3.Ldtc4 

890  1  l-i-l  ootion 

900  ELSE 

910  BEEP 

920  INPUT  "6IVE  THE  NAME  OF  THE  EXISTING  DATA  FILE'  ,D2_f 1  lei 

930  PRINT  USIN6  "lBX,**Old  file  name:  ' *  .  1 4A • 1 D2_f 1  lei 

940  ASSI6N  9File2  TO  D2_file« 

950  ENTER  •File2iNrun 

960  ENTER  •File2iOold» 

970  PRINT  USIN6  "16X,--This  data  sat  tanen  on  :  " *  ,  1 4A"  iDoldS 

980  ENTER  »File2iLdtcl  ,Ldtc2  .Ldtc3  ,^dtc4 

990  IF  Ldtc1>0  OR  Ldtc2>0  OR  Ldtc3>0  OR  Ldtc4>0  THEN 

1000  PRINT  USIN6  • 16X  ,' ' Thermocouples  yere  defective  at  locations: ** ,*<3D  ,4X  )* I 

Ldtcl  ,Ldtc2  ,Ldtc3  ,Ldtc4 

1010  END  IF 
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1020 

1030 

1040 

1050 

1060 

I0EI 

1062 

1070' 

1080 

1090 

100 

I  10 

120 

130 

140 

150 

160 

170 

180 

190 

1200 

1210 

1220 

1230 

1240 

1250 

260 

1270 

1280 

290 

300 

1310 

1320 

1330 

340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

420 

1430 

440 

1450 

460 

1470 

480 

1490 

500 

1510 

520 

530 

540 

1550 

1560 

1570 

1580 

590 

1600 

1610 

1620 

1630 

1640 

1650 


ENTER  SFnoZiltt 

END  IF 

Idtc-0 

IF  Ldtc1>0  THEN  Idtc-Idtc+1 

IF  Ldtc2>0  THEN  Idtc-Idtc+1 

IF  Ldtc3>0  THEN  Idtc-Idtc+1 

IF  Ldtc4>0  THEN  Idtc-Idtc+1 

IF  In-0  AND  Ihm-1  THEN  1595 

BEEP 

INPUT  "UANT  TO  CREATE  A  PLOT  FILE''  (  0-N  ,  1  -Y  ) "  ,  Ip  lot 

IF  Iplot-1  THEN 

BEEP 

INPUT  -GIVE  NAME  FOR  PLOT  FILE'  ,P_f i leS 

CREATE  BDAT  P_fileS,4 

ASSIGN  ©Plot  TO  P_fileS 

END  IF 

IF  Ihm-I  THEN 

BEEP 

INPUT  'WANT  TO  CREATE  Uo  FILE7  < 0-N  .  1 -Y  )  ■  ,  Iuf 

IF  Iuf- 1  THEN 

BEEP 

INPUT  "ENTER  Uo  FILE  NAME",UfileS 

CREATE  BDAT  UfileS.4 

ASSIGN  eufile  TO  UfileS 

END  IF 

BEEP 

INPUT  -WANT  TO  CREATE  Re  FILE7  < 0-N  ,  1 -Y  )"  ,  Ire 

IF  Ire-1  THEN 

BEEP 

INPUT  "ENTER  Re  FILE  NAME"  ,Ref i leS 

CREATE  BDAT  Ref  ileS,10 

ASSIGN  9Refile  TO  RefileS 

END  IF 

END  IF 

PRINTER  IS  1 

IF  Im-0  THEN 

BEEP 

PRINT  USING  '4X 

IF  Ihm-0  THEN 

PRINT  USIN6  "6X  ,-"0 

PRINT  USING  "6X  .""  1 

PRINT  USIN6  "6X  ,""2 

PRINT  USIN6  "6X  ,"3 

PRINT  USIN6  "6X  ,"4 

PRINT  USIN6  "6X  ,'"5 

PRINT  USIN6  "6X  ,""6 

PRINT  USIN6  '6X  ,"7 

PRINT  USING  "SX  .-'8 

PRINT  USIN6  'EX  ,-'9 

PRINT  USIN6  "6X 

PRINT  USING  'SX 

PRINT  USIN6 

PRINT  USIN6 

ELSE 

PRINT  USIN6 

PRINT  USIN6 

PRINT  USING 

PRINT  USIN6 

PRINT  USIN6 

END  IF 

INPUT  Itt 

OUTPUT  IFilal iltt 

END  IF 

PRINTER  IS  701 

IF  1 1 1 < 1 0  THEN  PRINT  USING  *  t SX 

IF  Itt>9  THEN  PRINT  USING  "1GX," 


Select  tube  number"" 

Smooth  4  inch  Ref""' 

Smooth  4  inch  Cu  (Press/Slide)' 

Soft  Solder  4  inch  Cu""" 

Soft  Solder  4  inch  HI6H  FLUX""1 

Uieland  Hard  8  inch""" 

HI6H  FLUX  8  inch"" " 

6EUA-K  40  Fins/ in""" 

6EUA-K  26  Fins/in""" 

6EUA-T  19  Fina/in"'" 

9  SEWA-T  OR  6EUA-TY  26  Fins/in""' 

10  THERh-OEXCEL-E'*" 

11  THERMOEXCEL-HE"  " 


6X 

•"12 

TURBO-B"" 

6X 

""13 

SEWA-K  19  Fin 

6X 

•■0 

Snoot h  tube  "  " 

6X 

" "  1 

High  Flux""  " 

EX 

"•2 

Turbo-B* " " 

6X 

•■3 

High  Flux  Hod 

6X 

••4 

Turbo-B  Hod" " 

""Tube  Number: 

'TuDe  Number : 


"  ,0'iltt 

.00"  :  Itt 
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1660  IF    inm-l    TutN   PRINT    usiNb      Isa,       moe    iype;  ,  i  an    i  mauui 

1670  BEEP 

1680  INPUT  'ENTER  OUTPUT  VERSION  ( 0-LON6 , 1 -SHORT  ,2-NONE  >"  .  Iov 

1690  BEEP 

1700  INPUT  "SELECT  ( 0-LIQ  ,  I -VAP  ,2-<  LIQ+VAP >/2  ) '  ,  Ilqv 

17101 

1720'  DIMENSIONS  FO  TESTED  TUBES 

1730'  ELECTRIC  HEATED  MODE 

17401  Dl"0ianeter  at  tnermocouo le  positions 

1750  DATA  .01 1 1 125,. 011 1 125,. 0M l 1 25  .  .01 29540  ,. 01 2446  .  .01 29540  ,. 0 1 00965 

1760  DATA  .0100965, .01 157, .01 157, .01 157, .01157, .01 157, .0100965 

1770  READ  Dla< • > 

1790  Dl-DlaUtt) 

1790' 

18001  D2«Diameter  of  test  section  to  the  base  of  fins 

1810  OAT A  .01 5875,. 01 5875,.  01 5875,. 01 5824  ,  .0 1 5875  .. 0 1 5824  ..01270 

1820  DATA  .0127, .0138, .0138, .0138, .0138. .0138, .0127 

1830  REAO  D2a( •  ) 

18401 

18501  Di»Inside  diameter  of  unenhanced  ends 

1860  OATA  . 0127, .0127, .0127, .0132, .0127, .0132  ,.011 1125. .011 1 125 

1870  DATA  .01 18  ,.01 18  ,.01 18, .01 18 ,.01 18 ,.01  I  I  125 

1880  READ  Dia( •  I 

1890' 

1900'  Do"0utside  diameter  of  unenhanced  ends 

1910  DATA  .0I5B75..015875..0I5B75,.015824,.015875..015824,.01270,.01270 

1920  DATA  .01331  ..01331  ,.01331  ,.01331  ..0158, .0127 

1930  READ  OoaC  I 

1940' 

1950'  L-Length  of  enhanced  surface 

1960  DATA  .1016, .1016, .1016. .1016  ..2032,. 2032  ,  .2032  .  .2032  .  .2032 ,.2032 ,.2032  ..20 
32. .2032. .2032 

1970  READ  La(» ) 

1980' 

1990'  Lu"Length  of  unenhanced  surface  at  the  ends 

2000  OATA  .0254  ..0254  ..0254  ,.0254,. 0762  ..0762  ..0762  ,.0762  ,.0762 ..0752 ,.0762 ..07 

52  ..0762  ,.0762 

2010  READ  Luai • > 

2020' 

2030'  Kcu-Thernal  Conductivity  of  tube 

2040  DATA  398.344,344,45.344.45.344,344,398.398.398,398.398,344 

2050  READ  Kcual ♦  ) 

2060  IF  Ihn-1  THEN 

2070i 

2080'  Data  statements  for  water  heating  node 

2090' 

2100  DATA  0.015875,0.015875,0.0169.0.0138.0.0169,0,0.0,0.0 

2110  READ  D2a< •  ) 

2120  DATA  0.0127,0.0127,0.0145.0.0127.0.0145.0.0.0,0.0 

2130  READ  Oia< •  I 

2140  DATA  0.015875.0.015875.0.0169.0.015875.0.0169.0,0.0.0.0 

2150  READ  Doa<«  > 

2160  DATA  0.3048.0.3048,0.3048.0.3048.0.3048.0,0.0,0,0 

2170  READ  La( • > 

2180  OATA  0.0254,0.0254,0.0254,0.0254,0.0254.0,0.0.0,0 

21 90  READ  Luaf •  ) 

2200  OATA  398,45,398,45,398,0.0.0,0.0 

2210  REAO  Kcua( •  ) 

2220  END  IF 

2230  D2-D2a(  Itt  ) 

2240  Di-Oia(Itt) 

2250  0o-0oa( Itt  > 

2260  L-La(Itt) 

2270  Lu-Lua(  Itt  ) 

2280  Kcu-Kcua(  Itt  ) 

2290  <n-.8 
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230d  Fr-.3 

2310  IF  Itt-0  THEN  Cf-1 .70E+9 

2320   IF  It t>B  THEN  Cf-3 . 7037E+ 1 0 

2330  A-PI»<Do*2-Ol*2 >/4 

2340  P-PI»Do 

2350  IF  I hm-1  THEN 

23E0  BEEP 

2370  INPUT  "TUBE  INITIATION  MODE.  ( 1 -HOT  WATER ,2-STEAM  ,3-COLD  WATER )"  ,1  tin 

2380  IF  Itin-I  THEN  PRINT  USIN6  '16X,"Tube  Initiate:  Hot  Water"" 

2390  IF  Itm-2  THEN  PRINT  USINS  "16X,*"Tube  Initiate:  Stean"" 

2400  IF  Itin-3  THEN  PRINT  USINS  "1BX,"Tube  Initiate:  Cold  Water"" 

2410  INPUT  "TEMP/VEL  MODE:  < 0-T-CONST  ,V-DEC i 1 -T-DEC  .V-CONST i  2-T-INC  .U-CONST  ) ' 

Itv 

2420  IF  Itv-0  THEN  PRINT  USIN6  " 1 6X  , " "Tenp/Uel  Mode:  Constant /Decrees ing"  "  " 

2430  IF  Itv-I  THEN  PRINT  USIN6  "  1  6X  ,  "Tenp/Uel  Mode:  Decreasing/Constant"" 

2440  IF  Itv-2  THEN  PRINT  USINS  " 1 6X  , " "Tenp/Vel  Mode:  Increasing/Constant""" 

24S0  INPUT  "WANT  TO  RUN  WILSON  PLOT?  (  1 -Y  ,0-N  ) "  ,  lui  1 

24S0  IF  Ihn-I  AND  Iwil-0  THEN 

2470  IF  Itt-0  THEN  Ci-.032 

2480  IF  Itt-1  OR  Itt-3  THEN  Ci-.059 

2490  IF  Itt-2  OR  Itt-4  THEN  Ci-.062 

2500  BEEP 

2510  INPUT  "ENTER  CI  ( DEF  :  WH- . 032  ,HF- .059  ,TB- . 062  ) '  ,Ci 

2520  PRINT  USING  " 1 6X  , " "Sieder-Tate  """ 

2530  PRINT  USING  *  1 SX  . " "   Constant    -  ",Z.4D'iCi 

2540  END  IF 

2550  END  IF 

2550  IF  Ihn-I  AND  In-1  AND  luil-l  THEN 

2570  IF  Itt-0  THEN  Ci-.032 

2580  IF  Itt-I  OR  Itt-3  THEN  Ci-.0S9 

2590  IF  Itt-2  OR  Itt-4  THEN  O-.062 

2S00  ASSIGN  9File2  TO  • 

2S10  CALL  WilsontCf  ,Ci  ) 

2S20  ASSIGN  9File2  TO  D2_fileS 

2530  ENTER  9F i le2 sNrun  .DoldS  ,Ldtc I  .Ldtc2  ,Ldtc3  ,Ldtc4  ,1 1 1 

2S40  END  IF 

2S50  Nsub-0 

2SS0  IF  Idp-4  THEN  Ihn-1 

2E70  IF  Ihn-1  THEN  NsuD-8 

2680  Ntc-6 

2690  IF  Ihn-0  THEN  Ntc-12 

2700  J-l 

2710  Sx-0 

2720  Sy-0 

2730  Sxs-0 

2740  Sxy-0 

2"'S0  Repeet:  i 

27S0  IF  In-0  THEN 

2780  Ido-2 

2',90  ON  KEY  1,15  RECOVER  2750 

2800  PRINTER  IS  I 

2810  PRINT  USING  "4X, "SELECT  OPTION'" 

2820  PRINT  USIN6  'SX,"0-TAKE  DATA""' 

2830  IF  Ihn-0  THEN  PRINT  USIN6  "6X,""1-SET  HEAT  FLUX"" 

2840  IF  Ihn-I  THEN  PRINT  USING  "6X.""l-SET  WATER  FLOW  RATE'" 

2850  PRINT  USIN6  "6X,"2-SET  Tsaf" 

2360  PRINT  USIN6  "4X,"N0TE:  KEY  I  -  ESCAPE'"" 

2870  BEEP 

2880  INPUT  Ido 

2890  IF  Ido>2  THEN  Ido-2 

2900  IF  Ido-0  THEN  4440 

29101 

29201  LOOP  TO  SET  HEAT  FLUX  OR  FLOWMETER  SETTING 

2930  IF  Ido-I  THEN 

2940  IF  Ihn-0  THEN 
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2950  OUTPUT  709i-AR  AF20  AL2l  VRS" 

2960  BEEP 

2970  INPUT  'ENTER  DESIRED  Qdp-  .Dodp 

2980  PRINT  U5IN6  '4X  ,'  'DESIRED  Oao   ACTUAL  Qdp'" 

2990  Err- I  000 

3000  FOR  I-l  TO  2 

3010  OUTPUT  709i'AS  SA* 

3020  Sun-0 

3030  FOR  J i-l  TO  S 

3040  ENTER  709 i E 

3050  Sun-Sunt-E 

3060  NEXT  Ji 

3070  IF  I-l  THEN  Volt-Sun/5 

3080  IF  1-2  THEN  Anp-Sun/5 

3090  NEXT  I 

3091  PRINT  'Vs:  '  .Volt 

3092  PRINT  'Is:'  .Ano 
3100  Anp-ABS(  ArnD»1  .9182  ) 
3110  Volt-ABS(Volt»25) 

3120  AQdp-Volt»Amp/(PI»D2»L > 

3130  IF    ABS< Aadp-Dadp >>Err   THEN 

3140  IF  AQap>Dqdp  THEN 

3150  BEEP  4000. .2 

3160  BEE3  4000, .2 

3170  BEEP  4000, .2 

3180  ELSE 

3190  BEEP  250.. 2 

3200  BEEP  250,. 2 

3210  BEEP  250.. 2 

3220  END  IF 

3230  PRINT  USIN6  '  4X .HZ . 3DE ,2X ,MZ . 3DE' lOqdp  .Agdp 

3240  UAIT  2 

3250  GOTO  3000 

3260  ELSE 

3270  BEEP 

3280  P»:nT  USIN6  '4X .MZ.30E.2X  .MZ.3DE" iDodp  ,AQdp 

3290  Err-500 

3300  WAIT  2 

3310  60T0  3000 

3320  END  IF 

3330  ELSE 

3340  BEEP 

3350  INPUT  "ENTER  FLOWMETER  SETTING* .Fma 

3360  60T0  2810 

3370  END  IF 

3380  END  IF 

3390' 

3400i  LOOP  TO  SET  Tsat 

3410  IF  Ido-2  THEN 

3420  IF  Ikdt-1  THEN  3470 

3430  BEEP 

3440  INPUT    "ENTER   DESIREO   Tsaf.Dtld 

3450i  PRINT  USIN6  '4X,"  DTaat   ATsat    Rate     Tv     Rate'" 

3460  Ikdt-1 

3470  Oldl-0 

3480  Old2-0 

3490  Nn-1 

3500  Nra-Nn  MOD  15 

3510  Nn-Nn+1 

3520  IF  Nrs-I  THEN 

3530  IF  Ihn-0  THEN  PRINT  USIN6  UX,"   Tsat     Tldl     Tld2     Tv     Taump'" 

3540  IF  I  hit- I  THEN  PRINT  USIN6  *4X,"   Tsat    Tldl     Tld2     Tv    Tsunp   Tinle 

t    Tp i le   Tout " " " 

3550  ENO  IF 

3560  IF  Ihm-0  THEN  OUTPUT  709i'AR  AF08  ALII  VR5' 

35"*0  IF  Ihn-I  THEN  OUTPUT  709;  "AR  AF0  AL5  VR5" 
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3530  FOR  i-l  TO  6 

3590  IF  Ihn-0  AND  I>4  THEN  3860 

3500  Suki-0 

3510  OUTPUT  709i  "AS  Sft" 

3620  FOR  Ji-1  TO  20 

3630  ENTER  709iEliq 

36*0  Sun-Sun-t-El  lq 

3650  NEXT  Jj 

3660  EliQ-Sun/20 

3670  Tld-FNTvsv(EliQ  ) 

3680  IF  I-l  THEN  Tldl-Tld 

3690  IF  1-2  THEN  Tld2-Tld 

3700  IF  1-3  THEN  Tv-Tld 

3710  IF  1-4  THEN  Tsunp-Tld 

3720  IF  1-5  THEN  Tmlet-Tld 

3730  IF  1-6  THEN  Tout-Tld 

3740  NEXT  I 

3750  IF  Ihn-I  THEN 

3760  OUTPUT  709) "AR  AF00  AL00  VR5" 

3770  OUTPUT  709i "AS  SAU 

3780  Sun-0 

3790  FOR  Kk-i  TO  20 

3800  ENTER  709iE 

3810  Sun-Sun+E 

3820  NEXT  Kk 

3830  Enf(7 >-ABS(Sun/20  ) 

3840  Tpile-Enf(7)/3.96E-4 

3850  END  IF 

3860  Atld-< Tldl+Tld2 >».5 

3870  IF  ABS( Atid-Dtld)>.2  THEN 

3880  IF    AtloVDUd   THEN 

3890  BEEP  4000,. 2 

3900  BEEP  4000,. 2 

3910  BEEP  4000, .2 

3920  ELSE 

3930  BEEP  250,.  2 

3940  BEEP  250,. 2 

3950  BEEP  250,. 2 

3960  END  IF 

3970  Errl-Atld-Oldl 

3980  Oldl-Atld 

3990  Err2-Tv-01d2 

4000  01d2-Tv 

4010  IF  Tldl>100.  THEN  4060 

4020  IF  Ihn-0  THEN  PRINT  USIN6  "  4X ,5(  NDDD . DD ,2X  ) " t Dt Id  ,Tldl  ,Tld2 ,Tv .Tsunp 

4030  IF  Ihn-1  AND  Idp-0  THEN  PRINT  USIN6  *4X ,7(  MOD . DO  ,2X  ) ' iDt Id  ,T  Id)  ,T ld2 ,Tv ,Ts 

unp  ,T  inlet  ,Tpi  le 

4040  IF  Ihn-1  AND  Idp-4  THEN  PRINT  USIN6  ' 4X  ,5( MDO  .  DD  ,2X  )  ,3( M3D .OD ,2X  ) ' i Dt Id  ,T1 

d  1  ,T1 d2  ,Tv  , Tsunp  ,T inlet ,Tpile .Tout 

4050  WAIT  2 

4060  60T0  3500 

4070  ELSE 

4080  IF  ABS(Atld-Dtld)>. 1  THEN 

4090  IF  Atld>Dtld  THEN 

4  100  BEEP  3000,. 2 

4110  BEEP  3000,. 2 

4120  ELSE 

4130  BEEP  800. .2 

4  140  BEEP  800,.  2 

4150  END  IF 

4160  Errl-Atld-Oldl 

4170  Oldl-Atld 

4180  Err2-Tv-01d2 

4190  01d2-Tv 

4200  IF  Ihn-0  THEN  PRINT  USIN6  " 4X ,5( MOOD . DD  ,2X >" iDt Id  .Tldl  ,Tld2  .Tv  , Tsunp 
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4210  IF  Ihti-l  THEN  PRINT  USING  •  4X  ,5(  MOD  .  DO  ,2X  >  ,3<  M3D  .  DD  , I  x  ) "  ;  Dt  i  d  ,  *  I  d  1  ,  T  1  d2  ,T  v 

,Tsunp, Timet  ,  <  pi  j.e  ,  i  out 

4220  WAIT  2 

4230  GOTO  3500 

4240  ELSE 

4250  BEEP 

42S0  Erri-Atld-Oldl 

4270  Oldl-Atld 

4280  Err2-Tv-01d2 

4290  01d2-Tv 

4300  IP  Ihw-0  THEN  PRINT  USING  ' 4X ,5( MODD .DO  ,2X  )  • i Dt Id  ,Tldl  ,Tld2  ,Tv  .Tsuno 

4310  IF  Ihn-I  THEN  PRINT  USING  ' 4X  ,8( MOO . DO  ,2X  )* ; Dt Id  .Tldl  ,Tld2 ,Tv .Tsunp  .T i niet 

,Tpi le  .Tout 

4320  UAIT  2 

4330  GOTO  3500 

4340  END  IF 

4350  END  IF 

43G0  END  IF 

4370'  ERROR  TRAP  FOR  I  do  OUT  OF  BOUNDS 

4380  IF  Ido;2  THEN 

4390  BEEP 

4400  GOTO  2810 

4410  END  IF 
4420' 

4430'  TAKE  DATA  IF  Im-0  LOOP 

4440  IF  Iko 1-1  THEN  4480 

4450  BEEP 

4460  INPUT  "ENTER  BULK  OIL  X'  .Boo 

4470  Ik.ol-1 

4480  IF  Ihn-0  THEN  OUTPUT  709i'AR  AF00  ALII  URS ' 

4490  IF  Ihn-1  THEN  OUTPUT  709i'AR  AF0  AL5  VRS" 

4500  IF  Ihn-0  THEN  Ntc-12 

4510  FOR  I-l  TO  Ntc 

4520  OUTPUT  709s "AS  SA" 

4530  Sun-0 

4540  FOR  Ji-1  TO  20 

4550  ENTER  709;E 

4560  Sum-Sum+E 

4570  IF  I-(S-NsuD)  OR  I-(10-Nsub)  THEN  Et(Ji-l)-E 

4580  NEXT  Ji 

4590  Kdl-0 

4600  IF    I-<9-Nsub>    OR    I-<10-N3ub>    THEN 

4610  Eave-5un/20 

4620  Sun-0. 

4630  FOR  Jk-0  TO  19 

4640  IF  ABS(Et(  Jk  >-EaveKS.0E-6  THEN 

4650  Sun-Sufi+Et<  Jk  > 

4560  ELSE 

4670  Kdl-Kdl+I 

4680  ENO  IF 

4690  NEXT  Jk 

4700  IF  I-(9-N3ub>  OR  I-(10-Naud>  THEN  PRINT  USING  ^X/'Kdl 

4710  IF  Kdl>l0  THEN 

4720  BEEP 

4-30  BEEP 

4^40  PRINT  USIN6  *4X,"Too  nuch  scattering  in  data  -  repeat  data  set 

4750  SOTO  2800 

4750  ENO  IF 

4770  END  IF 

4780  Emf  (  I  >-Supi/<  20-Kdl  ) 

4790  NEXT  I 

4800  IF  Ihn-1  THEN 

4810  OUTPUT  709i*AR  AF00  AL00  UR5 " 

4820  OUTPUT  709i 'AS  SA- 

4830  Sun-0 

4840  FOR  Kk-1  TO  20 


• ,DD' iKdl 
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4353  5N"^    709;E 

4553  3un»Sun*E 

4878  NEXT   Kk 

4830  E«f<7>-ABS(Sun>/20 

4830  END    IF 

4900  IF  Ihn-0  THEN 

4310  OUTPUT  709  )  "  AR  AF20  AL2  1  KIRS' 

4920  FOR  1-1  TO  2 

4930  OUTPUT  709i  "AS  SA" 

4940  Sum-0 

4950  FOR  Ji-l  TO  5 

49S0  ENTER  709iE 

4970  Sum-Sum+E 

4980  NEXT  Ji 

4990  IF  1-1  THEN  Vr- Sum/5 

5000  IF  1-2  THEN  Ir-Sum/5 

5010  NEXT  I 

5020  END  IF 

5030  ELSE 

5040  IF  Ihm-0  THEN  ENTER  9F i le2 i Boo  ,To IdS  ,Emf ( •  )  ,Ur  ,  Ir 

5050  IF  Ihm-1  THEN  ENTER  9Fi le2 sBop  ,Told$  .Emf ( •  )  ,Fms 

5060  ENO  IF 

5070i 

5080'  CONVERT  e«f'S  TO  TEMP  ,UOLT  .CURRENT 

5030  Twa-0 

5100  FOR  1-1  TO  Ntc 

51 10  IF  latc>0  THEN 

5120  IF  I-Ldtcl  OR  I-Ldtc2  OR  I-ldtc3  OR  I-Ldtc4  THEN 

5130  T(I>  —  93.39 

5140  60T0  5240 

5150  END  IF 

5160  END  IF 

5170  IF  Itt<4  AND  Ihm-0  THEN 

5180  IF  I>4  AND  I<9  THEN 

5190  T(I>  —  99.99 

5200  GOTO  5240 

5210  END  IF 

5220  ENO  IF 

5230  T( I  >-FNTvsv(Emf <  I  )  ) 

5240  NEXT  I 

5250  IF  Itt<4  THEN 

5251'    UARNIN6:  4  BAD  THERMOCOUPLES  ARE  NOT  POSSIBLE  FOR  ITT<4 

5260  FOR  I-l  TO  4 

5270  IF  I-Ldtcl  OR  I-Ldtc2  THEN 

5280  Tua-Tua 

5290  ELSE 

5300  Twa-Tua+T(  I  ) 

5310  END  IF 

5320  NEXT  I 

5330  Tua-Twa/(4-Idtc > 

5340  ELSE 

5350  IF  Ihn-I  THEN  5450 

5360  FOR  I-l  TO  8 

5370  IF  I-Ldtc!  OR  I-Ldtc2  OR  I-Ldtc3  OR  I-Ldtc4  THEN 

5380  Twa-Twa 

5390  ELSE 

5400  Tua-Twa*T( I  ) 

5410  END  IF 

5420  NEXT  I 

5430  Tu-Tua/( 8-Idtc > 

5440  END  IF 

5450  Tld-T( 9-Naub  ) 

5460  Tld2-T( 10-Nsub  ) 

5470  Tlda-<Tld*Tld2 >«.S 

5480  Tv-T( 1 1-Nsub  > 

5490  IF  Itt<3  AND  Ihm-0  THEN 
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5500 

55)0 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5500 

55)0 

5620 

5630 

5640 

5650 

5660 

5670' 

5680' 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820' 

5830' 

5340 

5850 

5860 

5870 

5880 

5890 

5900 

5910' 

5920' 

5930 

E-IB) 

5940' 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

5020 

6030 

5040 

6050 

6060 

6070 

6080 

6090 

6100 

61  10 

6120 

6130 

6140 


Tld2— 99.99 

Tv-i T( 10  )  +  T( 1  I  )  )/2 

END  IF 

Tsunp-T( 12-Nsub  > 

IF  Ihn-0  THEN  5570 

Tinlet-T< 13-NauD ) 

Tout-T( 14-Neub ) 

IF  Ihn-0  THEN 

Anp-ABS(  Ir* I  .9182  ) 

Yolt-ABS< Vr  )»25 

Q«Volt»A«p 

ENO  IF 

IF  Itt-0  AND  Ihm-0  THEN 

Kcu-FNKcu(Tu ) 

ELSE 

Kcu-Kcua< Itt  ) 

END  IF 

FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 

IF  Ihn-0  THEN  Tu-Tw-Q.|_OG(D2/Dl  >/ ( 2»PI  *Kcu»L  ) 

IF  Ilqv-0  THEN  Taat-Tlda 

IF  Ilgv-1  THEN  Tsat-'Tlda+Tv  )«.5 

IF  I1qv-2  THEN  Tsat-Tv 

IF  Ihn-1  THEN 

Tavg-T inlet 

6rad-37.9853+. 104388'Tavg 

Tdrop-ABS(Emf (7  )  )•!  . E*6/(  10»6rad> 

Tavgc-Tinlet-Tdrop».5 

IF  ABS(Tavg-Tavgc >> .01  THEN 

Tavg-( Tavg  +  Tavgc  >•  .5 

60T0  S7S0 

END  IF 

COMPUTE  WATER  PROPERTIES 
IF  Ihm-1  THEN 
Kw-FNKul  Tavg) 
rtuwa-FNMuw*  Tavg  ) 
Cow-FNCpwi Tavg  ) 
Pry-FNPrwf  Tavg  ) 
Rhow-FNRhow( Tavg  ) 
Tu i-Tavg 


.23688E-7-Fns»4. 31897 
Mdot»( 1 .0365-T inlet •<  1  . 96644E-3-T inlet »S  .  252E-6 >>/1 .0037 


Compute  MOOT 

Mdot-3.9657E-3+Fna»(3.61955E-3-Fns«(8.82006E-6-FnaM 
)  ) 

"dot 
Kdt-0 

Q-rtdot •Cpw»Tdrop 

Ln td-T drop /LOG ( (Tinlet-Tsat)/( T inlet -Tdrop-T sat >  ) 
Uo-Q/(PI*Oo*L*Lmtd ) 
Ru-0o»L06(Do/Di  )/( 2.»Kcu  ) 
Tu-T9at*Fr»Lmtd 
Uw-Hdot/(Rhow*P I«Dl*2/4  ) 
Rew-Rhow«Vu»Di/Muwa 

Hi-Cl»Ky/Dl«Rew" . 8»Prw" l  1/3.  >•<  Muwa/FNHuw( Tyi  )  >" .  14 
Tmic-Tavg-Q/(PI»Do»L»Hi  ) 
IF  ABS(Tui-Tuic  >>.0I  THEN 
Tui-( Tui*Tuic  )•  .5 
60T0  6030 
END  IF 

Tui-( Tu i +Tuic  )• .5 
Ho- I /(  I /Uo-0o/( Di*Hi  )-Ru  ) 
END  IF 
END  IF 

IF  Ih*-I  THEN 
Thetao-Q/iHo»PI»Do«L ) 
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6150 

s  1 6a 

6170 

51  80 

5:90 

6200 

6210 

6220 

6230 

62*0 

62501 

6260> 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380> 

6381 

6382 

6383 

6384 

6385 

6386 

6387 

6388 

6390' 

54001 
6410 

6420 
6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

651  l 

6520' 

6530' 

6540 

6550 

6560 

65701 

6580' 

6590 

6600 

6610 

66201 

6630' 

5640 

5650 

6660 

6670 

6680' 

6690' 

6700 

6710 


Tw-Tsat*Tnetao 

ELSE 

Thetab-Tn-Tsat 

END  IF 

IF  Thetab<0  THEN 

BEEP 

INPUT  -TWALL<TSAT  (0-CONTINUE 

IF  Iev-0  THEN  SOTO  2760 

IF  Iev-l  THEN  7220 

END  IF 

COMPUTE  VARIOUS  PROPERTIES 
Tf i ln-( Tw+Tsat  )•  .5 
Rho-FNRho( Tf t In  ) 
Mu»FNMu<  Tf  i  In) 
K-FNK( Tf iln > 
Cp-FNCp< Tf iln  ) 
Beta-FN8eta( Tf iln  ) 
Hfg-FNHfg( Tsat  ) 
Ni-Hu/Rho 
Alpha-K/(Rho»Cp  > 
Pr-Ni/Alpha 
Psat-FNPsaUTsat  ) 


I -END)'  .lev 


PRINT  "Rho: 
PRINT  "Mu:" 


,Rho 
Mu 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


,K 


"K  : ' 

'Cp: " ,Cp 

"Tf  iln: "  ,Tf  iln 

*Hfg: "  ,Hfg 

■Psat : ■  ,Psat 
PRINT  "Tsat :  "  ,Tsat 

COMPUTE  NATURAL-CONUECTIUE  HEAT-TRANSFER  COEFFICIENT 
FOR  UNENHANCED  END(S) 
HOar-190 

Fe«( Hoar«P/(Kcu»A ) )" .5«Lu 
Tanh-FNTanh( Fe  ) 
Theta»Tnetab»Tanh/Fa 

Xx«( 9.81 •Beta»Thetao»Do"3»Tann/( Fe«Ni»Alp ha))". 166667 
Yy(  l+(  .559/Pr  )*(  9/16  ))"!  8/27  ) 
Hbarc-K/Do»l  . 6+  .387«Xx/Yy  )"2 
IF  ABS( <Hbar-Hbarc  )/Hbarc )>.00l  THEN 
Hbar-( Hbar+Hbarc  )• .5 
SOTO  6420 
END  IF 
PRINT  "NAT  CONU  h:",Hbar 

COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  END(S) 

01"<Hbar«P»Kcu»A)".5»Thetab»Tanh 

Qc-0-2«Ql 

As-PI«D2»L 

COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

Qdp-Qc/As 

Htube-Qdp/Thetab 

Caf"(Cp«Thetab/Hfg)/( Qdp/(Hu»Hfg  >•<  .014/(9.81 »Rho  )" . 5  )* ( 1 /3 .  )*Pr*  1  .  7  ) 

RECORD  TIME  OF  DATA  TAKIN6 
IF  In-0  THEN 
OUTPUT  709TTD" 
ENTER  709iTold$ 
END  IF 

OUTPUT  DATA  TO  PRINTER 
PRINTER  IS  701 
IF  Iov-0  THEN 
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5728  PRINT 

S7301  PRINT  USIN6  "10X."'Oata  Set  Numoer  -  " "  ,DDD  ,2X  ," "8ulk  Oil  X  •  "  _  ,OD  .  D  ,SX  .  I 

4A" i J  ,Bop  .ToldS 

673!  PRINT  USING  *10X."D«ta  Set  Number  •  ■  •  ,DDD  .2X  , "  'Bulk  Oil  X  -  "".DD.D'iJ.B 

op 

6732  PRINT  ■  TIME : '  , TIMES* TIMEDATE  ) 

6740  IF  IhKi-8  THEN 

6750  PRINT  USIN6  "I0X,""TC  No :     I       2      3      4      5      6      7 

6760  PRINT  USIN6  "10X,""Tenp  :  ' '  ,8(  1  X  ,MDD  .  OD  )  ■  i  T(  I  )  ,T(  2  >  ,T(  3  >  ,T<  4  )  ,T<  S  )  ,T(  6  )  ,T( 

7)  ,T(8  ) 

6770  PRINT  USIN6  "10X,**  Tya     Tliqd    Tliad2   Tvapr   Psat    Tsunp""" 

6780  PRINT  USIN6  ' I 0X ,2< MOO .00  , 1  X  >  .  I  X  .MDD . DD  ,  I  X  ,2(  I  X .HDD .OD ) ,2X ,MDO .D" ; Tw  ,T  Id  ,T 

ld2 ,Tv  ,Psat  ,Tsunp 

6790  PRINT  USING  *10X,"  Thetab   Htuoe      Qdp""" 

5800  PRINT  USING  " I 0X  ,MDD  .  3D  ,  1 X  ,MZ . 3DE  .  1 X  ,MZ . 3DE" ; Thetab .Htube .Qdp 

6810  ELSE 

5820  PRINT  USING  " I 0X . ' "  Pns     Uy    Tsat   Txnl  Tdrop  Thetab     q         Uo 

Ho""" 

6830  PRINT  USING  " 1 0X  ,*< 20. 00, IX  )  .Z.3D.1X  .DD.DD  .IX  ,3<MZ.30E  ,1X  )" ;Fms  ,Uw  .Tsat  ,Ti 
nlet .Tdroo .Thetab ,Qdp  ,Uo  ,Ho 

6840  END  IF 

6850  END  IF 

6850  IF  Iov-I  THEN 

6870  IF  J-l  THEN 

5880  PRINT 

6890  IF  Ihm-0  THEN 

6900  PRINT  USING  "  1 0X , " *  RUN  No   OilX   Tsat     Htube      Qdp       Thetab'"" 

6910  ELSE 

6920  PRINT  USING  " 1 0X  , " "   FMS    OILX   TSAT    HTUBE      QDP       THETAB""" 

5930  END  IF 

6940  END  IF 

6950  IF  Ihw-0  THEN 

6960  PRINT  USING  " 1 2X  ,3D  ,4X  ,D0  ,2X  ,MDD . 00  ,3( 1 X  ,MZ . 30E  ) " t J  .Bop .Tsat .Htude .Qdp .The 
tab 

5970  ELSE 

6980  PRINT  USING  " I 2X  .3D  ,4X  .DO  ,2X  ,M00 . 00  ,3(  I  X ,MZ . 3DE  ) " ;Fms ,Bop .Tsat .Htuoe ,Qoo ,T 
netao 

6990  END  IF 

7000  END  IF 

7010  IF  In-0  THEN 

7020  BEEP 

7030  INPUT  "OK  TO  STORE  THIS  DATA  SET  ( 1-Y ,0-N )?"  ,0k 

7040  ENO  IF 

7050  IF  Ok-I  OR  I««1  THEN  J-J+l 

7060  IF  Ok-I  AND  In-0  THEN 

"'070  IF  Ihw-0  THEN  OUTPUT  BF  i  le  I  iBoo  .ToldS  ,Enf  (  •  )  ,Ur  .  lr 

7080  IF  Ihn-I  THEN  OUTPUT  9Fi le I i Boo  .ToldS  ,Enf ( •  )  ,Fms 

7090  END  IF 

7)00  IF  Iuf-I  THEN  OUTPUT  •UfiletUw.Uo 

7110  IF  Ire-I  THEN  OUTPUT  9Ref i le iFns  .Rew 

7120  IF  (Iw-1  OR  Ok-I)  ANO  Iplot-l  THEN  OUTPUT  9P lot iQdp  .Thetab 

7130  IF  In-0  THEN 

7140  BEEP 

7150  INPUT  "WILL  THERE  BE  ANOTHER  RUN  i  I -Y  ,0-N  P"  ,6o_on 

7160  Nrun-J 

7170  IF  6o_on-0  THEN  7220 

7180  IF  6o_on<>0  THEN  Repeat 

7190  ELSE 

^200  IF  J<Nrun*l  THEN  Repeat 

7210  END  IF 

7220  IF  In-0  THEN 

7230  BEEP 

7240  PRINT  USIN6  "10X.""NOTE:  "".ZZ.""  data  runs  were  stored  in  file  "".l0A"iJ- 
I  ,02_filet 

7250  ASSI6N  CFilel  TO  • 
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7263 
7270 
7280 
7290 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
plot 
7450 
7460 
7470 
7480 
7490 
7500 
7510 
7520 
7530 
75401 
7550' 
7560 
75701 
7580 
7590 
7600 
7610 
7620 
7621 
7622 
76301 
7640 
7650 
7651 
7652 
7655' 
7557 
7658 
7659 
7650 
7670 
7680 
7681 
7682 
7690' 
7700 
7710 
7720 
7721 
7722 
7723' 
7725 
T"5 
7726 
7727 
7730 


X-Y  pairs  u/ere  stored  in 


OUTPUT  8Fi ie2iNrun-l 

ASSIGN  BFiial  TO  Dl_f i leS 

ENTER  •Filet ;0atel.Ldtc1  ,Ldtc2.Itt 

OUTPUT  »File2iDateS  .Ldtcl  ,Ldtc2,Itt 

FOR  1-1  TO  Nrun-I 

IF  Ihn-0  THEN 

ENTER  ©File  I  iBop  .ToldS  ,Enf  (  •  )  ,\Jr  ,  Ir 

OUTPUT  SFile2iBop  ,  ToldS  ,Emt < • ) ,Vr ,Ir 

ELSE 

ENTER  OF  i  lei  ;8op  ,Toid$  ,Enf  (  •  )  .Fits 

OUTPUT  »File2iBoo  , ToldS  ,Emf ( • ) ,Fms 

ENO  IF 

NEXT  I 

ASSI6N  9Filel  TO  • 

PUR6E  "DUMMY" 

ENO  IF 

BEEP 

PRINT 

IF  Iplot-1  THEN  PRINT  USING  "IBX/'NOTE:  "".Z2 

data  file  " "  .  1  0A"  i  J-  1  ,P_f  i  leS 

ASSIGN  9File2  TO  • 

ASSIGN  ©Plot  TO  • 

IF  Iuf-I  THEN  ASSIGN  SUfile  TO  • 

IF  Ire-1  THEN  ASSIGN  SRefile  TO  • 

CALL  Stats 

BEEP 

INPUT  "LIKE  TO  PLOT  DATA  ( 1 -Y  ,0-N  )?"  ,Ok 

IF  Ok- 1  THEN  CALL  Plot 

SUBEND 


CURVE  FITS  OF  PROPERTY  FUNCTIONS 

DEF  FNKcu<  T  ) 

OFHC  COPPER  250  TO  300  K 

Tk-T+273. 15      'C  TO  K 

K-434-. 1 12»Tk  . 

RETURN  K 

FNENO 

DEF  FNMu(T  ) 

COM  /Iprop/  I  ft 

IF  Ift-0  THEN 

170  TO  360  K  CURUE  FIT  OF  YISCOUSITY 

Tk-T+273. 15      iC  TO  K 

Mu-EXP( -4.4636  +  <  101 1 .47/Tk  >  >•! .0E-3 

END  IF 

IF  Ift-I  THEN 

223  TO  373  K 

Mu-3.4S93-4.2GE-2«T+3.9485E-4»T-2-4. 1 93E-6»T"3+2 .0709E-8»T"4 

Mu-Mu*l .0E-4 

END  IF 

RETURN  Mu 

FNEND 

OEF  FNCp(T  ) 

COM  /Iprop/  Ift 

IF  Ift-0  THEN 

180  TO  400  K  CURUE  FIT  OF  Cp 

Tk-T+273. 15      'C  TO  K 

Cp-.40I88+! .65007E-3»Tk+l . 5 1 494E-5+Tk "2-6 . 67853E- 1 0»Tk "3 

Cp-Cp»1000 

END  IF 

IF  Ift-1  THEN 

223  TO  373  K  CURVE  FIT  Cp(R124) 

Cp-1  .0542  +  2. I405E-3  +  T+1  .0709E-5«T'2-6 . 472 1 E-8»T*3-1 . 4324E-9»T"4+4 . 1 36E- I  1* 


Cp-Cp»l000 
ENO  IF 
RETURN  Cp 
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7~40  FNEND 

7750  OEF  FNRhol T  ) 

7751  COn  /Iprop/  [ft 

7752  IF  Ift-0  THEN 

7750  Tk-T+273. 15  ! C    TO   K 

7770  X-1-< 1 .8«Tk/7S3.95>       IK    TO   R 

7780  Ro-36. 32*51 . US41 4»X*( I /3>+IG.*180l5«X*17.47B838»X*.5+l . I 138Z8»X*2 

7790  Ro-Ro/. 062428 

7791  END  IF 

7792  IF  Ift-I  THEN 

7793'  223  TO  373  K  (R-124  > 

7796  Ro-l434.8-2.8619»T-5.7267E-3»T-2-7.2852E-5»T-3 

7797  ENO  IF 
7800  RETURN  Ro 
7810  FNEND 

7820  DEF  FNPr(T) 

7830  Pr-FNCp< T  )»FNMu(T  )/FNK(T  ) 

7840  RETURN  Pr 

7850  FNENO 

7860  DEF  FNK(T) 

7861  COM  /Iprop/  Ift 

7862  IF  Ift-0  THEN 
78701  T<360  K  WITH  T  IN  C 

7880  K-.071-.000261»T 

7881  END  IF 

7882  IF  Ift-1  THEN 

7883'  223  TO  373  K   (R-124) 

7885  K-7.5191E-2-3.5436E-4.T-I . 954SE-7.T'2+3 . 1 835E-9.T"3 

7886  END  IF 
7890  RETURN  K 
7900  FNEND 

7910  DEF  FNTanh(X  ) 

7920  P-EXP(X) 

7930  O-l/P 

7940  Tanh-( P-Q >/(P+Q  ) 

7950  RETURN  Tanh 

7960  FNEND 

7970  DEF  FNTvsv(U) 

7980  COM  /Cc/  C(7)  .leal 

7990  T-C(0) 

8000  FOR  I-l  TO  7 

8010  T-T+C(I>»V*I 

8020  NEXT  I 

8030  IF  Ical-1  THEN 

8040  T-T-6.7422934E-2+T.(9.0277043E-3-T.(-9.3259917E-5)  ) 

8050  ELSE 

8060  T-Tf8.626897E-2+T»(3.76l99E-3-T.5.0689259E-5) 

8070  END  IF 

8080  RETURN  T 

8090  FNEND 

8100  OEF  FN8eta<T  ) 

8110  Rop-FNRhol T+. 1  ) 

8120  Ron-FNRho<T-.  1  ) 

81  30  Beta--2/<  Rop*Ron  >»<Rop-Roi>i  >/  .2 

8140  RETURN  Beta 

8150  FNEND 

8160  DEF  FNHfglT  ) 

8161  COM  /Iprop/  Ift 

8162  IF  Ift-0  THEN 

8170  Hfg-l  .374!344E+S-T.(3.3094361E-2->-T.i  .2165143  ' 

8171  END  IF 

8172  IF  Ift-I  THEN 

8173'  223  TO  373  K  CURUE  FIT  (R124) 

8176  Hfg-I59.93-.4609»T-1 . 458E-3'T '2- 1 .3715E-5»T*3 

8177  Hfg-Hfg«1000 

8178  END  IF 
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8180 

8190 

3200 

8201 

8202 

8210' 

8220 

8230 

8240 

8241 

8242 

8243' 

8245 

8247 

8248 

8250 

8250 

8270 

8280 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

8360 

8370 

8380 

8390 

8400 

84  10 

8420 

8430 

8440 

8450 

8450 

8470 

3480 

8490 

8500 

8510 

8520 

8530 

8540 

8550 

8560 

8570 

8580 

8590 

8600 

8610 

8620 

8630 

9640 

8650 

3560 

8670 

8680 

8690 

8700 

8710 

8720 

8730 

8740 

8750 


RETURN  HfQ 

FNEND 

OEF  FNPsattTc ) 

COM  /Iprop/  Ift 

IF  Ift-0  THEN 

0  TO  80  deg  F  CURVE  FIT  OF  Psat 

Tf-1  .8«Tc  +  32 

Pa-5.94552S+Tf •( . 1 S3S2082*Tf •<  1  . 4840963E-3+T f »9 .6 1 5067 1 E-6 ) ) 

Pg-Pa-14.7 

END  IF 

IF  Ift-I  THEN 

223  TO  373  K   (R-124) 

Pa-<  162.68+5.946«Tc*9.2081E-2»Tc-2+6.9359E-4.Tc"3  >•. 14504 

Pg-Pa-14.7 

ENO  IF 

IF  Pg>0  THEN       !  +  -PSI6,  — in  Hg 

Psat-Pg 

ELSE 

Psat-Pg«29. 92/14.7 

END  IF 

RETURN  Psat 

FNEND 

DEF  FNHsmooth< X  ,Bop  ,Isat  ) 

DIM  A<  5  >  ,B(5  )  ,C(5  )  ,0(5  ) 

DATA  .20526  ,.25322  .. 31 9048  .  .55322 ,.79909.1 .00258 

DATA  .745 15,. 72992  ,  .73 1 89  , . 71 225 , . 58472 ..64197 

DATA  .41092  ,.17726  ..25142  ,.54806  ,.81916.1 .0845 

DATA  .7I403..72913,.72565,.696691  , .665867  ,  .6 1 889 

READ  A(  •  )  ,B(  •  )  ,C(  •  )  ,D(  •  I 

IF  Bod<6  THEN  I-Bop 

IF  Bop-6  THEN  1-4 

IF  Bop-10  THEN  1-5 

IF  Isat-1  THEN 

Hs-EXP(  A(  I )+B( I  )»L06( X  )  ) 

ELSE 

Ha-EXP(C(  I  )-t-D(  I  )»L06(  X  )') 

END  IF 

RETURN  Ha 

FNEND 

DEF  FNPoly( X  ) 

COM  /Cply/  A(10,10>.C(I0),B(5)  .Nop .Iprnt  .Opo  ,  I Ioq  ,  I f n  ,  I  join  ,N join 

Xl-X 

Poly-8(0> 

FOR  1-1  TO  Nop 

IF    Ilog-1    THEN   X1-L06U  ) 

Poly-Poly+B<  I  >»X  1  ~  I 

NEXT    I 

IF  Ilog-1  THEN  Poly-EXP( Poly > 

RETURN  Poly 

FNEND 

SUB  Poly 

DIM  R< 10)  ,S< 10)  ,Sy( 12  >  ,Sx< 12  >  ,Xx(  100  >  ,Yy<  100  ) 

COM  /Coly/  A(  10,10)  ,C<  10)  ,B<5)  ,N  ,  Iprnt  ,  Opo,  Hog.  I  f  n  ,  I  join  ,N  join 

COM  /Xxyy/  Xp( 25  )  ,Yp( 25  ) 

FOR  1-0  TO  4 

B(I  )-0 

NEXT  I 

BEEP 

INPUT  "SELECT  (0-FILE,  I -KEYBOARD, 2-PR06RAM  )"  ,In 

Im-Im+1 

BEEP 

INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS'. Np 

IF  In-1  THEN 

BEEP 

INPUT  'ENTER  DATA  FILE  NAME',0  file* 

BEEP 
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8760  INPIT  'LIKE  TO  EXCLUDE  DATA  PAIRS  < l-Y  ,0-N >?"  ,ied 

8770  IF  ied-l  THEN 

8780  BEEP 

8790  INPUT  'ENTER  NUMBER  OP  PAIRS  TO  BE  EXCLUDED-  ,  Ipex 

8800  END  IF 

8810  ASSI6N  |Filfl  TO  D_file« 

8820  ElSE 

8830  BEEP 

8840  INPUT  "UANT  TO  CREATE  A  DATA  FILE  < I -Y  ,0-N  )? '   ,Yes 

8850  IF  Yes-l  THEN 

88E0  BEEP 

8870  INPUT  "GIVE  A  NAME  FOR  DATA  FILE"  ,D_f  l let 

8880  CREATE  BOAT  D_fileS,5 

8890  ASSIGN  IFile  TO  D_file« 

8900  END  IF 

8910  END  IF 

8920  BEEP 

8930  INPUT  "ENTER  THE  ORDER  OF  POLYNOMIAL"  ,N 

8940  FOR  1-0  TO  N»2 

8950  Sy(I)-0 

8960  Sx(I>-0 

9970  NEXT  I 

8980  IF  led-!  ANO  Irt-1  THEN 

8990  FOR  I-l  TO  Ipex 

9000  ENTER  »FileiX  ,Y 

9010  NEXT  I 

9020  ENO  IF 

9030  FOR  I-l  TO  Np 

9040  IF  In-1  THEN 

9050  IF  Opo-2  THEN  ENTER  9File;X,Y 

9060  IF  Opo<2  THEN  ENTER  IFilajY.X 

9070  IF  Opo-I  THEN  Y-Y/X 

9080  IF  Ilog-1  THEN 

9090  IF  Opo-2  THEN  Xt-X/Y 

9100  X-L061X) 

9110  IF  Opo-2  THEN  Y-LOG(Xt) 

9120  IF  Opo<2  THEN  Y-t_06<  Y  ) 

9130  END  IF 

9140  END  IF 

9150  IF  In-2  THEN 

9160  BEEP 

9170  INPUT  "ENTER  NEXT  X-Y  PAIR",X,Y 

9180  IF  Yes-l  THEN  OUTPUT  SFileiX.Y 

9190  ENO  IF 

9200  IF  In<3  THEN 

9210  Xx(I>-X 

9220  fy( I  )-Y 

9230  ELSE 

9240  x-Xp<  I-l  ) 

9250  Y-Yp( I-l  ) 

9250  END  IF 

9270  R(0)-Y 

9280  Syl0)-Sy(0>  +  Y 

9290  S( 1  >-X 

9300  Sx<  I  >-Sx<  1  )+X 

9310  FOR  J-l  TO  N 

9320  R( J  )-R( J-l  >»X 

9330  Sy( J >-Sy<J >+R< J > 

9340  NEXT  J 

9350  FOR  J-2  TO  N»2 

9360  S( J  )-5( J-l  >»* 

9370  Sx( J )-Sx(J  )+S( J  ) 

9380  NEXT  J 

9390  NEXT  I 

9400  IF  Yes-l  ANO  Iot-2  THEN 

94  10  BEEP 
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9420 

3430 

9440 

9450 

9460 

9470 

9480 

9490 

9500 

9510 

9520 

9530 

9540 

9550 

9550 

9570 

9560 

9590 

9600 

9610 

9620 

9630 

9640 

9650 

9660 

9670 

9680 

96  90 

9700 

9710 

9720 

EPANC 

9730 

9740 

9^50 

9760 

9770 

9780 

9790 

9800 

9810 

9820 

9830 

9840 

9850 

9860 

9870 

9880 

9890 

9900 

9910 

9920 

9930 

9940 

9950 

9960 

9970 

9980 

9990 

10000 

10010 

10020 

10030 

10040 

10050 

10060 


PRINT  uSIN6  "IZX.UU,""  *~Y  Pairs  were  storeo  in 

END  IF 

Sx(0  >-Np 

FOR  1-0  TO  N 

C( I  >-Sy(I  ) 

FOR  J-0  TO  N 

A(  I  ,J  >-Sx<  I+J  ) 

NEXT  J 

NEXT  I 

FOR  1-0  TO  N-l 

CALL  Divide!  I  ) 

CALL  Subtract* 1  +  1  ) 

NEXT  I 

B(N  )-C(N  )/A(N  ,N  ) 

FOR  1-0  TO  N-1 

B<  N-1 -I  )-C(N-l-I  ) 

FOR  J-0  TO  I 

B(N-I-I  >-B<  N-1 -I  )-A(N-1-I  ,N-J  >»B(N-J  ) 

NEXT  J 

B(N-1 -I  )-B(N-1 -I  >/A(N-1-I  .N-l-I  ) 

NEXT  I 

'PRINTER  IS  701 
i PRINT  B( •  ) 
'PRINTER  IS  705 

IF  Iprnt-0  THEN 

PRINT  USIN6  -12X  ,■ "EXPONENT     COEFFICIENT" "  ■ 

FOR  1-0  TO  N 

PRINT  USIN6  -15X  ,DD,5X  ,MD.70E" il  ,B(I  ) 

NEXT  I 

PRINT  •  " 

PRINT  USIN6  •12X,'"DATA  POINT      X 


i tf*  ;Nd  ,U  r  i  leS 


Y(CALCULATED)  DISCR 


FOR  1-1  TO  No 

Yc-B(0) 

FOR  J-l  TO  N 

Yc-Yc+B( J  >»Xx<  I  )" J 

NEXT  J 

0-Yy(  I  )-Yc 

PRINT  USIN6  -15X ,3D ,4X ,4(M0.SDE ,1X )' il ,Xx( I >  ,Yy( I  )  ,Yc  ,D 

NEXT  I 

END  IF 

ASSIGN  »File  TO  • 

SUBEND 

SUB  Divide(M) 

COM  /Cply/  A<  10,10 > ,C( 10 ) ,B(5 ) .N.Iprnt ,Opo  .Ilog.Ifn  .Ijoin  ,Njoin 

FOR  I-M  TO  N 

Ao-A( I  ,M  > 

FOR  J-fl  TO  N 

A(  I  ,J  )-A(  I  ,J  )/Ao 

NEXT  J 

Cl  I  )-Ct I )/Ao 

NEXT  I 

SUBEND 

SUB  Subtracts  ) 

COM  /Cply/  A(  10,10  )  ,C(  10)  .6(5  )  ,N  ,Iprnt  ,Opo  .Iloo.Ifn .Ijoin.Njoin 

FOR  I-K  TO  N 

FOR  J-K-l  TO  N 

A(  I  ,J  )-A(K-l  ,J  )-A(  I  .J  ) 

NEXT  J 

C(  I  )-C(K-l  )-C< I  ) 

NEXT  I 

SUBEND 

SUB  Plin 

COM  /Cply/  A( 10  ,10)  ,L( 10  )  ,B(5)  ,N  .Iprnt  ,Opo  , Iloo.Ifn  .Ijoin.Njoin 

COM  /Xxyy/  Xx<25  )  ,Yy<25) 

PRINTER  IS  70S 
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'0070  BEEP 

10080  INPUT  'UANT  TO  PLOT  Uo  vs  Urn'  U -Y  ,0-N  ) '  ,  Iuo 

10090  IF  Iuo-0  THEN 

10100  BEEP 

10110  INPUT  'SELECT  <0-h/h01  sane  tuba . l-h<HF  )/h< an )"  ,Irt 

10120  BEEP 

10130  INPUT  'SELECT  h/h  RATIO  <  1 "FILE ,0-COMPUTED > "  .Ihrat 

10U0  IF  Ihrat-0  THEN 

10150  BEEP 

10160  INPUT  'WHICH  Tsat  (  1 -6 . 7  ,0— 2 . 2 >'  ,laat 

10170  ENO  IF 

10180  Xnin-0 

101 90  Xnax-10 

10200  Xstep-2 

10210  IF  Irt-0  THEN 

10220  Ynin-0 

10230  Ynax-I  .4 

10240  Ystao-.2 

10250  ELSE 

10250  Ymin-0 

10270  Ynax-15 

10280  Ystep-5 

10290  END  IF 

10300  ELSE 

10310  Opo-2 

10320  Ymin-0 

10330  Ymax-12 

10340  Ystep-3 

10350  Xnin-0 

103S0  Xmax-4 

10370  Xstep-I 

10380  ENO  IF 

10390    IF    Ihrat-I     THEN 

10400    Ynun-0 

10410  Ynax-18 

10420  Ysteo-3 

10430  Xr-un-0 

10440  Xnax-9 

10450  Xstep-2 

10450  END  IF 

10470  BEEP 

10480   PRINT    "INiSPldP    2300.2200.8300.58001- 

10490   PRINT    'SC   0.100.0.100;TL    2.0i' 

10500    Sf  x-l00/(Xmax-Xi>un  > 

10510    Sfy-l00/( Ymax-Yeun  ) 

10520    PRINT    "PU   0.0    PO' 

10530  FOR  Xa-Xnin  TO  Xmax  STEP  Xatep 

10540  X-( Xa-Xnin  )»Sfx 

10550  PRINT  'PA'  iX  .'  .01  XTi  ' 

10560  NEXT  Xa 

10570  PRINT  'PA  100,0iPUi  ' 

10580  PRINT  'PU  PA  0.0  PD' 

10590    FOR    Ya-Yeun    TO    Ynax    STEP    Ystep 

10600    Y-(  Ya-Ymn)»Sfy 

10610  PRINT  'PA  0," iY  ,'YT" 

10520  NEXT  Ya 

10630  PRINT  "PA  0.100  TL  0  2" 

10640  FOR  Xa-Xnin  TO  Xnax  STEP  Xstep 

10650  X-( Xa-Xnin  >»Sfx 

10560  PRINT  'PA' tX  ,'  ,  1001  XT' 

10670  NEXT  Xa 

10680  PRINT  'PA  100.100  PU  PA  100.0  PD' 

10690   FOR   Ya-Ymn    TO    Ymax    STEP    Ystep 

10700  Y-( Ya-Ymin  )»Sfy 

10710  PRINT  'PD  PA  I00,',Y.'YT' 

10720  NEXT  Ya 
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0730  PRINT  'PA  100,100  PU" 

07*0  PRINT  "PA  B, -I  Sft  1.5.2" 

0750  FOR  Xa-Xmin  TO  Xnax  STEP  Xstop 

07S0  X-( Xa-Xnin >*Sfx 

0770  PRINT  "PA" iX  ,"  .01  ' 

0780  IF  Iuo-0  THEN  PRINT  "CP  -2  ,-t iLB' iXat ' " 

0790  IF  Iuo-I  THEN  PRINT  *CP  -  I . 5  ,- 1  iLB" ; Xa i  " " 

0800  NEXT  Xa 

0810  PRINT  -PU  PA  0  ,0" 

0820   FOR   Ya-Ynin   TO   Yrtax    STEP    Ystep 

0830    IF    ABS<YaKl.E-5    THEN    Ya-0 

0840  Y-< Ya-Ynin  )»Sfy 

0850  PRINT  -PA  0,"  jY,"" 

0860  IF  Iuo-0  THEN  PRINT  "CP  -4  ,- . 25 ; LB- i Ya  ;  "  ■ 

0870  IF  Iuo-I  THEN  PRINT  "CP  -3  ,-  ,25tLB"  i  Yat " " 

0880  NEXT  Ya 

0890  Xlabel$-"Oil  Percent" 

0900  IF  Iuo-0  THEN 

0910  IF  Irt-0  THEN 

0920  Ylabel»-"h/h0X" 

0930  ELSE 

0940  Ylabel»-"h/hsmooth" 

0950  END  IF 

0960  PRINT  "SR  l.5,2(PU  PA  50,-10  CP " i -LEN( X  label* >/2 i "0tLB" i XlabelS i " " 

0970  PRINT  -PA  -11,50  CP  0,";-LEN(YlabelS  )/2»S/S;"DI  0  , 1 iLB" i YlabelSi" " 

0980  PRINT  "CP  0,0" 

0990  ELSE 

1000  PRINT  "SP0iSP2" 

1010  PRINT  "SR  1 .2  ,2.4iPU  PA  -8  ,35:01  0,1 iLBUiPR  1  ,0.5iLBoiPR  -1  ,0.5; LB  ! kU/n 

1020  PRINT  "PR  -1  ,0.5iSR  1,1.5iLB2;SR  1.5,2iPR  .5,.5iLB.iPR  .5,0iLBK>" 

1030  PRINT  'PA  42,-10iOI  1 ,0;LBViPR  .4,-liLButPR  I , .SiLB(n/s )' 

1040  PRINT  "SP0;SP1 " 

1050  END  IF 

1060  Ipn-0 

1070  BEEP 

1080  INPUT  'WANT  TO  PLOT  DATA  FROM  A  FILE  (  1 -Y ,0-N >?"  ,Okp 

1090  Icn-0 

1  100  IF  Dkp-I  THEN 

1110  BEEP 

1120  INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE"  ,D_f i let 

1 130'  IF  Iuo-0  THEN 

1 140  BEEP 

1150  INPUT  "SELECT  (0-LINEAR,  1 -LOS( X  ,  Y  )"  ,  I  log 

11 S01  END  IF 

1170  ASSISN  eFile  TO  D_file$ 

1 180  BEEP 

1190  INPUT  "ENTER  THE  BESINNIN6  RUN  NUMBER", Md 

1200  BEEP 

1210  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED"  .Npairs 

1220  IF  Iuo-0  AND  Ihrat-0  THEN 

1230  BEEP 

1240  INPUT  "ENTER  DESIRED  HEAT  FLUX'.Q 

1250  END  IF 

1260  BEEP 

1270  PRINTER  IS  1 

1280  PRINT  USIN6 

1290  PRINT  USING 

1300  PRINT  USIN6 

1310  PRINT  USIN6 

1320  PRINT  USING 

1330  PRINT  USIN6 

1340  INPUT  Syn 

1350  PRINTER  IS  705 

1360  PRINT  "PU  DT 

1370    IF    Syn-I    THEN   PRINT    "SM 
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4X 

""Se 

lee t  a  synbol  : " " " 

4X 

" "  1 

Star    2   Pius  sign" "  " 

4X 

""3 

Circ le  4   Square  "  "  " 

4X 

""5 

Ronous " "  " 

4X 

""6 

Right -s ide-up  triangle" 

4X 

•"7 

Up-side-doun  triangle"" 

1380  IF  Syw-2  THEN  PRINT  "SM+" 

1390  IF  Syn-3  Tn£N  rniNi   sno 

1400  Nn-4 

1410  IF  I  log- 1  THEN  Nn-1 

1420  IF  Hd>1  THEN 

1430  FOR  1-1  TO  (Md-1  ) 

1440  ENTER  CFileiXa  ,Va 

MS0  NEXT  I 

1460  END  IF 

1470  IF  Ihrat-0  THEN 

1480  01-0 

1490  IF  I  loo-'  THEN  Q-LOS(Q) 

1500  ENO  IF 

1510  FOR  1-1  TO  Npairs 

1520  IF  Iuo-0  AND  Ihrat-0  THEN 

1530  ENTER  9FileiXa.B<»  ) 

1540  Ya-B(0) 

1550  FOR  K-l  TO  Nn 

1550  Ya-Ya+B(K  )»0"K 

1570  NEXT  K 

1580  END  IF 

1590  IF  Iuo-l  OR  Ihrat-1  THEN 

1500  ENTER  9FileiXa,Ya 

1610  IF  Iuo-l  THEN  Ya-Ya/1000 

1620  END  IF 

1630  IF  Iuo-0  AND  Ihrat-0  THEN 

1640  IF  I  log- 1  THEN  Ya-EXP( Ya  ) 

1650  IF  Ilog-0  THEN  Ya-Ql/Ya 

1660  IF  Irt-0  THEN 

1670  IF  Xa-0  THEN 

1680  Yo-Ya 

1690  Ya-1 

1700  ELSE 

1710  Ya-Ya/Yo 

1720  END  IF 

1730  ELSE 

1740  Hsrt-FNHanooth(  Q  ,Xa  ,Isat  ) 

1750  Ya-Ya/Hsn 

1760  END  IF 

1770  END  IF 

1780  Xx(  1-1  >-Xa 

1790  Yy(  1-1  )-Ya 

1800  X-(Xa-Xnin  )»Sfx 

1810  Y-( Ya-Ynin  >»Sfy 

1820  IF  Sym>3  THEN  PRINT  "Sfl" 

1830  IF  Syn<4  THEN  PRINT  "SR  1.4,2.4" 

1840  PRINT  "PA"  ,X,Y  ,*" 

1850  IF  Sym>3  THEN  PRINT  "SR  1.2.1.6" 

1860  IF  Sym-4  THEN  PRINT  "UC2  .4  ,99  .0  ,-8  ,-4  .0  ,0  ,8  .4  ,0  .  i  " 

1870  IF  Sym-S  THEN  PRINT  "UC3  ,0  ,99  ,-3  .-6  .-3  ,6  ,3  ,6  .3  .-6 i ' 

1880  IF  Syn-6  THEN  PRINT  "UC0  .5  .  3  ,99  ,3  ,-8  ,-6  .0  ,3  ,9  i  " 

1890  IF  Syn-7  THEN  PRINT  *UC0  ,-S . 3  .99  ,-3  .8  ,6  .0  .-3  . -8 i  " 

1900  NEXT  I 

1910  BEEP 

1920  ASSI6N  »File  TO  • 

1930  END  IF 

1940  PRINT  "PU  SM' 

1950  BEEP 

I960  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  (  1 -Y  ,0-N  P "  ,Okp 

1970  IF  Okp-1  THEN 

1980  BEEP 

1990  PRINTER  IS  I 

2000  PRINT  USIN6  "4X, "-Select  line  tyoe:""" 

2010  PRINT  USIN6  "6X,""0     Solid  line""' 

2020  PRINT  USIN6  "SX  ." '  1      Dashed"'" 

2030  PRINT  USIN6  "6X.""2.,,5  Longer  line  -  aasn"** 
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12040  INPUT  Ipn 

itDsd  PRINTER  la  70S 

12060  BEEP 

12070  INPUT  "SELECT  ( 0-LINEAR  ,  I -LOG< X  ,Y  )  )"  ,  I  log 

12080  Iprnt-I 

12090  CALL  Poly 

12100  IF  Iuo-I  THEN 

12110  BEEP 

12120    INPUT    "DESIRE    TO   SET    X   Lower   and   Upper    Limit    (  1 -Y  ,0-N >? '  .  Ix  1  in 

12130    IF    IxUm-0    THEN 

12U0    xii-0 

12150   Xul-7 

12150   END    IF 

12170  IF  Ixlim-I  THEN 

12180  BEEP 

12190  INPUT  "ENTER  X  Lower  Limit", Xll 

12200  BEEP 

12210  INPUT  "ENTER  X  Upper  Limit", Xul 

12220  END  IF 

12230  END  IF 

12240  FOR  Xe-Xll  TO  Xul  STEP  Xstep/2S 

12250  Icn-Icn+1 

122S0  Ya-FNPoly(Xa) 

12270  IF  luo-1  THEN  Ya-Ya/1000 

12280  Y-< Ya-Ymin >»Sfy 

1 2290  X-< Xa-xmin  )»Sfx 

12300  IF  Y<0  THEN  Y-0 

12310  IF  Y>100  THEN  SOTO  12410 


12320 

Pu-0 

12330 

IF  Ipn-1  THEN  Idf-Icn  MOD  2 

12340 

IF  Ipn-2  THEN  Idf-Icn  MOD  4 

12350 

IF  Ipn-3  THEN  Idf-Icn  MOD  8 

12350 

IF  Ipn-4  THEN  Idf-Icn  MOD  16 

12370 

IF  Ipn-5  THEN  Idf-Icn  MOD  32 

12380 

IF  Idf-1  THEN  Pu-1 

12390 

IF  Pu-0  THEN  PRINT  "PA'.X.Y, 

■PD" 

12400 

IF  Pu-I  THEN  PRINT  "PA"  ,X  ,Y,' 

■PU" 

12410 

NEXT  Xa 

12420 

PRINT  "PU" 

12430 

60T0  1 1070 

12440 

END  IF 

12450 

BEEP 

I  2460 

INPUT  "UlnNT  TO  QUIT  (  l-Y  ,0-N )?"  ,1 

12470 

IF  Iouit-1  THEN  12490 

12480 

60T0  1 1070 

12490 

PRINT  "PU  SP0" 

12500 

SUBEND 

12510 

SUB  Stats 

12520 

PRINTER  IS  701 

12530 

J-0 

12540 

K-0 

12550 

BEEP 

125S0 

IF  Iplot-I  THEN  ASSIGN  SFile 

TO  P 

12570 

BEEP 

12580 

INPUT  "LAST  RUN  No't 0-OUIT  )  " 

,Nn 

12590 

IF  Nn-0  THEN  12950 

12600 

Nn-Nn-J 

12610 

Sx-0 

12620 

Sy-0 

12630 

Sz-0 

12640 

Sxa-0 

12650 

Sys-0 

12660 

Szs-0 

12670 

FOR  I-l  TO  Nn 

12680 

J-J+l 

12690 

ENTER  •FileiQ.T 

file$ 
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;-00  m-q/t 
2718  Sx-Sxt-Q 

2720  Sxa-SxsfQ-2 

2730  Sy-Sy+T 

2740  SyS-Sye*T"2 

2750  Sz-Sz+H 

27B0  Sz9-Szs*H*2 

2770  NEXT  I 

2780  Qave-Sx/Nn 

2790  Tave-Sy/Nn 

2800  Have-Sz/Nn 

2810  Sdevq-SOR( ABS<  (Nn»Sxs-Sx"2  >/(Nn»(Nn-l  )  )  )  > 

2820  Sdevt-SQR( ABSl  (Nn»Sys-Sy"2  )/(Nn»(Nn-1  )  ) )  ) 

2830  Sdevh-SQR( ABS<  (Nn»Szs-Sz"2  >/(Nn«(Nn-l  )  )  )  > 

2840  Sh-l00«Sdevn/Have 

2850  Sq-100«SdevQ/Qave 

2860  St-100«Sdevt/Tave 

2870  IF  K-l  THEN  12930 

2880  PRINT 

2890  PRINT  USIN6  'MX, ""DATA  FILE  :  "  "  ,  I  4A"  ;F  i  leS 

2900  PRINT 

2910  PRINT  USIN6  'MX.'-RUN  Htube     SdevH    Odp        SdevO    Thetab  SdevT'" 

2920  K-l 

2930   PRINT   USIN6    ■ 11 X  ,DD  ,2( 2X  ,0 . 3DE  ,  I  X  ,30 .20  )  ,2X  ,OD  .3D  .  1 X  ,3D . 2D" i J .Have ,Sh  .Qave 

Sq  .Tave  .St 

2940  SOTO  12570 

2950  ASSIGN  •File!  TO  • 

2960  PRINTER  IS  1 

2970  SUBEND 

2980  SUB  Coef 

2990  COM  /Cply/  A<  1  0 , 1  0 ) ,C(  1 0 > ,B< 5 ) ,N .Iprnt ,0po  ,  I  log  ,  I f n  ,  I  join  .Njoin 

3000  BEEP 

3010  INPUT  "GIVE  A  NAME  FOR  CROSS-PLOT  FILE"  ,CpfS 

3020  CREATE  BOAT  CpfS  .2 

3030  ASSIGN  OFile  TO  CpfS 

3040  BEEP 

3050  INPUT  'SELECT  ( 0-lINEAR  ,  1 -L06( X ,Y  )  ) •  , I  log 

3060  BEEP 

3070  INPUT  "ENTER  OIL  PERCENT  i - 1 -STOP  ) "  ,Bop 

3080  IF  Bop<0  THEN  13120 

3090  CALL  Poly 

3100  OUTPUT  CFileiBop  ,B( • > 

31  10  60T0  13060 

3120  ASSIGN  eFile  TO  • 

3 '30  SUBEND 

3140  SUB  UilsonfCf  ,Ci  > 

3150  COM  /Uil/  D2  ,Di  ,Do  .L.Lu  .Kcu 

3160  OIM  Emf  < 12  ) 

31701  ULISON  PLOT  SUBROUTINE  OETERMINE  CF  AND  CI 

3180  BEEP 

3190  INPUT  "ENTER  DATA  FILE  NAME',File» 

3200  BEEP 

3210  PRINTER  IS  I 

3220  PRINT  USIN6  '4X, "Select  option:*** 

3230  PRINT  USIN6  "4X,'-   0  Wary  Cf  and   Cl**' 

3240  PRINT  USIN6  '4X,--   1  Fix   Cf  Wary  Cl*** 

3250  PRINT  USIN6  -4X,--   2  vary  Cf  Fix   Cl"*" 

3260  INPUT  -ENTER  OPTION"  .  Ic f u 

3270  PRINTER  IS  701 

3280  IF  Icf  ix-0  THEN  13320 

3290  IF  lcfix>0  THEN  BEEP 

3300  IF  Icfix-1  THEN  INPUT  'ENTER  Cf".C»f 

3310  IF  Icfix-2  THEN  INPUT  "ENTER  CI'.Ci 

3320  PRINTER  IS  1 

3330  INPUT  "Uant  To  Vary  Coef f '<  I -Y  ,0-N  )*  , Iccoef 
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13340  IF  Iccoef-1  THEN  INPUT  "ENTER  COEFF".R 

13358  PRINTER  IS  7a  1 

13360  IF  Icfix-0  OR  Icfix-2  THEN  Cfa-.004 

13370  IF  Icfix-1  THEN  Cfa-Csf 

13380  Cia-Cl 

13390  Xn-.8 

13400  Fr-.3 

13410  Jj-0 

13420  Rr-3. 

13430  IF  Iccoef-1  THEN  Rr-R 

13440  PRINTER  IS  1 

134S0  PRINT  Do.Di  .Kcu 

134S0  ASSIGN  SFile  TO  File* 

13470  ENTER  »F i la iNrun .OateS  ,Ldtc 1  ,Ldt c2  , 1 t 1 

13480  Ry-0o«L06(Do/0i  )/C2»Kcu  ) 

13490  Sx-0 

13500  Sy-0 

13510  Sxy-0 

13520  Sx2-0 

13530  Sy2-0 

13540  FOR  1-1  TO  Nrun 

13550  ENTER  «Fi le iBoo .ToldS  ,Enf ( •  )  ,Fms 

13560'  CONVERT  EMF'S  TO  TEMPERATURE 

13570  FOR  J-1  TO  5 

13580   T(J  >-FNTvsv<Emf( J  )  ) 

13590   NEXT    J 

1 3600   Tsat-( T( 1  >+T<2  >  >•  .5 

13610    Tavg-T(5) 

13620   6rad-37.98S3+. 104388»Tavg 

13630    Tdrop-Emf<7  )•!  .E+6/( 10.»Grad> 

13640  Tavgc-T<5 >-Tdrop».5 

13650  IF  ABS<Tavg-Tavgc >>.01  THEN 

13660  Tavg-<  Tavgt-Tavgc  )•  .5 

13670  60T0  13620 

13680  END  IF 

136901 

137001  Compute  properties  of  water 

13710  Kw-FNKu(Tavg  ) 

13720  Muwa-FNMuw( Tavg  ) 

13730  Cpu-FNCpwCTavg  ) 

13740  Pru-FNPrw(Tavg  ) 

13750  Rhou-FNRhow< Tavg  ) 

13760' 

137701    Compute   properties    of    Freon-114 

13780   Lmtd-Tdrop/L06( ( T( 5  )-Tsat  )/ ( T( 5  >-Tdrop-Tsat  )  ) 

13790    IF    Jj-0   THEN 

13800  Tu-Tsat*Fr»Lntd 

13810  Thetab-Ty-Tsat 

13820  Jj-I 

13830  END  IF 

13840  Tf-(T«+Tsat  >».5 

13850  Rho-FNRho(Tf  ) 

13860  Mu-FNMu(Tf  ) 

13870  K-FNK(Tf  ) 

13880   Cp-FNCpUf  ) 

13890  Beta-FN8eta<Tf  ) 

13900  Hfg-FNHfg( Tsat  ) 

13910  Ni-Mu/Rho 

13920  Alpha-K/(Rho»Cp  ) 

13930  Pr-Ni/Alpha 

139401 

13950>  Analysis 

13960!  COMPUTE  MOOT 

13970  A-PI«(Do*2-Di*2)/4 

13980  P-PI.Oo 

13990  Mdot-3.9657E-3+Fms«(3.61955E-3-Fms»(8.82006E-6-Fma»( 1 . 23688E-7-Fms«4 . 31 897 
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E-l 0  > >> 

4000  o-Mdot«cpw«Torop 

401 0 i  COMPUTE  NATURAL-CONUECTIUE  HEAT-TRANSFER  COEFFICIENT 

4020'  FOR  UNENHANCED  END(S) 

4030  Hbar-190 

4040  Fe-(Hbar«P/(Kcu»A  )  >* .5»Lu 

4050   Tann-FNTanMFe  ) 

4060  Theta»Thetad»Tanh/Fe 

4070  Xx-(  9.81  ♦8eta»Thetao»Do"3»Tann/( F««Ni 'Alpha  )  )* . 1  ESSE 7 

4080  Yy-( !♦( .559/Pr  r<  9/16  )>"< 8/27  ) 

4090  Hbarc-K/Do»( .6+.387«Xx/Yy  )"2 

4100  IF  ABS( (Hbar-Hbarc  )/Hbar >>  301  THEN 

4110  Hbar-< Hbar+Hbarc  )•  .5 

4120  GOTO  14040 

4130  END  IF 

4  140> 

4150'  COMPUTE  HEAT  LOSS  RATE  THR0U6H  UNENHANCED  ENDS 

41 S0  Ql-(Hbar»P»Kcu»A)".5»Thetao»Tanh 

4170  Qc-Q-2»Q1 

4180  As-PI»D2»L 

4190'  COMPUTE  ACTUAL  HEAT  FLUX 

4200  QdD-Oc/As 

4210  IF  Icfix-a  OR  Icftx>l  THEN  Csf-1 /Cf " (  1  . /Rr  ) 

4220  Thetab-Csf /Co»Hfg»<Odo/(Mu»Hfo  )•<  . 0  I 4/( 9 . 8 1 »Rho  )  >* .5  ) " ( 1 /Rr )»Pr*1 .7 

4230  Ho-Ooo/Thetab 

4240  Onega-Ho/Cf 

4250  Uo-Q/(PI«Do»L«Lntd> 

42S0   Vu-Mdot/<  Rhow»PI»Di*Z/4  ) 

4270   Rey"Rhou»Uu«Oi/Muyia 

4280    Tyi-TufQ»Rw/(PI»Oo«L ) 

4290    Gama-Kui/Oi»Rew* . 8»Pru"(  1/3.  >♦< Muwa/FNMuu( Tui  )  )* . 14 

4300'  PRINTER  IS  1 

4310  Yy-(  1 ./Uo-Ru)»Onega 

4320  Xu-Omega«Do/(Gama»Di  ) 

4330  Sx-Sx-t-Xy 

4340  Sy-5y*Yw 

4350  Sxy-5xy+Yy»Xw 

4360  Sx2-Sx2+Xw»Xy 

4370  Sy2-Sy2+Yy«Yy 

4380  NEXT  I 

4390  ASSIGN  "File  TO  • 

4400  M«( Sx»Sy-Nrun*Sxy  )/( Sx •Sx-Nrun«Sx2  ) 

4410  C-<Sy-Sx»M)/Nrun 

4420  IF  Icflx-0  OR  Icfix-3  OR  Icfix-4  THEN 

4430  ClC-l/M 

4440  Cf C- I /C 

4450  END  IF 

4460  IF  Icfix-1  THEN 

4470  Cic-1 /M 

4480  Cfc-Cf 

4490  END  IF 

4500  IF  Icfix-2  THEN 

4510  Cic-Ci 

4520  Cfc-1/C 

4530  END  IF 

4540  IF  A8S( (Ci-Cic >/Cic  )>.00l  OR  ABSl ( Cf-Cf c  )/Cf c  )>  .001  THEN 

4550  Ci-(Cl+Clc  )».5 

4560  Cf-( Cf+Cfc  )».5 

4570  PRINTER  IS  I 

4530  PRINT  USIN6  '2X,"    Csf  •  " "  .MZ.30E  ,2X  ,' -  Cx  -  ' '  ,M2 . 3DE- i Csf .Ci 

4590  PRINTER  IS  701 

4600  60T0  13460 

4610  END  IF 

4620  PRINT 

4630  PRINTER  IS  701 

4640  PRINT  USIN6  "Z3X,--   Cf  Ci"* 
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14653  PRINT  USING  "3X  ,"  "ASSUMED      " " ,MZ . 3DE  ,  3x  ,MZ . 30E ' ; C ca  ,  Z i a 

14660  PRINT  USIN6  "8X  ,"  "CALCULATED    " "  ,P1Z  .  3DE  ,3X  ,n:  .  3DE  "  ;  Cs  f  ,  Z  : 

14570  PRINT 

14680  Sum2-Sy2-2«M»5xy-2«C»Sy+M"2»Sx2+2»M«C»Sx+Nrun«C*2 

14690  PRINT  USIN6  *10X,"Sun  of  Squares  -  " " ,Z . 3DE" s Sum2 

14700  PRINT  USIN6  '10X, "Coefficient  -  ".D.DDD'iRr 

14710  SUBEND 

14720  DEF  FNMuulT  ) 

14730  A-247.8/(T+l33.l5  > 

14740  riu-2.4E-5»l0"A 

14750  RETURN  Mu 

14760  FNEND 

14770  DEF  FNCpw(T) 

14780  Cpu-4.21 l208SB-T.<2.2S826E-3-T.<4.42361E-5+2.7l428E-7.T  )  ) 

14790  RETURN  Cpu*1000 

14800  FNEND 

14810  DEF  FNRhoy(T  ) 

14820  Ro-999.52946+T»<  . 01 269-T«( 5 . 48251 3E-3-T* 1 .234147E-5)  ) 

14830  RETURN  Ro 

14840  FNEND 

14850  DEF  FNPrulT  ) 

14860    Prw-FNCowlT  >»FNMuw< T  >/FNKw< T  ) 

14870  RETURN  Pru 

14880  FNEND 

14890  DEF  FNKw(T  ) 

14900  X-(T+273.IS)/273.  15 

14910  Kw— .92247+X»(2.8395-X»(  I .8007-X»(  .52577- . 07344* X  )  )  ) 

14920  RETURN  Kw 

14930  FNENO 

14940  SUB  Plot 

14950   COn    /Cply/    A<  10,10)  ,C( 10)  ,B<5> , Nop ,Iprnt  ,0po  .Ilog.Ifn  .Ijoin.Njoin 

14960   DIM   Bs(3) 

14970  INTE6ER  Ii 

14980  PRINTER  IS  1 

14990  Idv-0 

15000  8EEP 

15010  INPUT  "LIKE  DEFAULT  UALUES  FOR  PLOT  (  1 -Y ,0-N  >?• , Idv 

15020  Ooo-0 

15030  BEEP 

15040  PRINT  USIN6  "4X, ""Select  Option:"" 

15050  PRINT  USING  *SX."0   q  versus  Oeita-T""" 

15060  PRINT  USIN6  "SX,"1   h  versus  delta-T"" 

15070  PRINT  USING  "6X,""2   h  versus  q""" 

15080  INPUT  Opo 

15090  BEEP 

15100  INPUT  "SELECT  UNITS  ( 0-SI  .  1 -ENGLISH )". Iun 

151 10  PRINTER  IS  705 

15120  IF  IdvOI  THEN 

15130  BEEP 

15140  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS", C* 

15150  BEEP 

15160  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS", Cy 

15170  BEEP 

15180  INPUT  "ENTER  MIN  X-UALUE  (MULTIPLE  OF  10)",Xnin 

15190  BEEP 

15200  INPUT  "ENTER  MIN  Y-VALUE  (MULTIPLE  OF  10)",Yrun 

15210  ELSE 

15220  IF  Opo-0  THEN 

15230  Cy-3 

15240  Cx-3 

15250  Xmin-. 1 

15260  Ymn-100 

15270  END  IF 

15280  IF  Opo-1  THEN 

15290  Cy-2 

15300  Cx-2 
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1 531  a 

15320 
15330 
153*0 
15350 
I53S0 
15370 
15380 
15390 
15400 
15410 
15420 
15430 

15440 

15450 
15460 
15470 
15480 
15490 
15500 
15510 
15520 
15530 
15540 
15550 
15550 
15570 
15580 
15590 
15500 
I  55  10 
15520 
15530 
15540 
15550 
155S0 
15570 
15680 
15690 
15700 
15710 
15720 
15730 
15740 
15750 
15760 
15770 
15780 
15790 
15800 
15810 
15820 
15830 
15840 
^5850 
15860 
15870 
15880 
15890 
15900 
15910 
15920 
15930 
15940 
15950 
15960 


Xmn-  . I 
ruin- 100 

END  IF 

IF  Opo-2  THEN 
IF  Iun-0  THEN 
Cy-3 
Cx-2 

Xnin-I 000 
Ymin-1 00 
ELSE 
Cy-3 
Cx-3 

Xnun-1 00 
Ynin-I 0 
END    IF 
END    IF 
END    IF 
8EEP 

PRINT  " IN; SP 1  : IP  2300  ,2200  .8300  .6800 i " 
PRINT  "SC  0.I00.0.100ITL  2  ,01 ■ 
Sfx-100/Cx 
Sfy-100/Cy 
BEEP 

INPUT  -WANT  TO  BY-PASS  CASE''  (1-Y.0-N)" 
IF  Ibyp-I  THEN  16790 
PRINT  'PU  0.0  PD" 
Nn-9 

FOR  1-1  TO  Cx+1 
Xat-X«un«  10"(  1-1  > 
IF  I-Cx+1  THEN  Nn-1 
FOR  J-1  TO  Nn 
IF  J-1  THEN  PRINT  "TL  2  0" 
IF  J-2  THEN  PRINT  "TL  1  0" 
Xa-Xat»J 

X-lST<  Xa/Xrun  "'Sf x 
PRINT  "PA" iX  ,*  .0;  XT  i  " 
NEXT  J 
NEXT  I 

PRINT  'PA  100  ,0iPUi  " 
PRINT  "PU  PA  0 ,0  PD* 
Nn-9 

FOR  I-l  TO  Cy+1 
Yat-Ymn«10*<  1-1  ) 
IF  I-Cy+1  THEN  Nn- I 
FOR  J-1  TO  Nn 
IF  J-1  THEN  PRINT  "TL  2  0' 
IF  J-2  THEN  PRINT  "TL  1  0" 
Ya-Yat«J 

Y-L6T( Ya/Ynin)«Sfy 
PRINT  "PA  J,' |Y  ,"YT" 
NEXT  J 
NEXT  I 

PRINT  'PA  0,100  TL  0  2' 
Nn-9 

FOR  1-1  TO  Cx+I 
Xat-Xnin»10*(  1-1  ) 
IF  I-Cx-H  THEN  Nn-I 
FOR  J-1  TO  Nn 
IF  J-1  THEN  PRINT  "TL  0  2" 
IF  J>1  THEN  PRINT  "TL  0  I" 
Xa-Xat«J 

X-L6T(Xa/Xnln  )»Sfx 
PRINT  "PA"iX  ,"  .1001  XT- 
NEXT  J 
NEXT  I 
PRINT  'PA  100,100  PU  PA  100,0  PD" 


■  Ibyp 
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15970 

Nn-9 

15330 

FOR  I- 

1  TO  Cyfl 

15990 

Yat-Y«in«10"( I- 

16000 

IF  I-C 

y+1  THEN  Nn-1 

16010 

FOR  J- 

1  TO  Nn 

16020 

IF  J-l 

THEN  PRINT  *TL 

0  2- 

16030 

IF  J>1 

THEN  PRINT  'TL 

0  1" 

16040 

Ya-Yat 

♦  J 

16050 

Y»l6T( Ya/Yrun  )«Sfy 

16060 

PRINT 

■PD  PA  100.'  .Y  ,' 

■YT" 

16070 

NEXT  J 

16080 

NEXT  I 

16090 

PRINT 

'PA  100,100  PU" 

16100 

PRINT 

•PA  0,-2  SR  1 .5 

,2- 

161  10 

Ii-LGTUmn ) 

16120 

FOR  I- 

1  TO  Cx+1 

16130 

Xa-Xmln« I0*< 1-1  ) 

16140 

X"LST( Xa/Xmin >»Sfx 

16150 

PRINT 

■PA-iX  .*  ,0i" 

16160 

IF  U>i 

■0  THEN  PRINT  "CP  -2 

!iLB10iPR  -2 

,2i 

LB"; 

Ell 

•  • 

16170 

IF  Ii<l 

9  THEN  PRINT  "CF 

'  -2, 

-2i 

LBI0 

iPR  0,2 

iLB 

Mi 

1  " ' 

16180 

Ii-I i+ 

1 

16190 

NEXT  I 

16200 

PRINT 

"PU  PA  0,0" 

16210 

Ii-L5T( Ymin ) 

16220 

Y 1 0-1 0 

16230 

FOR  I- 

1  TO  Cy+1 

16240 

Ya»Ynin»10"( 1-1 ) 

16250 

Y-L6T( Ya/Ynin)»Sf y 

16260 

PRINT 

•PA  0  ,"  iY  ,•" 

16270 

PRINT 

'CP  -4,-.25tLB1« 

hPR 

-2 

,2.LB 

*  «  I  A  !  "  " 

16280 

Il-Il+ 

1 

16290 

NEXT  I 

16300 

BEEP 

16310 

INPUT  ' 

•UANT  USE  DEFAULT  LABELS  (1 

■Y  ,0-N  ) 

?• 

Idl 

16320 

IF  IdlOl  THEN 

16330 

BEEP 

16340 

INPUT 

•ENTER  X-LABEL- 

,Xlab 

ell 

16350 

BEEP 

16360 

INPUT 

'ENTER  Y-LABEL' 

.YlaDelS 

15370 

END  IF 

16380 

IF  Ooo' 

C2  THEN 

16390 

PRINT  ' 

'SR  1  ,2iPU  PA  48 

1,-14 

1  " 

16400 

PRINT  ' 

'L8(TtPR  -1  .6,3 

PO  PR  1 

.2.0 

PUlPR 

•  5, 

-4tLBuoiPR  .5,1 i' 

16410 

PRINT  ' 

'LB-TiPR  .5,-1 iLBsat 

iPF 

!  .5. 

1  i' 

16420 

IF  Iun> 

■0  THEN 

16430 

PRINT  ' 

'LB)  <K>" 

16440 

ELSE 

16450 

PRINT  ' 

'LB)  (FT 

16460 

END  IF 

16470 

END  IF 

15480 

IF  Opo-2  THEN 

15490 

IF  Iun-0  THEN 

16500 

PRINT  ' 

'SR  1  .5.2iPU  PA 

40.- 

I4i 

LBg 

[U/niSR 

1  , 

1 .SiPR 

0.5.1 iL62iSR 

0.5.-1 

iLB  )* 

16510 

ELSE 

16520 

PRINT  ' 

'SR  1  .5.2iPU  PA 

34.- 

14, 

LBq 

(Btu/hnPR 

.5, 

.5. 

LB.iPR  .5.-. 

16530 

PRINT  ' 

'LBftiPR  .5.1 iSF 

!  1  ,1 

.5. 

LB2;' 

5R  1  .5, 

2iPR  .5 

,-  1 

iLB)i" 

15540 

END  IF 

16550 

END  IF 

16560 

IF  Opo-0  THEN 

16570 

IF  Inn- 

•0  THEN 

15580 

PRINT  ' 

'SR  1  .S,2iPU  PA 

-12. 

401 

01  0 

,1 iLBq 

(U/ 

niPR 

-1 

,0.SiSR  1  ,1 . 

.5.2(PR  l . .Si 

iLB  >" 

16590 

ELSE 

16600 

PRINT  ' 

'SR  1  .5,2iPU  PA 

-12. 

32i 

DI  0 

.1  1  LBq 

(Bt 

u/hr 

iPF 

1  -.5.. St  LB. 1 

1 .5.2.PR 


5. 


StLB2iSR    1 


iPR    .5.  .Si 


135 


66'0  =  RINT  "LB'tiSR  i  ,; .5iPR  -! ,.5il32;PR  I,.SiSR  l.5,2»LB)" 

ss:a  end  if 

5630  END  IF 

5640  IF  Opo>0  THEN 

6S50  IF  Iun-0  THEN 

6  582   PRINT  " SR  I  .5,2|PU  PA  - 1  2. 38 1 01  0.1 iLBh  (U/riPR  -I  ,.5iSR  I  ,I.S|LB2|SR  I 

5.2iPR  .  5  , .  5  i  " 

6G70  PRINT  *SR  1.2,2.4iPU  PA  -12  .37 1 01  0.1il_BhiPR  I  ,0.SiLBo»PR  -1.0.5tLB  (U/m 

SB30  PRINT  "PR  -l,.5iSR  1.1.5  i  LB2  i  SR  I.S.2iPR  . S , . S ;LB. iPR  .S.0ILBK)" 

6690  ELSE 

S700  PRINT  "SR  l.5.2:PU  PA  -12,28iDI  B.liLBh  (Btu/hnPR  -.5 ,  .SiLB. iPR  .5..5T 

S710  PRINT  "LBftiPR  -i,.5iSR  1.1.5iLB2iSR  1.5,2iPR  .5,.5iLB.iPR  .5,.SiLBF>" 

5720  ENO  IF 

6730  END  IF 

6740  IF  Idl-0  THEN 

5750  PRINT  "SR  1.5,2;PU  PA  50,-16  CP" i -LEN(  X  labels  )/2 : "0 ;LB" i XlabelS i * " 

5760  PRINT  -PA  -14,50  CP  0.- i-LEN( Y labels >/2»5/ 6 i "01  0  , 1 iLB" : YlabelS i " " 

6770  PRINT  "CP  0.0  DI" 

6780  END  IF 

6790  Ipn-0 

6800  Xll-l .E*B 

6810  Xul  — l  .E+6 

6820  Icn-0 

5830  Ifn-0 

6840  Ijoin-I 

6850  BEEP 

S860  INPUT  "WANT  TO  PLOT  OATA  FROM  A  FILE  C  1 -Y ,0-N >?' ,0k 

6870  IF  Ok-1  THEN 

6880  BEEP 

68S0  INPUT  'ENTER  THE  NAME  OF  THE  OATA  FILE"  ,D_f i leS 

6900  ASSIGN  »File  TO  D.fileS 

6910  BEEP 

6920  BEEP 

6930  INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER"  ,Md 

5940  BEEP 

6950  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED"  .Npairs 

69601  BEEP 

6970"  INPUT  "CONNECT  DATA  WITH  LINE  (  1  -Y  ,0-N  )"> "  ,  Ic  1 

6980  BEEP 

6990  PRINTER  IS  1 

7000  PRINT  USING  *4X. ""Select  a  synool:""" 

7010  PRINT  USING  "6X."-1   Star      2   Plus  sign""" 

7020  PRINT  USING  "6X.""3   Circle    4   Square""" 

7030  PRINT  USING  "6X.""5   Ronous""" 

7040  PRINT  USING  *SX,""6   Right-side-up  triangle""" 

7050  PRINT  USING  "6X,""7   Up-aide-doyn  triangle""" 

7060  INPUT  Sym 

"•070  PRINTER  IS  70S 

''080  PRINT  "PU  01" 

7090  IF  Sym-I  THEN  PRINT  "SM»" 

7100  IF  Syn-2  THEN  PRINT  "SM+" 

7110  IF  Syn-3  THEN  PRINT  "SMo" 

7120  IF  Md>l  THEN 

7130  FOR  I-l  TO  (Md-I  ) 

7140  ENTER  eFileiYa.Xa 

7150  NEXT  I 

7160  END  IF 

7170  FOR  I-l  TO  Npairs 

"•180  ENTER  CFileiYa.Xa 

7190  IF  I-l  THEN  Ql-Ya 

7200  IF  I-Npairs  THEN  02-Ya 

7210  IF  Op o-l  THEN  Ya-Ya/Xa 

7220  IF  Opo-2  THEN 

7230  O-Ya 

7240   Ya-Ya/Xa 
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I -250  »a-Q 

17260  £N0  IF 

17270  IF  Xa<Xll  THEN  Xll-Xa 

17280  IF  Xa>Xul  THEN  Xul-Xa 

17290    IF    lun-l    THEN 

17300  IF  Opo<2  THEN  Xa-Xa»1.8 

17310  IF  Opo>0  THEN  Ya-Ya».1751 

17320  IF  Opo-0  THEN  Ya-Ya«.317 

17330  IF  Opo-2  THEN  Xa-Xa».317 

17340  END  IF 

17350  X-LGT( Xa/Xnin >«Sfx 

17350  Y-LST( Ya/Yntn  )»Sfy 

17370  Kj-0 

173B0  CALL  Symb<X  ,Y  ,Syn  ,Icl  ,K j  ) 

17390  60T0  17520 

17*00  IF  Syn>3  THEN  PRINT  'S*T 

17410  IF  Syn<4  THEN  PRINT  "SR  1.4,2.4* 

17420  IF  Icl-0  THEN 

17430  PRINT  "PA"  ,X  ,Y,U  " 

17440  ELSE 

17450  PRINT  "PA"  .X  ,Y ,"PD" 

17460  ENO  IF 

17470  IF  Sym>3  THEN  PRINT  'SR  1.2,1.6' 

17480  IF  Sym-4  THEN  PRINT  "UC2 .4 ,99 ,0 ,-8 ,-4 ,0 .0  ,8  ,4  ,0i  ■ 

17490  IF  Syn-5  THEN  PRINT  * UC3 .0 . 99 . -3 ,-6 ,-3 ,S  ,3  ,6  ,3  ,-6  ■ 

17500  IF  Sym-6  THEN  PRINT  'UC0 .5. 3 .99 ,3 .-8 ,-6 ,0  .3  ,8  I ' 

17510  IF  Syn-7  THEN  PRINT  'UC0  ,-5 . 3  ,99  .-3  ,8  ,6  .0  .-3  ,-8i " 

17520  NEXT  I 

17530  PRINT  "PO" 

17540  BEEP 

17550  INPUT  "WANT  TO  LABEL?  ( 1 -Y  ,0-N  )  ■  ,  I  lab 

17560  IF  I  lab- 1  THEN 

17570  PRINT  'SP0iSP2' 

17580  BEEP 

17530  IF  Klab-0  THEN 

17600  Xlab-5 

17610  Ylab-85 

17620  INPUT  "ENTER  INITIAL  X  ,Y  LOCATIONS'  ,X lab  ,Ylab 

17630  Xtt-Xlab-5 

17640  Ytt- Ylab+8 

17650  PRINT  'SR  1.1 .5" 

17650  PRINT  'SMiPA-  ,Xtt  .Ytt  ,"LB 

17670  Ytt-Ytt-3 

17680  PRINT  -PA-  ,Xtt  .Ytt  ,"LB      Oil   Flux 

17690  IF  Syn-1  THEN  PRINT  "SM«' 

17700  IF  Syn-2  THEN  PRINT  'SM+' 

17710    IF    Syti-3   THEN   PRINT    *SMo- 

17720  Klab-1 

17730  END  IF 

17740  Kj-I 

17750  CALL  Synb(Xlab  ,Ylab  .Syn.Icl  .Kj  ) 

17760  PRINT  'SR  1,1 .5iSM" 

17770  IF  Syn<4  THEN  PRINT  *PR  2.0" 

17780  BEEP 

17790  INPUT  'ENTER  BOP'  , Bop 

17800  IF  Bop<10  THEN  PRINT  'PR  2,0iLB'iBopi 

17810  IF  Boo>9  THEN  PRINT  'PR  .  5  ,0iLB' iBop i 

17820  Ihf-0 

17830  IF  Q1>02  THEN  Ihf-1 

17840  IF  Ihf-0  THEN  PRINT  "PR  4,0iLBInc" 

17850  IF  Ihf-1  THEN  PRINT  'PR  4,0iLBDec' 

17860  PRINT  "PR  2  ,0iLB ' iD_f l la$ i ' ' 

17870  PRINT  'SP0iSP1iSR  1.5,2" 

17880  Ylab-Ylab-5 

17890  END  IF 

17900  BEEP 


Haat   File" 


Name" 
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7910  ASSI5N  9F.ie  TO  • 

7920  Kll-Xil/ 1 .2 

7938  Kul"Xul»1 .2 

7940'  SOTO  8040 

7950  END  IF 

7960  PRINT  "PU  Sh" 

7970  BEEP 

7980  INPUT  'UANT  TO  PLOT  A  POLYNOMIAL  ( 1-Y ,8-N )?"  ,6o_on 

7990  IF  5o_on-l  THEN 

3000  BEEP 

8010  PRINTER  IS  1 

8020  PRINT  USIN6  '4X, "Select  line  type:'*' 

8030  PRINT  USIN6  "6X ,**0     Solid  line*" 

8040  PRINT  USIN6  "BX  .*  *1      Dashed"" 

8050  PRINT  USIN6  *BX,"*2,,,5  Longer  line  -  dash'" 

8050  INPUT  Ion 

3070  PRINTER  IS  705 

8080  BEEP 

8090  INPUT  "SELECT  (  0-LIN , 1 -LOG( X  ,  Y  )  ) '  ,  I  log 

8100  Iprnt-1 

81 10  CALL  Poly 

8120  IF  Ifn-0  THEN 

8130  BEEP 

8140  INPUT  'ENTER  NUMBER  OF  FILES  TO  JOIN?'  ,Njoin 

8150  END  IF 

81G0  Ijoin-0 

8170  IF  Ifn<Njoln  THEN  Ijoin-1 

8180  IF  Ifn-0  OR  Ijoin-1  THEN 

8! 90  FOR  Ij-0  TO  3 

8200  Bs<  Ij  )-Bs(Ij  )+B(Ij  ) 

8210  NEXT  Ij 

8220  Ifn-Ifn+1 

8230  END  IF 

8240  IF  Njoin-Ifn  THEN 

8250  FOR  Ij-0  TO  3 

8250  B( Ij  )-Bs( Ij  >/Njoin 

8270  Bs(  Ij  )-0 

8280  NEXT  Ij 

8290  Ifn-0 

8300  ElSE 

8310  60T0  15850 

8320  END  IF 

8330  BEEP 

8340  INPUT  'ENTER  Y  LOWER  AND  UPPER  LIMITS' ,Y1 1  ,Yul 

3350  FOR  Xx-0  TO  Cx  STEP  Cx/200 

8350  Xa-Xmn*  l0*Xx 

8370  IF  Xa<Xll  OR  Xa>Xul  THEN  18640 

8380  Icn-Icn+I 

8390  Pu-0 

8400  IF  lon-l  THEN  Idf-Icn  MOO  2 

8410  IF  Ipn-2  THEN  Idf-Icn  MOD  4 

8420  IF  Ion-3  THEN  Idf-Icn  MOO  8 

8430  IF  Ipn-4  THEN  Idf-Icn  MOD  15 

8440  IF  Ipn-5  THEN  Idf-Icn  MOO  28 

8450  IF  Idf-i  THEN  Pu-I 

8460  IF  Ooo-0  THEN  Ya-FNPoly( Xa  ) 

8470  IF  Ooo-2  ANO  Ilog-0  THEN  Ya-Xa/FNPoly ( Xa  ) 

8480  IF  0oo-2  ANO  Ilog-I  THEN  Ya-FNPol y( Xa > 

8490  IF  Opo-1  THEN  Ya-FNPoly ( Xa  I 

8500  IF  Ya<Ynin  THEN  18540 

8510  IF  Ya<YU  OR  Ya>Yul  THEN  18540 

8520  IF  Iun-1  THEN 

8530  IF  0po<2  THEN  Xa-Xa«1.8 

8540  IF  Opo>0  THEN  Ya-Va«.l7S1 

8550  IF  Opo-0  THEN  Ya-Ya».317 

8560  IF  0oo-2  THEN  Xa-Xa«.317 
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16570 

loSorf 

18590 

16600 

18610 

19520 

18630 

18640 

18660 

18660 

18670 

18680 

18690 

18700 

18710 

18720 

18730 

18740 

18750 

18760 

18770 

18780 

18790 

18800 

18810 

18820 

18830 

18840 

18850 

18860 

18870 

18880 

18890 

189C0 

18910 

18920 

18930 

18940 

18950 

18960 

18970 

18980 

18990 

19000 

19010 

19020 

19030 

19040 

19050 

19060 

19070 

19080 

19090 

19100 

191  10 

19120 

19130 

19140 

19150 

19160 

19170 

19180 

19190 

19200 

19210 

19220 


END  IF 

f  ™i_6T  (.  ra/  tnin  /»Sr  y 

X-L6T(Xa/Xmn >«Sfx 

IF  Y<0  THEN  Y-0 

IF  Y>100  THEN  60T0  18640 

IF  Pu-0  THEN  PRINT  "PA*  ,X  .  Y  , -P0' 

IF  Pu-I  THEN  PRINT  "PA"  ,X  ,Y , "PIT 

NEXT  Xx 

PRINT  'PIT 

60T0  16850 

END  IF 

BEEP 

INPUT  "WANT  TO  PLOT  REILLY'S  DATA?  (  I -Y  ,0-N  ) '  .  Ir ly 

IF  Opo-0  OR  Opo-1  THEN 

XI 1-3 

Xul-20 

ENO  IF 

IF  Opo-2  THEN 

Xll-10000 

Xul-100000 

END  IF 

IF  Irly-1  THEN 

Yl  1-20 

Yul-70 

BEEP 

INPUT  "ENTER  LOWER  AND  UPPER  Y-LIMITS  FOR  PL0TTIN6"  , Yl 1  , Yul 

FOR  Xx-0  TO  Cx  STEP  Cx/200 

Xa-Xrun»10"Xx 

IF    Xa<XH    OR   Xa>Xul    THEN    18980 

Xl-L06(Xa> 

IF   Opo-0   THEN    Y1--2.5402837E-1+X1»(4.9720151-X1»2.5134787E-1  ) 

IF   Opo-1    THEN   Y1--2.5402837E-1+X1.(3.97201S1-X1«2.5I34787E-I  ) 

IF   Opo-2    THEN   Yl— 3.7073801E-H-Xl»(8.72S9l90E-!-Xl.6.8826842E-3> 

Ya-EXP( Yl  ) 

Y-LST(  Y«/Ynin.>»Sfy 

X-L6T(Xa/Xnin  )»Sfx 

Ipu-0 

IF  Y<Y11  THEN  Ipu-1 

IF  Y>Yul  THEN  60T0  18980 

IF  Ipu-0  THEN  PRINT  "PA "  ,X  ,  Y  .  "PD" 

IF  Ipu-1  THEN  PRINT  "PA-  ,X  ,Y  ,  #PU" 

NEXT  Xx 

PRINT  -PU" 

END  IF 

BEEP 

INPUT  "WANT  TO  PLOT  R0H5EN0U  CORRELATION?  ( 1 -Y  ,0-N  )*  ,  Irohs 

IF  Irohs-I  THEN 

YU-15 

Yul-80 

BEEP 

INPUT    "ENTER    Tsat    (Deo   CT.Tsat 

Csf-.O040 

BEEP 

INPUT  'ENTER  Csf  ( DEF-0 . 004  )  ■ ,Csf 

Tf-Tsat+2 

FOR   Xx-0    TO   Cx    STEP    Cx/200 

Xa-Xi<iin»10AXx 

IF    Xa<Xll    OR   Xa>Xul    THEN    19420 

Xl-L06(Xa> 

IF   Opo<2    THEN   Tf-Tsat+Xa/2 

Rho-FNRho( Tf  > 

K-FNK( Tf  ) 

Mu-FNMuUf  ) 

Cp-FNCp(Tf) 

Hfg-FNHfQ(T3at  ) 

Ni-Hu/Rho 
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I  3230  P--Cp •lu ' * 

'9140  Gmega-C  3  f  »Hf  g/Co«  '  ^  •  0'  *''  '  3  •  5  '  *'no  ■   .5   "l.«^c- 

19250  IF  Opo-0  THEN  Ya-< xa/Omega  >*3 

13268  IF  Ooo-l  THEN  Ya-( Xa/Onega  )"3/Xa 

19270  IF  Opo-2  THEN  Ya-Xa" ( 2 . /3  >/0mega 

19290  IF  Opo-2  THEN 

19290  Tfc-Tsat*Xa/Ya«.S 

19300  IF  ABSi  Tf-Tf c  )>.01  THEN 

19310  Tf-< Tf+Tfc  '..5 

19320  GOTO  19170 

19330  END  IF 

19340  END  IF 

19350  Y-L6T( Ya/Ynin  )«Sf y 

19350  X-L6T(Xa/Xmn  i*Sfx 

19370  Ipu-0 

19390  IF  Y<Y11  THEN  Ipu-1 

19390  IF  Y>Yul  THEN  19420 

19400  IF  Ipu-0  THEN  PRINT  "PA ■ , X , Y , "PO " 

19410  IF  Ipu-I  THEN  PRINT  "  PA  "  ,X  .  Y  ,  •  PU  ' 

1 9420  NEXT  Xx 

19430  PRINT  *PU' 

19440  ENO  IF 

19450  BEEP 

19460  INPUT  -UANT  TO  QUIT  (  I -Y  ,0-N  ) "  .  Igt 

19470  IF  Iqt-1  THEN  19490 

19480  60T0  16800 

19490  PRINT  "PU  PA  0,0  SP0" 

19500  SUBEND 

19510  SUB  SynP(X  ,Y  ,Syn  ,Icl  ,Kj  ) 

19520  IF  Sym>3  THEN  PRINT  "SW 

19530  IF  Syei<4  THEN  PRINT  "SR  1  .4.2.4" 

19540  Yad-0 

19550  IF  Kj-1  THEN  Yad-.8 

19560  IF  Icl-0  THEN 

19570  PRINT  "PA*  ,X  ,Y  +  Yad  ,"  " 

19580  ELSE 

19590  PRINT  *PA'  ,X ,Y*Yad  ,"PD' 

19600  END  IF 

19610  IF  Syn>3  THEN  PRINT  "SR  1.2.1.6" 

19620  IF  Sy^-4  THEN  PRINT  "UC2  .4  .99  .0  ,  -3  ,-4  ,0  .0  ,8  .4  ,0  ;  " 

19630  IF  Sym-5  THEN  PRINT  ■  UC3  ,0  ,  99  . -3  ,-S  . -3  .6  .3  ,6  .3  ,-6  ;  ' 

19640  IF  Syn-6  THEN  PRINT  ■ UC0 ,5 . 3 . 99  ,3 ,-8  ,-6  ,0  .3  .8 i " 

19650  IF  Syn-7  THEN  PRINT  "UC0  .-5 . 3  . 99  .-3  .8  .6  .0  , -3  ,-8 ; • 

19660  IF  Kj-1  THEN  PRINT  'SMiPR  0,-.8" 

19670  SUBEND 

19680  SUB  Purg 

19690  BEEP 

19700  INPUT  "ENTER  FILE  NAI1E  TO  BE  DELETED-  ,F  i  ie$ 

19710  PURGE  File* 

19720  GOTO  19690 

i 9730  SUBEND 

19740  SUB  Tdcn 

19750  COM  /Cc/  C<7),Ical 

I  9760  OIM  Enf<  1  ) 

19770  OATA  0.  10086091  ,25727. 94369. -"'57345. 8295. 7802S595. 81 

19780  DATA  -9247486589  .6 . 97688E* M  .-2 . 56  I 92E* I  3  .3 . 94078E* I  4 

19790  REAO  C< •  ) 

19800  BEEP 

19810  INPUT  "6IVE  A  NAME  FOR  FILE  TO  BE  CREATED"  ,F i  le$ 

19820  BEEP 

19830  INPUT  'SELECT  TUBE  ( 0-UH  .  I -HF  ,2-UT ) "  . 1 1 t 

19840  BEEP 

19850  INPUT  'SELECT  THERMOCOUPLE  TYPE  ( 0-NEU  ,  1 -OLD  )",  leal 

19660  IF  Itt<2  THEN  Di-.0127 

19870  CREATE  BDAT  File*. 4 

19880  ASSI6N  IF  lit  TO  FileS 
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:  9890  OUTPU'  BFiiei Itt 

19900  J-0 

19910  BEEP 

19920  INPUT  'ENTER  MONTH,  DATE  AND  TIME  ( MM: DD :HH :MM : SS > "  .DatoS 

19930  OUTPUT  709i'TD'iDat«S 

19940  OUTPUT  709i*TO" 

19950  ENTER  709iDate* 

19960  PRINTER  IS  1 

19970  PRINT 

19980  PRINT  '  Month,  date  and  time:  "iDateS 

19990  PRINT 

20000  PRINT  USIN6  "tax,""  Fns       Tin       Tev       Vu'2      Tdrop'" 

20010  IF  K-0  THEN 

20020  PRINTER  IS  701 

20030  PRINT 

20040  PRINT  ■  Month,  date  and  tine:  " iDateS 

20050  IF  Itt-0  THEN  PRINT  USING  "10X,""Tube  Type:  Uieland  Snooth""- 

200S0  IF  Itt-I  THEN  PRINT  USIN6  "10X,-"Tube  Type:  High  Flux""" 

20070  IF  Itt-2  THEN  PRINT  USIN6  "10X,"*Tube.  Type:  Turbo-B"" 

20080  PRINT 

20090  PRINT  USIN6  "I0X,""  Fna      Tin      Tev      Uw*2      Tdrop"" 

20100  PRINTER  IS  t 

201 10  K-1 

20120  END  IF 

20130  SEEP 

20140  INPUT  "ENTER  FLOUMETER  READING". Fns 

20150  OUTPUT  709i"AR  AF0  AL4  UR1" 

20160  FOR  L-0  TO  4 

20170  OUTPUT  709i  "AS  SA" 

20180  IF  L>0  AND  L<4  THEN  20250 

20190  S-0 

20200  FOR  1-0  TO  9 

20210  ENTER  709.E 

20220  S-S+E 

20230  NEXT  I 

20240  IF  L-0  THEN  Emf ( 0  >-A8S< S/ 1 0 > 

20250  IF  L-4  THEN  Enf (  I  )-ABS( S/ 1 0  ) 

20260  NEXT  L 

20270  OUTPUT  709s "AR  AF00  AL00  UR 1 " 

20280  OUTPUT  709i  "AS  SA" 

20290  Etp-0 

20300  FOR  1-0  TO  9 

20310  ENTER  709iEt 

20320  Etp-Etp+Et 

20330  NEXT  I 

20340  Etp-Etp/10 

20350  Tin-FNTvsv(Emf< 1  )  > 

203G0  Tev-FNTvsvl Enf (0)  ) 

20370  6rad-37.98S3+. l04388»Tin 

20380  Mdot-3.9657E-3*F*s»<3.61955E-3-Fms»<8.82006E-6-Fms»( 1 . 23688E-7-Fns«4 . 3  I  897 

E-10))) 

20390  Uu-Mdot/( 1000«PI»Oi*2  )»4 

20400  Tdrop-Etp»1 .E+6/( 10»6rad  ) 

20410  PRINT  USIN6  "10X ,3(DD.0D,4X  )  ,1X  ,Z.DD,4X  .MZ.40" iFns  , Tin, Tev  ,Vu*2  , Tdrop 

20420  BEEP 

20430  INPUT  "WANT  TO  ACCEPT  THIS  DATA  SET'  (  1 -Y  ,0-N  )"  ,0k 

20440  J-J+l 

20450  IF  Ok-0  THEN 

20460  J-J-1 

20470  6OT0  20130 

20480  ELSE 

20490  OUTPUT  •FlleiFFi8,E(xf(»).Etp 

20500  PRINTER  IS  701 

20510  PRINT  USIN6  "  t  0X  ,3<  OD  .00  ,4X  )  ,  1  X  ,Z  .00  ,4X  ,MZ  .  40"  iFria  .Tin  .Tev  ,^"2  .Tdrop 

20520  PRINTER  IS  I 

20530  BEEP 
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205*0  :nput  will  ~^e=e  be  anctne=  oa~*  =e"~   -r.a-N   s-  3n 

20550  '.r    3o_on-i  THEN  20'  20 

20560  ENO  IP 

20570  ASSI6N  IFila  TO  • 

20580  PRINTER  IS  701 

20590  PRINT 

20600  PRINT  USIN6  '10X.'*NOTE:  '•.ZZ,"  data  seta  are  stored  in  file  "".ISA'iJ.F 

lleS 

20610  PRINTER  IS  I 

20620  SUBEND 

20630  SUB  Uoprt 

20640  PRINTER  IS  1 

20650  BEEP 

20660  INPUT  "Enter  Uo  File  Nane"  ,F i le* 

20670  BEEP 

20680  INPUT  "Number  of  Data  Runs",Nrun 

20690  INPUT  'Do  You  Want  a  Plot  F i le'< 1 -Y ,0-N > " , Ipiot 

20700  BEEP 

20710  IF  Iplot-1  THEN 

20720  INPUT  "6ive  Plot  File  Nana"  ,P_f l leS 

20730  CREATE  BOAT  P_file$,4 

20740  ASSIGN  SPlot  TO  P_file« 

20750  ENO  IF 

20750  PRINTER  IS  701 

20770  PRINT 

20780  PRINT 

20790  PRINT  USIN6  " 1 0X  ,  ■ "    Water  Vel        Uo*"" 

20800  ASSI6N  »File  TO  FileS 

20810  IF  Iplot-I  THEN  ASSI6N  9Filel  TO  P.fileS 

20820  FOR  1-1  TO  Nrun 

20830  ENTER  ©FilesUu  ,Uo 

20840  IF  Iplot-I  THEN  OUTPUT  •FileliUw.Uo 

20850  PRINT  USING  " 1 5X  ,0 . DO  ,6X  ,MZ . 3DE" i Vy ,Uo 

20860  NEXT  I 

20870  ASSIGN  SFile  TO  • 

20880  ASSIGN  *Filel  TO  ♦ 

20890  PRINT  USIN6  "  1  0X , " "NOTE  :  ' "  ,22 , ' V  data  sets  art   stored  in  file  " * .ISA" iNru 

n. FileS 

20900  IF  Iplot-I  THEN 

20910  PRINT  USIN6  "I0X,""NOTE:  "".ZZ.--  X-Y  Pairs  are  stored  in  file  "",l5A"iNru 

n  ,P_f ileS 

20920  END  IF 

20930  PRINTER  IS  1 

20940  SUBEND 

20950  SUB  Select 

20960  COM  /Idp/  Idp 

20970  BEEP 

20980  PRINTER  IS  I 

20990  PRINT  USIN6  "4X,""Select  option:"" 

2'000  PRINT  USIN6  "6X,"  0  Taking  data  or  re-process  ing  previous  data""" 

21010  PRINT  USIN6  '6X,""  I  Plotting  data  on  Log-Log  """ 

21020  PRINT  USING  "6X,""  2  Plotting  data  on  Linear""" 

21030  PRINT  USIN6  "5X.""  3  Mane  cross-plot  coefft  file""" 

21040  PRINT  USIN6  "6X.""  4  Re-circulate  water""" 

21050  PRINT  USING  "6X,""  5  Purge""" 

21060  PRINT  USING  "SX,""  6  T-Orop  correct  ion" "  " 

21070  PRINT  USING  '6X,"  7  Print  Uo  File""" 

21080  PRINT  USIN6  "6X,""  8  Modify  X-Y  file""" 

21090  PRINT  USIN6  "6X,""  9  Move""" 

21100    PRINT    USING    "GX,""10    Contj/F  ixup  '  *  " 

21110  INPUT  Idp 

2 1120  IF  Idp-0  THEN  CALL  Main 

21  130  IF  Idp-1  THEN  CALL  Plo* 

21140  IF  Idp-2  THEN  CALL  Plin 

21150  IF  Idp-3  THEN  CALL  Coef 

21160  IF  Idp-4  THEN  CALL  Main 
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:ao*5  TnEN  CALi.  -urg 
2 '  1 83  IF  Iap»6  ~rt£U    CAL^  Tacn 
21190  IF  Idp-7  THEN  CALL  Uoprt 
21208  IF  Idp-8  THEN  CALL  Xynod 
21210  IF  Idp-9  THEN  CALL  Move 
21220  IF  Idp-10  THEN  CALL  Coup 
21230  SUBEND 
21240  SUB  Xynod 
21250  PRINTER  IS  1 
21260  BEEP 

21270  INPUT  'ENTER  FILE  NAME",File$ 
212B0  ASSI6N  CFilel  TO  Filet 
21290  BEEP 

21300  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS". Np 
21310  BEEP 

21320  INPUT  "ENTER  NEU  FILE  NAME".File2S 
21330  CREATE  BOAT  File2S,5 
21340  ASSISN  0File2  TO  File2$ 
21350  BEEP 

213B0  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS  TO  BE  DELETED"  , Ndel 
21370  IF  Ndel-0  THEN  21410 
21380  FOR  I-l  TO  Ndel 
21390  BEEP 

21400  INPUT  "ENTER  DATA  SET  NUMBER  TO  BE  DELETED"  ,Nd(  I  ) 
21410  NEXT  I 
21420  FOR  J-1  TO  Np 
21430  ENTER  iFilel iX  ,Y 
21440  FOR  I-t  TO  Ndel 
21450  IF  Nd(I  )-J  THEN  21490 
21460  NEXT  I 
21470  OUTPUT  0File2iX  ,Y 
21480  PRINT  J  ,X  ,Y 
21490  NEXT  J 
21500  PRINTER  IS  701  . 
21510  ASSI6N  CFilel  TO  • 
21520  ASSI6N  0File2  TO  • 
21530  SUBENO 
21540  SUB  Move 
215501  FILE  NAME:  MOVE 
215601 

21570  DIM  Bop (66  )  ,A(66  )  .8(66  )  ,C(66  )  .0(66  )  ,E(66  )  ,F(66  )  .6(66  >  ,H(66  )  ,J<66  )  ,K(66  )  ,L( 
66  )  ,M(66  ) 

21580  DIM  Told$(66)C14]  ,N(66)  .Ur<66>  .Ir<66> 
21590  BEEP 

21600  INPUT  "OLD  FILE  TO  MOUE"  ,D2_f ile» 
21510  ASSI6N  «File2  TO  D2_flle$ 
21620  ENTER  «File2iNrun  .Date*  .Ldtcl  .Ldtc2  ,1 1 1 
21630  FOR  I-l  TO  Nrun 
21640  ENTER  CFi  le2 t Bop( I >  ,Told$< I > 

21650  ENTER  •File2|A(  I ) ,B( I ) ,C( I ) ,D( I ) ,E( I ) ,F( I  )  ,6( I  )  ,H( I  )  ,J( I  )  ,K< I  )  ,L( I  )  ,M< I  )  ,N 
(I  ) 

21660  ENTER  «File2 tVr< I  )  ,Ir< I  ) 
21670  NEXT  I 

21680  ASSISN  0File2  TO  • 
21690  BEEP 

21700  INPUT  "SHIFT  OISK  AND  HIT  CONTINUE". Ok 
21710  BEEP 

21720  INPUT  "INPUT  BDAT  SIZE", Size 
21730  CREATE  BDAT  02_f l let  .Size 
21740  ASSI6N  »File1  TO  D2_filei 
21750  OUTPUT  IFilel tNrun  .Date*  .Ldtc 1  .Ldtc2  . 1 1 1 
21760  FOR  I-l  TO  Nrun 
21770  OUTPUT  »Filel ;Bop< I >,Told»< I > 

21780  OUTPUT  •Fllel |A(I  )  ,d(  I  )  ,C( I  )  ,D( I >  ,E( I  )  ,F( I  )  ,S(  I  )  ,H( I  ).J( I  )  ,K( I >  ,L( I  )  ,M(  I  )  , 
N(  I  ) 
21790  OUTPUT  »Filel iVr( I  )  ,Ir(I > 
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:'B0e  next  : 

2' 310  ASSIGN  9Filal  TO  • 

21820'  RENAME  'TEST*  TO  D2_fileS 

2'830  BEEP  2000,. 2 

21840  BEEP  4000,. 2 

21850  BEEP  4000,. 2 

21860  PRINT  'DATA  FILE  MOUED' 

21870  SUBEND 

21880  SUB  Comb 

21890'  FILE  NAME:  COMB 

21900' 

21910  OIM  Emf (  12  ) 

21920  BEEP 

21930  INPUT  "OLD  FILE  TO  FIXUP-  ,D2_f l le* 

21940  ASSIGN  »File2  TO  D2_fileS 

21950  D1_file»-'TEST- 

21950  CREATE  BOAT  O1_file$,20 

21970  ASSI6N  SFilel  TO  Dl.filet 

2'980  ENTER  9Fi le2 iNrun  .DateS  ,Ldtc I  ,Ldtc2  . 1 1 1 

21990  Nrunm-20 

22000  IF  K-0  THEN  OUTPUT  9F i le 1  i Nrunm , DateS ,Ldtc  1  ,Ldtc2 . 1 1 t 

22010  FOR  1-1  TO  Nrun 

22020  ENTER  •Fil«2iBop .Toldt.Emf ( • ) ,Vr .Ir 

22030  OUTPUT  9Filel tBoo  ,ToldS,Emf< » > ,Vr,Ir 

22040  NEXT  I 

22050  ASSIGN  eFile2  TO  • 

22060'  RENAME  'TEST'  TO  D2_fileS 

22070  BEEP  2000,. 2 

22080  BEEP  4000, .2 

22090  BEEP  4000, .2 

22100  BEEP 

22H0  INPUT  'UANT  TO  AOO  ANOTHER  FILE  (  1  -Y  ,0-N  P "  ,Oka 

22120  IF  Oka-I  THEN 

22130  K-i 

22140  BEEP 

22150  INPUT  "6IUE  NEW  FILE  NAME'.NfileS 

22160  ASSI6N  9File2  TO  NfileS 

22170  GOTO  21980 

22180  END  IF 

22190  ASSI6N  9File2  TO  ♦ 

22200  SUBEND 

Date  :     19  Jan  1993 

NOTE:  Program  name  :  DRPGB 
Disk  number   »  00 
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